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l. INTRODUCTION 


Tue theory and analysis of F, and parental data from a set of diallel 
crosses has been described and illustrated using the results of our 
own crosses between eight varieties of WV. rustica and published work 
on other species (Jinks and Hayman, 1953; Jinks, 1954a, 19546, 
1955; Hayman, 1954). The method, which is an extension of that 
developed by Mather (1949), permits the estimation of parameters 
for additive, dominance and environmental effects and allows the 
recognition of non-allelic interaction. The analysis has now been 
extended to cover the F, and backcross generations of a diallel set 
of crosses and is illustrated by the same diallel cross, the F, analysis 
of which was described in the earlier paper (Jinks, 19545). 

Before proceeding with the analyses of the later generations we 
will briefly restate the conclusions reached from the analysis of the 
parents and F,s. Of the two main characters followed, the time of 
flowering showed incomplete but significant dominance and a complete 
absence of non-allelic interaction for the two seasons 1951 and 1952. 
The other character, final height, appeared to show overdominance 
but this was traced to a spurious inflation of the dominance component 
by non-allelic interaction of a type comparable with the comple- 
mentary genes of classical genetics. After omitting from the analyses 
the crosses showing significant non-allelic interaction, all of which 
occurred in arrays 1, 2 and 4, re-analysis revealed only complete 
dominance. 


2. MATERIAL AND METHOD 


The 1951 and 1952 experimental layout has already been described in detail. 
In including F,s and backcrosses as well as parents and F,s the 1953 experiment 
was identical in design with that of 1952 so that we have the parents and F;s of the 
diallel cross for three consecutive seasons and the F, and backcross generations 
for the last two. 


3. THE ANALYSIS OF MEANS 
(i) Genotype-environment interaction 
In the first paper of this series (Jinks, 1954+) one method of 
detecting genotype-environment interaction was discussed, namely 
the within plot variance of non-segregating parental and F, families, 
i.e. E,. The analysis showed that for the two seasons 1951 and 1952 
I A 
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there was no significant difference between the mean E, for parents 
and that for the F,s. Furthermore, the significant heterogeneity of 
the E,s over families, i.e. both parents and Fs, was attributable to 
the varying reaction of the different genotypes to environmental 
differences. Since then estimates of the E,s for 1953 have been 
obtained and a second genotype-environment interaction, the variance 
of family means over the three seasons, has been investigated. The 
analysis of these two statistics has been presented and discussed by 
Jinks and Mather (1955) in relation to the stability in development 
of homozygotes and heterozygotes. 


(ii) Scaling tests 

The expected means of each generation following a cross between 
two inbred lines in the absence and presence of non-allelic interaction 
are given by Mather (1949) and Hayman and Mather (1955) 
respectively. In the absence of interaction the following relationships 
hold :— 

F, = }P,+2P2.+3F; 
B, = 3P,+2F, 
B, = $P,+3F, 

In the presence of non-allelic interaction, however, these equalities 
no longer hold. Although this test can be made on the three equalities 
independently as described by Mather (1949) and Mather and Vines 
(1952), a more convenient method which combines these three tests 
has been developed by Cavalli (1952) and used for the detection of 
non-allelic interaction by Jinks (19544, 1955). This test consists of 
estimating by weighted least squares the three parameters 2d, Zh 
and M (the mid-parent) taking as weights the reciprocals of the 
squared standard errors of the generation means. The expectations 
of the generations in terms of these three parameters are as follows : 


P, = M+2d (2d and 2h here refer to the balance of 
P, = M—2d genes in opposition, ¢e.g. Xd is the net balance 
F, =M+2h of the alleles with a positive contribution 
F, = M+42h to the mean (+d) minus the contribution 
B, = M+42d+432h of the alleles with a negative effect (—d)). 
B, = M—42d+42h 


The estimated parameters can be tested for consistency over 


generations by comparing the observed and expected generation 
means, in this particular case as a x? for three degrees of freedom. 


(a) Height. Our earlier analyses of height using the regression of array covariance 
on array variance suggested that we should find non-allelic interaction in arrays 1, 
2 and 4, i.e. crosses involving lines 2, 5 and 12 as one parent. They further suggest 
that while parent 2 interacts with 5 and 12 the latter do not interact with one 
another. The scaling tests fully confirm these expectations. For both the 1952 and 
1953 seasons all crosses, apart from the cross 3 x 8, exhibiting significant interaction 
on the scaling test, have one of these three lines as a parent (tables 1a and 1d). 
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One difference in detail over the two seasons is that two crosses 5 x2 and 5x4 
which show no significant interaction in 1952 (P = 0:1-0-05 and 0°7-0°5 respectively) 
do so in 1953. 

Having separated all the crosses into two groups, the interacters and non- 
interacters, we can now return to the individual scaling tests to examine the ways 
in which the three equalities given earlier have failed. For this purpose we have 


TABLE 1 


The incidence of significan' non-allelic interaction for height in the 8x8 diallel in 1952 
and 1953 and flowering time in 1953 as revealed by the joint scaling test 





Height Flowering time 
Ia 1952 | 1b 1953 I¢ 1953 

| 7 7 i oo ie ee a 
Arrays | 2 3 4 5 67 8] 2 38 5 © > Si Se 4h 6 7 6 
i ee eee ee ee ee elie ie: ile Lint 
2 Nee ene Seaer oe -F— + + =— + + a fm + 
3 -——=— + -—-—-—-+ —— - — + 
4 a Seema Po wep Se ae She ies 
5 = a i De ee 
6 -- -- - + 
y | —_ — _— 














independently pooled all interacting and non-interacting crosses over reciprocals, 
blocks and seasons for each generation (table 2). 

The expected values for the F, and backcross generation means have been 
calculated according to the three equalities given earlier. Examination of the 
observed-expected deviations in table 2 shows that the most marked difference 


TABLE 2 


The three scaling tests carried out independently for the non-interacting and interacting crosses 
pooled over reciprocals, blocks and seasons for each generation for the character height 








] | 4 
| Generation | P, P, Fr; |. + Be bos 
: | 
| Item | | 
= 
| | = 
Pooled | Observed . | 40°18 51°94 60-98 53°23 53°10 | 53°32 
4 | 53°3 
interacting Expected i sees eee ses 53°52 | 50°58 56°46 
crosses | Deviation . joss en ce —0'29 +252 | —3Z14 
| | | 
Pooled non- | Observed. | 38-76 45°01 45°93 | 43°47 | 42°54 | 44°28 
interacting | Expected ‘ | eas hes i 43°91 | 42°35 | 45°47 
crosses (| Deviation | —0'44 | +0719 | + 1°19 
| } | 











between the crosses showing significant interaction and those showing no interaction 
lies in the two backcross means (B, and B,). Thus the means of the backcross 
families showing non-allelic interaction are virtually the same irrespective of whether 
the F, is backcrossed to the smaller or the larger parent, i.e. B, and B, respectively. 

(b) Flowering time. For flowering time the scaling tests revealed no crosses 
showing non-allelic interaction for the season 1952. This includes the cross 1 x3 
which gave rise to a deviation from the F, Wr/Vr regression in a manner suggesting 
non-allelic interaction. It will be recalled that the latter interpretation of this 
deviation was discounted at the time and an alternative explanation based on the 
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pre-planting treatment in this season was put forward. Contrasting with this we 
find a number of crosses showing significant interaction in the scaling tests for the 
1953 season as shown in table rc. 

This raises a problem that we have not so far discussed, namely, that certain 
types of interaction may escape detection in the F, Wr/Vr regression analysis. 
Duplicate genes, for example, might well escape detection by this method and yet 
be picked up by the scaling tests. Further discussion of this matter will be postponed, 
however, until the results of other methods of detecting interactions have been 
described. 

One thing is certain, even at this stage, namely, that the non-allelic interactions 
for flowering time in 1953 are a distinct system from those that appear in the 
inheritance of height and, as we shall see later, for leaf length. Apart from the 
differences in behaviour in their respective F, Wr/Vr regression, they also have a 
different distribution amongst the families and give rise to different types of upsets 
in the individual scaling tests. Thus of 12 crosses showing interaction for height 


TABLE 3 


The three scaling tests carried out independently for the interacting and non-interacting crosses 
pooled over reciprocals and blocks for leaf length 








Generation P, P, F, F, B, B, 
Item 
Pooled Observed .| 18-60 | 24°42 | 23°81 2160 | 22:29 | 2298 
interacting Expected. hs ae ae 22°66 21°21 24°11 
crosses Deviation . nen ae we —1:06 | +1:08 | —1r-13 
Pooled non- (| Observed .| 18:19 23°73 22°87 22°02 20°90 23°87 
interacting Expected. aes es fre 21*92 20°53 23°30 
crosses Deviation . wie ca ei +o-10 | +0:397 | +0°57 
































in 1953 and the 14 crosses showing interaction for flowering time in the same season 
only 5 show interaction for both characters—the expectation on a random basis 
being 6. Furthermore the effect of the interaction in flowering time is in all cases 
to give F, and backcross families which flower earlier (i.e. have a lower mean) than 
one would expect from their parental and F, flowering times. 

(c) Leaf length. While the character leaf length failed to show any significant 
non-allelic interaction on the F, Wr/Vr regression analysis for the two seasons over 
which it was tested, there were some indications of interaction. For example, the 
regression coefficient was always lower than expected for no interaction (= o-7) 
but not significantly so because of the large standard error. Examination of the 
graphs of Wr against Vr shows that the large error and low regression coefficient 
are in all cases due to the point for array 1, all other points being a good fit with 
the expected slope of one. The joint scaling test for 1952 (leaf length was not scored. 
in 1953) shows that non-allelic interaction is indeed present and is mainly con- — 
centrated in array 1. Two crosses show significant interaction at the 1 per cent. 
level, namely 1X2 and 3x4, while a further six crosses, three of which occur in 
array 1, are significant at the 5 per cent. level. All crosses showing non-allelic 
interaction for leaf length also show interaction for height at a higher level of 
significance. It would seem, therefore, that we are dealing with the same group of 
interacting loci with two different levels of effect. 

The observed and expected generation means for the averages of interacting 
and non-interacting crosses are given in table 3. 

As with height, the backcross generation means for leaf length for the backcross 
to the smaller and larger parent are more alike than expected in the crosses exhibiting 
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non-allelic interaction. Unlike the height means, however, the deviation of the 


observed F, from its expectation is of the same order of magnitude as that of the 
backcrosses. 


The existence of non-allelic interaction in certain crosses raises the problem of 
rescaling. This has been investigated not only in the NV. rustica data but also in 
re-analyses of published maize diallels (Jinks, 1955). Suitable changes of scale 
such as log transformations will remove interaction from all crosses exhibiting this 
phenomenon. Unfortunately in a diallel set of crosses where only a proportion 
of the crosses show non-allelic interaction, any change of scale that successfully 
removes this interaction results in the appearance of significant disturbances in 
previously non-interacting crosses. 


(iii) Heterosis 

The relationship between the parental and F, means is the same 
for all three seasons. Thus the overall F, mean for height is significantly 
greater than that of the overall mid-parent, while the F, mean flowering 
time is significantly earlier than that of the average parent. An 
analysis of variance to test the consistency of the magnitude of the 
difference between the F, mean and the parental mean over seasons 
is given in table 4. 

TABLE 4 


The analysis of variance of the relationship of the 
overall parental and F, mean over seasons 


























Character Flowering time Height 

Item N MSS MSS 
PvE; . ‘ . . ‘ I 599°7 2151°4 
P v F, X seasons ‘ ; : 2 38°3 152°5 
Duplicate exrur ; ; P 192 3°9 5°9 
Reciprocal differences. ‘ 84 | 24°1 14°9 








Two error variances are available for testing the significance of 
the parental and F, differences and their consistency over seasons. 
These are the duplicate error derived from the SS of differences 
between identical parental and F, families in the two blocks for the 
three seasons and the SS of differences between reciprocal F, families 
in the three seasons. The latter SS is significant for both characters 
against the duplicate error ; but this significance must be interpreted 
with caution. As pointed out in previous papers, the design of the 
experiment was such that while reciprocal families were independently 
randomised prior to sowing the seed, blocks were not independently 
randomised until planting out into the field. This inadequacy of the 
design has now been remedied, but for the present results the reciprocal 
difference SS must be regarded as the better estimate of error in the 
experiment as a whole. 

For flowering time the F, heterosis is significant and consistent 
in magnitude over seasons. For height, on the other hand, there is 
a significant overall heterosis over the three seasons, but the magnitude 
of this heterosis varies significantly over seasons. 


A2 
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4. SECOND DEGREE STATISTICS 
(i) Analysis of F,s and parents for 1953 
Before proceeding to the F, and backcross analyses we will first 
consider the analysis of the 1953 F, and parental means and compare 
them with those given earlier for 1951 and 1952. 
In all essential details the 1953 results agree with the earlier ones 
(table 5). 
The character height again shows significant overdominance 


(= >1) and significant non-allelic interaction (by,;y,>1) while 


flowering time again shows incomplete but significant dominance 


D 
Analysis of array 7 for the character height, 7.e. the only array showing 
no significant interaction in the scaling test (table 15) gave only 
complete dominance. As in previous seasons, therefore, the high 
dominance ratio can be related to an inflation due to non-allelic 
interaction. 


(o< Ay <1) and no indication of non-allelic interaction (by,;y, = 1). 

















TABLE 5 
The analysis of height and flowering time for 1953 
Height 
Statistic Flowering time 
Complete Array 7 

Dp . : ‘ . 64°6635 64°6635 14.7°9608 
Mm ' . ° 231 ripe 76°4384 83:0117 
is « : ; R 158-644) ae 69°2976 
ef X “ ‘ —133°5268 —139°4684 +40°4673 
. A ; ° 3°7551 1*1821 0'5619 
uv. ‘ : ; 0°1726 ee 0°2084 
bWr/Vr_ . ‘ : 0°5626-+-0°1134 sa 0°9154-+0-0560 














Other points of similarity over seasons are the product uv for 
height, which has never varied from 0-17, and the distribution of the 
array points on the Wr/Vr graphs for both characters. As the latter 
comparison has been analysed in detail by Allard (in preparation), 
it will not be discussed further here except in reference to the F, 
and backcross Wr/Vr graphs. 


Certain statistics show differences over the three seasons. Thus 


the genetical components D, H,, H, and F are greater in magnitude 
for both characters in 1953 than in either of the previous seasons. 
For flowering time this change has apparently not influenced the 
dominance ratio, but for height the ratio is definitely higher than 
in previous seasons. On the whole the height statistics have remained 
much more constant over seasons than those for flowering time. 
Further discussion of this subject will, however, be delayed until 
least squares estimates of the statistics have been obtained. 








ae oe 


a See — 


aa 











—— - eer 
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(ii) The expectations for the F, and backcross analyses 


The expected statistics for the F, generation are of the same 
general form as those of the F, except that the contribution of h is 
halved by the one generation of inbreeding. For this reason the 
coefficients of H, and H, are } of those of the F, statistics, while the 
coefficient of F is halved, being second and first degree statistics in 
h respectively (table 6). 


TABLE 6 


F, expected means and variances 
































Original parent lines | Genotype . ; AA aa Mean of 
Male | Frequency ‘ Ua Va array 
Female Mean : ‘ da —da 
AA. ua. da. Genotype of F, . AA 4AA : Aa: faa 
Frequency - Ua? UaVa 
Mean ; j da tha uada+}vaha 
Variance . : o $dg*+tha? 4vada?+4vaha? 
aa. Ua. —da. Genotype of F, . | }AA: $Aa: faa aa 
Frequency , UaVa Ua? 
Mean : tha —da —Vada+4uaha 
Variance . : $da?+}ha? o Suada?+}uahg? 
Parental mean (Uua—Va)da Overall mean of progenies (Ua—Va)da+UaVaha 


Overall mean variance of progenies uaVada?+ $uavaha?. 


The composition of the F, variances and covariances are as 
follows :— 





| 


ee | Many independent 
F, statistic | One gene | genes 








Mean variance of arrays = UaVada?+}uavaha?— = $D+,7,H,—iF+E, 
UaVa(Ua—Va)daha 


Mean covariance of arrays = 2UaVada?—UuaVa(UaVa(Ua—Va)daha = $D—3F +~ E, 


Variance of array means = UaVada?-+}uaVaha?—ua*ve*ha?— = $D+7,H,—7;H.— 
UaVa(Ua—Va)daha 4F+1 E, 
n 
Mean family variance = UaVada?+ fuavaha” = $D+4H,+E, 











The analysis of the means of reciprocal backcross families is divided 
into two parts. The first makes use of the fact that in any cross of a 
diallel set of crosses the following relationship holds between the 
F, family mean and the means of the two reciprocal backcrosses 
provided non-allelic interaction is not present. F, = 4(B,+B,) 
(table 7). 

Thus the expected statistics for the analysis of the means of 
reciprocal backcross families of a diallel set of crosses are identical 
with those given for the analysis of F, family means. 
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The second part of the analysis of backcrosses utilises the statistics 
obtained from the differences between the means of reciprocal back- 
cross families. These terms are entirely composed of the additive 
component of variation, 1.e. d (table 7). 

The expected variance and covariance for the analysis of the 
differences between reciprocal backcross family means are as follows :— 








Ree | Many independent 
Statistic | One gene genes 
Mean variance of arrays = UaVada?” - tD+E, 
Mean covariance of arrays = +2uavada? ss +4D+* E, 
Variance of array means = UaVada?” = 1D+* E, 











The mean covariance of arrays will be positive or negative 
according to whether the analysis has been arranged such that B, 
is the backcross of the F, to the larger (or the smaller) of the two 
parents. 

One further statistic is available from the backcross generation, 
namely the mean summed reciprocal backcross family variance. For 
the one gene model the expectation of this variance is :— 


UaVadg?+UaVaha? 


which for many independent genes becomes 


+D+}H,+2E. 


(iii) Non-allelic interaction 


The contribution of non-allelic interaction to the family means 
differs over the generations under consideration here (Hayman and 
Mather, 1955). Since the contribution of non-allelic interaction is 
taken up by the four genetic parameters D, H,, H, and F, and as 
this contribution does not appear explicitly in our analyses, these 
parameters will be homogeneous over generations only in the absence 
of interactions (Mather, 1949; Mather and Vines, 1952). The 
homogeneity of least squares estimates of the four parameters over 
generations will therefore provide us with a further test of non-allelic 
interaction. 

(iv) Linkage 

Linkage, though not affecting family means in the absence of 
non-allelic interaction and hence the second degree variances and 
covariances derived from them, shows its effect in the within family 
variances of segregating generations (Mather, 1949). 

The simplest case of two linked loci, A-a and B-b, independently 
distributed in the parental lines but showing a recombination value 
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g° yearn a 2° pean . . . . 
oO . . . . 
®y®atn-+ *p(®a—*n) . . ° ® 


aouwsIeA Ajturey sassor9yoeq peooidi901 pouluNs uvoUI [[VI2ZAO 





S9SSOIDAIEG [eIOIdI9II JO VDUIIOYIP [[VIZAG 
: S9SSOIDYIVq [eIOIdI991 Jo ULIUI [[eIZAO 






































2°y’n§g + .* pent =e 2° 4PAE+ .*ptaz | . * gouvlied sAviie 
®p*n+ } *y®ng+*pta— wai tpta— tyPag+ Bptn si | uvayq jo uvsyy 
*p("a—*n) 
| uvoul [eJUsIVg 
o 2(*4—*p)F *q | 
2 aout ep) "g ali 
= q P at 7 
0 *p— *p— *p— “ut "ut-+"PE— | A | i 
eee | eee 22h eee eee BAH Aguanbe1q *tp— “tA ‘eR 
a 2(*4+*P)t 0 + 8g) 
2 (eH =e) } svi ° ’ tg amy A 
YU pi-— Pp ‘ 8g ’ 
Be “uf "ut+"PE | oo *p *p : 
eee | eee Bath eee eee "n Aguanbar,q ‘ep -®n “yy 
@a—‘a) | Ca+'at | @a—"a) | Cat'adt 
tp— | ®p . uvoyy so[euld,y 
| BA tn Aduonbergy | soyeyy 
ee | VV + addjouexy | sour quared peurs1uoO 








Samuviiva puv suvaut paisadxa SSOLIYWY 


L alaVviL 











10 J. L. JINKS 


p(=1—q), makes the following contribution to the mean variance 
of F, and summed reciprocal backcross families from a set of diallel 
crosses. Mean family variance of F, 


= UaVada?+UpVndp?+3uaVaha?+dupvphp? 
+ 2UaVallpVp (I —2p) *hahp 


Mean summed reciprocal backcross variance 


= UgVada?+Upvndp?+UaVaha?-+-uUpvphp? 
+4UgVallpVp(1 —2p)hahp 


If we define two linkage parameters such that 


H, = 162uaVapVvy(1 —2p) *hahp 
and H, = 162uavVaupvy(1 —2p)hahp 


Then for many genes, some of which show linkage, the expectations 
for the two mean family variances are 
F, +D+$H,+}3H;+E, 
backcross }D+}H,+}H,+2E, 


Both linkage parameters can take sign. Thus if the linked genes 
show reinforcing dominance, 1.e. the dominance deviations have the 
same sign, then H, and H, will be positive. If, on the other hand, 
dominance is in opposition, 7.e. the dominance deviations have opposite 
signs, then H, and Hy, will be negative. Linkage between genes 
showing reinforcing dominance will, therefore, lead to an inflation 
of the within family variance, while linkage between genes showing 
opposing dominance will lead to a deflation. 

Although complete specification of simultaneous linkage and non- 
allelic interaction is now possible (Hayman and Mather, 1955) the 
limited number of statistics available in the present analyses is not 
sufficient for their estimation. The question now arises as to how far 
tests of heterogeneity of the specified components D, H,, H,, F and 
E,, over statistics, can distinguish between the two sources of disturb- 
ance. In so far as components for non-allelic interaction appear 
in all our statistics (except for the mean variance of F, families), 
it follows that non-allelic interaction can lead to heterogeneity of 
components whether these be derived from an analysis of family 
means or family variances of segregating generations of a diallel set 
of crosses. Linkage, on the other hand, effects only family variances 
and so can lead to heterogeneity only if this type of statistic is included 
in the estimation of components. Any heterogeneity of components 
over the statistics excluding family variances must, therefore, be 
ascribable to non-allelic interaction. Should such heterogeneity exist, 
then any further heterogeneity introduced by including the family 
variances in the estimations can obviously no longer be unambiguously 
ascribed to linkage. Our linkage test, therefore, is only valid where 
non-allelic interaction plays no significant part in the inheritance 
of the character under consideration. 
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(v) Non-random distribution of alleles in the parent lines 


One further source of disturbance may lead to heterogeneity of 
the components of variation as estimated from the analysis of family 
means and those obtained from the mean variance of families of 
segregating generations. Such disturbances may arise from non- 
random distribution of alleles in the original parental lines. Since 
these do not affect the homogeneity of the components of variation 
obtained from the diallel analyses of family means, they will not be 
confounded with disturbances arising from non-allelic interactions. 
They will, however, be confounded with linkage. A full account of 
the possibilities and magnitudes of such confounding is outside the 
scope of this paper and only pertinent results will be mentioned here. 

Association of alleles combined with reinforcing dominance will 
inflate the components of variation obtained from family means but 
will not affect those obtained from the mean family variances. It 
will thus mimic linkage between loci exhibiting opposing dominance. 
The alternative situation, i.e. dispersion of alleles and opposing. 
dominance, in so far as it will deflate the components of variation 
from family means while not affecting those from mean family 
variances, will mimic linkage between loci exhibiting reinforcing 
dominance. Combinations of association and opposing dominance 
or dispersion and reinforcing dominance will also mimic linkage, 
but the magnitude of the effect will be much less than that arising 
from the previous two combinations, while the linkage phase it mimics 
will depend on the relative magnitude of the deviation from random 
distribution and the dominance ratio. 

Two considerations make the confounding between correlated gene 
distribution and linkage less serious than it might appear at first. 
Firstly, non-random distribution of alleles can be detected before the 
estimation of the effects due to linkage is undertaken. The detection 
depends once again on the fact that the regression of Wr on Vr is 
a line of slope one only if our hypothesis of independence of the genes 
is true. It can therefore be used to detect not only non-allelic inter- 
actions but also non-random distribution of the alleles. Thus associa- 
tion gives a curve which is convex upwards while dispersion gives a 
curve which is convex downwards. Although a regular curvature 
of the line is unlikely to be mistaken for non-allelic interaction we 
have at least two methods of discrimination should such difficulty 
arise, namely, the scaling test and the test of homogeneity of the 
components of variation over statistics derived from family means, 
both of which detect only disturbances arising from non-allelic inter- 
actions. In the one case where difficulty may arise in interpreting 
the Wr/Vr regressions, i.e. where both sources of disturbance are 
present, the estimation of linkage is in any case invalidated by the 
presence of non-allelic interactions. 

Secondly, in some cases, for example where the inbred lines used 
in the diallel crosses are derived from a common mating pool, any 








12 J. L. JINKS 


non-random distribution of alleles will probably be due to the linked 
genes themselves. In such cases confounding between the two, while 
it might disturb the magnitude of the linkage effect, will not lead 
to the drawing of wrong conclusions regarding the presence or absence 
of linkage. 

In our own data, where a valid estimate of linkage effects is 
confined to the 1952 flowering time results, there is no suggestion of a 
downward convexity of the Wr/Vr regressions, i.e. dispersion of alleles, 
which would give rise to the significant reinforcing linkage (see later). 


(vi) The regression of array covariance on array variance 

The relationship between array covariance and array variance 
described for the F, analyses also holds, within the same limitation of 
independence of the genes, in the F, and averaged reciprocal backcross 
analyses, t.e. the regression of array covariance (Wr) on array variance 
(Vr) has a slope of one. There is one difference, however, in that 
the point of interaction of the regression line with the Wr axis is no 
longer Wr—Vr = }(D—H,) but }(D—}H,). 


TABLE 8 
The regression coefficients of Wr on Vr combined over blocks 
and seasons for the three generations 




















Character Height Flowering time 
| a : ware 
Generation | Complete data | 2, 4 and 7 diallel Complete data 
ecticonctescamat rite a ACES TAT ae eEy 
F, . " : 0:598--0°173 | 1°057--0°081 0943-0020 
ae ; : 0:802+-0°054 1°017-+0:062 1:083-++-0°155 
Backcross : 0:677-L0-219 | 1:004-+-0°088 0'999-E0°0590 








Interaction between non-allelic genes may lead to deviation from 
the expected slope in all three generations, although in our own data 
the deviation is less marked in the F,s and backcrosses. The regression 
coefficients for the three generations obtained from a combined 
analysis over blocks and seasons are given in table 8. Also included 
for comparison are the regression coefficients of a 3 x 3 diallel extracted 
from the height data which contains no cross showing non-allelic inter- 


action (table 14, 15). This small diallel consists of the three parents 


of arrays 2, 4 and 7 and their F,, F, and backcross progeny means. 

For flowering time and the 3 x 3 diallel extracted from the height 
data the regression coefficients do not differ significantly from one, 
so that there is no suggestion of non-allelic interaction in any of the 
three generations for either season. The regression coefficients for 
the complete height data, on the other hand, always give a regression 
coefficient less than one, but it is only significantly so for the F, 
generation. 
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In the absence of interaction, the Wr/Vr graphs for the F, and 
backcross generations should be identical, but in the presence of 
interaction this is not necessarily true. For both characters we find 
that the two regressions do not differ significantly either in slope or 
mean. While the joint regression for flowering time (fig. 1) does 
not differ significantly from a slope of one, the joint regression for 
height is significantly less than one (b = 0:722+0°130), indicating 
significant non-allelic interaction for height in these two generations. 

A further relationship exists between the regressions of Wr on Vr 
for the F, and F, (or backcross) generations. If for each array we 
draw a line through the Wr/Vr coordinates for its F, and F, array 


50 


Wr 
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Fic. 1.—The regression of array covariance (Wr) on array variance (Vr) for the F, (full 
circles) and backcross (crosses) generations for flowering time. The array points of 
both generations fall on the same straight line of unit slope. This indicates absence 
of non-allelic interactions but presence of dominance. The 8 arrays fall in the same 
order along this line in both generations. This order reflects the varying proportions 
of dominant and recessive alleles in the common parents of the arrays. Those with 
most dominant alleles have low variances and covariances, é¢.g. 5, 8 and 7, while those 
with most recessive alleles have high variances and covariances, e.g. 1, 3 and 6. 


points, they will all converge on a common point of intersection 
whose own coordinates vary characteristically with the degree of 
dominance. The Wr coordinate for the point of intersection is 
independent of the dominance relations as it contains only terms in d, 
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being 4D = S2uvd*. The Vr coordinate, however, does vary with 
the degree of dominance and for a one gene model is given by the 
expression 


VaVadaha — uaVaha? 
UaVadaha 





UaVa(da—ha) 2 -- . QUaVadaha. 
This expression has a maximal value of }D when h = 0. 

The point of intersection for other generations, e.g. F, and Fs, 
can be obtained by halving the contribution of dominance in the 
above expression to allow for the additional generation of inbreeding. 


Thus for the nt» and (n-+1)t® generation of inbreeding the coordinates 
of intersection will be 


Wr = 2UaVada2 and Vr = uaVa (a.— ha) 
n 


I I 2 
- UaVadaha—#uavVa (: ha) 
n n 





2 
. — UgVadaha 


oe 


I 
— UaVadaha 
n 


where n is large the Vr coordinate approaches the value of 


4D => DUaVade?. 


Therefore with a low initial degree of dominance or a large number 
of generations of selfing the point of intersection approaches the 
coordinates for the Wr/Vr array points of an F, diallel showing no 
dominance. 

These relationships break down in the presence of non-allelic 
interaction in such a way as to allow us to detect the interacting 
arrays. This may be illustrated by the following two gene models. 
The first model, which shows no interaction, consists of four parental 
lines AABB, AAbb, aaBB and aabb such that d, = dp = hg = hp = 2. 
The second model is the same except that the allele A in the presence 
of B contributes 4 units instead of 2 to the mean. In both these models 
the points on the Wr/Vr graph for arrays AAbb and aaBB are identical 
leaving us with only three distinct points on the graphs. These are 
given in figs. 2 and 3. For the first model the lines shown through 
the F, and F, points of each array intersect at the point expected for 
a system showing complete dominance, namely Wr = 4D, Vr = 3D. 
In the second model, the line drawn through the array points for 
arrays AAbb and aaBB passes below the expected point of intersection 
while the other two arrays pass slightly above it but intersect very 
close to it. 

Illustrations of both types of model can be found in the NV. rustica 
data. For this purpose blocks and seasons have been combined, 
and the F, and backcross Wr/Vr points pooled to give the graphs of 
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Fic. 2.—The above graph shows the F, (full circles) and F, (crosses) points of the Wr/Vr 
regression for a two gene model showing complete dominance. The lines drawn 
through the F, and F, array points for each array intersect at a common point. The 
Wr coordinate of this point is 4D and is independent of the degree of dominance. 
The Vr coordinate on the other hand increases with decreasing dominance to a 
maximum value of }D. For the complete dominance model Vr = 4D. 


20 








“10 ° 10 20 
Ver 


Fic. 3.—The above graph shows the F, (full circles) and F, (crosses) points of the Wr/Vr 
regression for a two gene model showing a complementary type of non-allelic inter- 
action. There is no common point of intersection when lines are drawn through the 
array points for the F, and F, generations for each array. This absence of a common 
point of intersection provides a means of detecting non-allelic interaction. The 
Wr/Vr regressions, of course, have a slope of less than 1. 
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figs. 4 and 5. For the complete height analysis there is no common 
point of intersection, as would be expected from the high incidence 
of non-allelic interaction. The lines drawn through the array points 
for the F, and combined F, and backcross generations fall roughly 
into three groups (fig. 4). The first group consists solely of array 7, 
which, it will be recalled, has never given any indication of the presence 
of non-allelic interaction in any of the tests so far applied. This line 
passes through the coordinates Wr = $D, Vr = 4D, which are those 
for a system showing complete dominance. The second group consists 
of arrays 2, 4 and 8 which intersect at a point whose coordinates 
are approximately Wr = 4D, Vr = }D. These three arrays have 
certain interaction properties in common. ‘Thus they all show 
significant non-allelic interaction in their crosses with the common 
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Fic. 4.—Illustrates the application of the method of detecting non-allelic interaction 
shown in figs. 2 and 3 for our height data. The F, array points are indicated by 
full circles and the F, and backcross generations have been combined to give the array 
points indicated by the crosses. The presence of non-allelic interactions leads to the 
absence of a common point of intersection. Only array 7, which contains no inter- 
acting crosses, intersects the Wr coordinate of the expected point of intersection 
(indicated by the dotted line), The remaining arrays fall roughly into two groups 
consisting of arrays 2, 4 and 8 and 1, 3, 5 and 6, a grouping that is confirmed by the 
scaling tests (tables 1a and 1b). The two Wr/Vr regression coefficients are less than 1. 


parents of arrays 1, 3 and 6 (table 1a, 14). The remaining four 
arrays I, 3, 5 and 6, each of which shows interaction with two or 
more of the common parents of arrays 2, 4 and 8, all lie above the 
expected point of intersection and only two of them, 1 and 3, intersect 
with one another where Wr = D and Vr = #D. 

If we now analyse the small 3 x 3 diallel, involving the interaction- 
free crosses of arrays 2, 4 and 7, in the same way we obtain a picture 
similar to our first model (fig. 5). There is only one point of inter- 
section and its coordinates are approximately Wr = $D and Vr = 35D. 
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After allowing for the environmental components of the Wr’s and 
Vr’s, this situation is compatible with a system showing just under 
complete dominance and absence of non-allelic interaction. 

When we apply the same test to the 1953 flowering time data 
we obtain unambiguous evidence of widespread non-allelic interaction 
as suggested by the earlier scaling tests. If we use the array points 
for the F, and joint F, and backcross Wr/Vr regressions we find 
that the lines drawn through only two arrays, 4 and 7, pass anywhere 
near the point of intersection expected from a genetical system showing 
partial dominance. This is in agreement with the scaling test which 
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Fic. 5.—Illustrates the application of the method of detecting non-allelic interaction 
shown in figs. 2 and 3 to a portion of the height data known to contain no crosses 
showing non-allelic interaction. This is identical with the model of fig. 2. There 
is only one point of intersection which falls on the expected Wr coordinate, i.e. Wr = $D 
(indicated by dotted line), while the Vr coordinate is approximately that expected 
of a system showing almost complete dominance, The F, Wr/Vr array points are 
again indicated by full circles and the combined F, and backcross points by crosses, 
The two Wr/Vr regressions have a slope of 1. 
suggests that these two arrays exhibit only a low incidence of crosses 
showing non-allelic interaction. The lines for the remaining arrays, 
while not showing a common point of intersection, all converge on a 
small area whose centre is approximately at the point Wr = 4D, 
Vr = }D. In the 1952 analysis we find no serious evidence of inter- 
action on this test, which is again in complete agreement with the 
scaling test and also the Wr/Vr regression analyses. Although three 
of the lines drawn through the array points, namely those for arrays 5, 
B 
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8 and 6, do not pass through the expected point of intersection, they 
do pass close enough for any deviation to be regarded as the con- 
sequence of error variation. 

While this method of detecting non-allelic interaction is not 
independent of the Wr/Vr regression method it does allow us to 
detect and in some cases classify the interacting arrays in a manner 
approaching the accuracy of the scaling test. 


5. ESTIMATION OF THE COMPONENTS OF VARIATION 

The process of estimation has been dealt with in detail by Mather 
(1949) and Mather and Vines (1952). In these analyses we depart 
from this method in only one respect, namely the treatment of the 
environmental components of variation. Of these only E, is 
independent of the design of the experiment, the rest, which are 
derivatives of E,, being dependent on the number of parental lines 
employed. Because of the obvious advantages of a general least 
squares solution applicable to all diallel crosses, it was decided to 
correct the observed statistics for all environmental components 
dependent on “n” prior to the estimation of the remaining para- 
meters. 

E,, which is included in the least squares estimation, and is only 
used in connection with the family variances of the’ segregating 
generations, is obtained as the mean within plot variance of the non- 
segregating parental and Fy, generations as described by Mather 
(1949). The E, component, however, is derived independently for 
each generation, as the mean variance of duplicate or reciprocal 
plot means, and is used only in connection with the generation from 
which it is derived. This is necessary because the variances and 
covariances of the diallel analysis of segregating generations contain 
both an environmental and a sampling component. By using an 
E, derived from the same generation, both these components are 
accommodated simultaneously, while an E, predicted from parents 
and F,s, as described above for the E,s, would depend only on the 
environmental portion. In many cases, the sampling error contributes 
the larger proportion of the non-inherited variation. This may be 
illustrated by reference to the E, components for the height analysis. 
For this purpose all E,s have been expressed on a single plot basis 
and averaged over blocks and seasons (table g). 


The first column of table 9 gives the E,s calculated from the mean - 


variance of duplicate and reciprocal plot means for each generation. 
The second column contains the expected E,s for the F, and backcross 
generations assuming no sampling error as $Py.+4Fyg2. The excess 
of the observed E,s over this predicted E,, which is ascribable to 


sampling error, is contained in the third column. For comparison - 


are included the expected sampling error calculated as : (mean plot 
n 


variance) where n is the number of plants per plot, 2.e. five. 
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Having corrected the observed statistics for all these non-inherited 
components of variation except E, we can now proceed with the 
least squares estimation. 








TABLE 9 
The environmental and sampling components of variation for height 
Estimated Estimated Expected 
Generation Observed E, environmental sampling sampling 
component component component 
| ar : : 4°0538 4°0538 0*0000 
F, . : 8+1960 8-1960 0*0000 aa 
ae ‘ ‘ 11°1847 6+1249 5°0598 7°9403 
B, and B, ‘ 17°5213 6+1249 11*3964 12°5843 























(i) Inclusive analysis 
Without making allowances for the possible effects of linkage, 
i.e. the inclusive analysis of Mather, five components of variation are 
involved, D, H,, H,, F and Ey. The sixteen basic equations available 


TABLE 10 
Inclusive ¢ matrix 


0°643932 CH,D 0°615684 CH,D —o0-257111 CFD 1*558234 CE,D —o0-144625 
0°615684 CH,H, 19°840244 CH,H, 10-411285 CFH, 9:224883 CE,H, —2-143653 
CDH, —0-257111 CH,H, 10-411285 CH,H, 26°573565 CFH, —2°054586 CE,H, —1:052370 

1°558234 CH,F 9°224883 CH,F —2-054586 CFF 12:021402 CE,F —1°155705 
—0'144625 CH,E, —2°143653 CH,E, —1°052370 CFE, —1°155705 CE,E,  0°408042 


for their estimation may be combined to give five equations yielding 
least squares estimates of the five components as described by Mather 
(1949.) The solution of these five sets of equations leads to a matrix 
of multipliers as shown in table 10. 


TABLE 11 
Exclusive ¢ matrix 


068508 CH,D 1*14917 CH,D 0-00000 CFD 1°85635 CE,D 000000 
114917 CH,H, 27°61940 CH,H, 14:22222 CFH, 13:43646 CE,H, 0-00000 
o-00000 CH,H, 14:22222 CH,H, 2844444 CFH, o-o0000 CE,H, 0°00000 
185635 CH,F 13°43646 CH,F 000000 CFF 14°32044 CE,F 0*00000 
0:00000 CH,E, 0-00000 CH,E, 000000 CFE, o0-00000 6=6(CE,E, 1*00000 
—2°51934 CH,H; —29°:91774 CH.H; —14:22222 CFH,; —17:14917 CE,H; —8-00000 
—1°83425 CH,H, —28-76857 CH,H, —14:22222 CFH, —15:29282 CE,H, —8-00000 


CH;D —2-°51934 CH,D —1-83425 
CH;H, —29°91774 CH,H, —28-76857 
CH;H, —14°22222 CH,H, —14:22222 
CH;F —17:14917 CH,F —15:29282 
CH;E,  —8-00000 CH,E, —8-00000 
CH;H, 162:95641 CH,H,; 96-43708 
CH;H, 96°43708 CH,H, 110°60283 


(ii) Exclusive analysis 
The analysis allowing for the presence of linkage, i.e. the exclusive 
analysis of Mather, involves two further components of variation, 
H, and H,, making seven in all. We have the same sixteen equations 
for their estimation and we can again obtain a matrix of multipliers 
(table 11). As three of the components H;, H, and E, are confined 
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to three equations, the latter must be a perfect fit with expectations. 
This provides an alternative method of exclusive analysis, whereby 
least square estimates of the four components D, H,, H, and F are 
obtained from the remaining thirteen equations. The other three 
components are then obtained by direct estimation from the observed 
statistics using the least squares estimates of D and H,. Both methods 
give identical results, but the seven by seven matrix of multipliers 
has the added advantage that the standard errors of all components 
can be obtained directly from it. 

From the observed statistics we can now find the least squares 
estimates of the components of variation in both the inclusive and 
exclusive analyses for each of the two blocks within the two seasons. 

















TABLE 12 
Analysis of variance of the height and flowering time data 
Height | Flowering time 
Item _ — a 
| | 
N MSS N MSS 
ee rn See ee ree 2S Sn ere Comer 
Linkage. 2 ‘ i ae Pe | 2 |. 134°66 
Residual interaction . ; II 150°39 9 28-68 
| | 
sista a Aer ee | we ea = = 
‘ Between seasons— 
Linkage. : ‘ ag sce 2 170°29 
£ Residual interaction . 11 31°41 9 88-94 
8 Components 5 261°77 5 5953°10 
Ee Se: ee Me! eer eee 
3 
2 Within seasons 
Ho Linkage . : : nen ‘se 4 14°96 
Residual interaction . 22 g'I4 18 | 4°59 
\ Components. ; 10 IT1*g2 10 43°36 











By substituting these estimates in the sixteen basic equations we can 
arrive at least squares estimates of the expected values of the observed 
statistics. The analysis of variance of the sum of squares of deviations 
of observed from expected and the consistency of these deviations 
over blocks and seasons have been dealt with in detail by Mather 
and Mather and Vines, so that only the final analyses are given here 
(table 12). 


As pointed out earlier, significant non-allelic interaction invalidates — 


our test for linkage, and therefore the linkage items are not given for 
the height analysis although they have been estimated and found 
to be non-significant. 


(ili) Consistency of components over seasons and blocks 


The above analysis reveals a marked difference in the stability 
of the components of variation over blocks and seasons for the two 
characters height and flowering time. Thus the components of 
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variation for height are not significantly more variable between 
seasons than they are between blocks in the same season. The com- 
ponents of variation for flowering time, on the other hand, are 
considerably more variable over seasons than they are within seasons. 
This one might expect if the change in the components in the former 
case had been brought about in response to differences in soil fertility 
presented by different blocks, while in the latter case they were induced 
by differences in the weather conditions in the two seasons. This 
explanation is not only adequate to explain our present analyses but 
has the added advantage of falling into line with the findings of 
Mather and Vines, in respect of the character height, and with current 
views on the physiology of the two characters we are investigating 
here. 
TABLE 13 


The components of variation of the exclusive analysis for 
flowering time 1952 and inclusive analysis for 1953 























Season 
Component —— 
| 1952 | 1953 
opens 

D. ; ; ; 20°0288-++ 3:1603 | 131°9830-+ 3°9340 
H, : : 4°4254+20°0658 75'077721 8369 
H, . . 2°1550-+ 20°3633 65°5492-+25'2722 
3 : F ; 5°2862-+ 14-4487 | 31 0838-++ 16-9979 
E, . - . ‘ 13°5321+ 3°8181 22°5840-+ 3°1316 
H; : : : 110°2602-++-48°7401 | ee 
H, ‘ : - | 132°8985-+-40°1544 





The above analysis of variance does not allow us to partition the 
heterogeneity of components over seasons between the additive and 
dominance components of variation. We can, however, arrive at a 
partial solution of this problem by comparing the least square estimates 
of the various components over the two seasons. The components 
for flowering time obtained from the exclusive analysis in 1952 and 
the inclusive analysis in 1953 are given in table 13. 

Since the analyses of the separate seasons reveals no significant 
disturbances after allowing for linkage in 1952, the standard errors 
of the components for this season have been calculated from the 
sum of squares of block differences of the observed statistics and so 
are based on sixteen degrees of freedom. 

In 1953 there are other significant disturbances which will be 
discussed later and for this reason the standard errors of the components 
include the significant sum of squares for heterogeneity over statistics 
within blocks (i.e. residual interaction) and are therefore based on a 
total of twenty-seven degrees of freedom. For the same reason the 
linkage components are not given as their estimation is invalided. 

Of the 1952 components of variation for flowering time H,, H, and 

B2 
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F are not significant, while only F is non-significant in 1953. That 
dominance plays a part in the inheritance of flowering time in 1952 
is, however, beyond dispute. Apart from the significant linkage 
components which are both terms in Zh, we have the significant 
Wr/Vr regressions for all three generations, and the significant heterosis 
for early flowering. But what is of more importance to our present 
discussion is that all components of variation apart from F are 
significantly different in the two seasons. Furthermore the dominance 
components H, and H, have changed by a greater factor than any 
of the other components of variation. This greater change, while 
large and consistent, is not significant. While, therefore, there is a 
suggestion that the dominance components are somewhat more 
variable over seasons than the additive components, in the present 
analyses they are not significantly so. 


TABLE 14 
The components of variation of the inclusive analysis of height 























Season 
Component | 
1952 1953 ' | 

block I block II combined blocks 

Dp : - | 37°1331- 3°7209|  44°4929-+ 5°2323|  66-0078-+ 5:9605|  55°2504-++ 5°0740 
mys } 76°7428--20°6538 | 217°4862+29'0433 | 228-4762+-33°0852 | 222:9812-++-28-1648 
H, . ‘ 61°7516-+23'9030 | 155°6763+33°6122 | 181+2954+ 38-2900] 168+4859-+-32°5956 
F ; . |—60-9226+ 16-0770 |—149°3846+-22°6074 |— 167-2167-+-25°7536 |—158-3007-+21°9296| 
E « : 7°2307-+ 2'9620 5*2820-+ 4°1651 10°7260-+ 4°7447 8-0040-+ 4°0391| 
} 





We find essentially the same situation for the components of 
variation of height. Owing to the presence of significant non-allelic 
interaction the standard errors of the components are based on the 
pooled sums of squares of block differences and residual interaction 
for a total of twenty-seven degrees of freedom. The components 
themselves, of course, are variously biased by the non-specified effects 
of the non-allelic interaction. They do serve, however, to illustrate 
the present discussion. 

Apart from E, all other components of variation are significantly 
higher in 1953 than in 1952. This increase in magnitude is again 
more marked in the case of the dominance components H, and H, 
than in the additive component D, but again there is no significant 
difference in this respect between the two types of components. Since 
the components show almost as large a difference between blocks 
within seasons as they do between seasons, the components for 1953 
have been separated into the two blocks. The standard errors of 
these components have been derived from the heterogeneity of com- 
ponents over the observed statistics within each block, and are hence 
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based on eleven degrees of freedom. Although these standard errors 
are inflated by disturbances arising from significant non-allelic inter- 
action, they are the only ones available and will therefore have to 
serve for our present purpose. The only component that shows a 
significant difference in the two blocks is the additive component D. 
We have therefore examples where the dominance components are 
apparently less stable over environmental fluctuations and others 
where the additive component seems the less stable. In conclusion 
we can only say that on the whole the two components of variation, 
the additive and the dominance, are equally susceptible to changes 
in the environment, the latter in the case of height being predominantly 
relatable to the soil conditions. 


TABLE 15 
The analysis of variance of flowering time in 1952 and 1953 




















1952 1953 
Item 

N MSS N MSS 
Linkage. , - ; 2 168-71 2 136-21 
Residual interactio ‘ . 9 15°87 9 101°76 
= a 
‘g | Between blocks— 
©! Linkage. . 2 13°35 2 16-58 
© Residual interaction 9 3°88 9 5°31 
= | Components 5 80-98 5 5°63 
eo 

















(iv) The linkage components 


The analysis of the two seasons separately shows that linkage plays 
a significant role in the inheritance of flowering time in both seasons 
(table 15). 

Not only is this linkage apparently significantly variable over 
seasons (tables 12 and 13) but in 1953 there is also significant residual 
interaction which, while confirming the findings of the scaling tests, etc., 
invalidates both the test for linkage and the least squares estimates 
of the linkage components. Further discussion of linkage must therefore 
be confined to the 1952 analyses where no further interactive disturb- 
ances have been revealed by any of the available tests. 

In 1952 the linkage components H, and H, are both positive and 
significantly different from zero. This suggests a preponderance of 
linkage between loci exhibiting reinforcing dominance. Since we only 
estimate the balance of the two types of dominance relations of the 
linked loci we cannot exclude the possibility that a proportion of the 
linked loci controlling flowering time shows opposing dominance 
relations. 
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From the ratio of the least squares estimates of the two linkage 
components H, and H, we can obtain an estimate of the linkage 
relations of the linked loci. 


Thus H, 162UaVaupVo(I —2p)*hahp 


H, 162 uaVauyVvn(1 —2p)hahp 





which when ha = hyp=h, and ua = up = u,, becomes 


en 


3 - 
ss = 1—2p 
4 


ae 


where p is the mean recombination frequency of the linked genes. 
If these conditions are not fulfilled then the mean recombination 
frequency will be weighted in favour of linked genes with larger 


TABLE 16 
The fall ratios 











Fall ratios Our comparable Observed fall ratios for 
(Mather, 1949) | estimate 1952. Flowering time 
tow | Ys 
\a*) | as C ‘0806 
1 Dr, | H,+-2H, 0:0197+0:089 
r (h2) | _H, eal 
< Hr, | H,+2H, 0°0164-++0°0777 
Dr,—Dr, | 2H,—2H, -1696-++0° 
3 Dr, H, 12H, 0°1676-+0°4702 














dominance effects and equal allele frequencies. For the 1952 flowering 
time data we find that p = 0-0852-++0-:1277. Owing to the relatively 
large standard error of this estimate it tells us very little except that 
there is significant linkage, i.e. p is significantly smaller than 0-50. 

With preponderant reinforcing linkage such as we have found 
in the 1952 flowering time analyses we can use the relative magnitude 
of H, and the two linkage components H, and H, to obtain some 
idea of the number of linked genes (Mather, 1949). The difference 
between these three components depends on the relative magnitude 
of hahp etc., on the recombination frequencies and the number of 
genes involved. ‘The estimation of the number of genes assumes 
equality of the hahp etc. and equal spacing along the genetical 
chromosome. Failure of either or both assumptions leads to an 
underestimation of the number of genes as does also the presence of 
a proportion of genes showing opposition linkage. The process of 
estimation depends on the observation that the maximum fall ratio 
is characteristic of the number of linked genes within the limitations 
of the above assumptions. Three fall ratios described by Mather can 
be adapted to our present components (table 16). 
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In equating our fall ratios to those of Mather we are making one 
further assumption, namely that for each pair of linked genes 
4UaVa+4UpVp = 32UaVallpVp- Non-equality will lead to a further 
reduction in the apparent number of linked genes, as calculated by 
all the above fall ratios. 

All three observed fall ratios are smaller than their own standard 
errors. In fact the third fall ratio in table 15 gives no evidence of 
linkage at all. The other two, in so far as they are both significantly 
smaller than one, show that there is at least significant linkage in 
the flowering time data. Furthermore they are both significantly 
smaller than the maximum possible fall ratios for four linked factors 
which are 0-40 and 0-54 respectively (P = 0-001). 


(v) Disturbances due to interaction 
As pointed out earlier in this account, interaction, though not 
appearing explicitly as components of variation, may nevertheless be 
detected by their effects in causing heterogeneity of the specified 


TABLE 17 
Exclusive matrix for single array analysis 
CDD 0*716050 CH,D 1°185185 CF; —1°925926 
CDH, 1°185185 CH,H, 28°444444 CF;H, —14°222222 
CDF; —1*925926 CH, fF; — 14°222222 CF;F; 15*111043 


components. In the case of height there is significant residual inter- 
action in both seasons. The joint analysis of 1952 and 1953 (table 12) 
confirms the significant residual interaction and, while showing 
significant heterogeneity of interaction over seasons, this heterogeneity 
is significantly smaller than the seasonal heterogeneity of the specified 
components of variation. In fact the residual interaction items are 
more stable over seasons and blocks within seasons than any other 
component of variation. 

A constant feature of the analyses of these experiments by the 
Wr/Vr regression method has been that, following the omission of 
all crosses showing non-allelic interaction from the data, we have 
never found a dominance ratio significantly greater than one. We 
can now test this finding on a joint analysis over generations and 
seasons. Two fragments of the height data have been consistently 
free from interaction on all the tests made so far, namely array 7 and 
the 3 <3 diallel between the common parents of arrays 2, 4 and 7. 
The latter can be analysed by the method already employed for the 
complete diallel, but the analysis of the single array requires a new 
matrix of multipliers. The equations appropriate for the analysis 
of single arrays have been given elsewhere (Jinks, 1955) and they 
provide nine statistics for the estimation of three components of 
variation, D, H, and F, (Fr = 8Zuvdh). The appropriate matrix 
of multipliers for the general solution, i.e. one omitting the environ- 
mental components that depend on the design of the experiment, is 
given in table 17. 
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We can now calculate the least squares estimates of the three 
parameters in the usual way and from them predict the expectations 
of the nine equations. The analysis of variance of the sums of squares 
of deviations of observed from expected is given in table 18. 

The only significant item in the above analysis of variance is the 
heterogeneity of components, i.e. D, H, and F;, over seasons. There 


TABLE 18 
Analysis of array 7 for height (exclusive) 





Item N MSS 





Residual interaction . ‘ , 6 80°61 





' Between seasons— 











2 Residual interaction . ‘ 6 121°22 

5 Components. : ‘ 3 2315°87 

0 

30. 

o ee 

3 | Within seasons— 

He Residual interaction . . 12 42°19 
Components . : . 6 136+52 

















is no significant residual interaction and the heterogeneity of residual 
interaction over seasons, which borders on significance, is probably 
merely a reflection of the significant heterogeneity of components. 
As the analysis of seasons separately reveals no significant disturbances, 
the standard errors of the components in the two seasons are based 
on the sum of squares of block differences of the observed statistics 
for nine degrees of freedom (table 19). 











TABLE 19 
The components of variation for height in array 7 in 1952 and 1953 
Component 1952 1953 
D .. : - | 44°5907- 3°2035 61-1961-+ 5°1225 
H, . 2 ; | 38-1536-+20°1910 94°4135+32°2854 
Fy . . , 39°3485-+14°7165 102°3600-+23°5314 

















Apart from H, in 1952, which is not significantly different from 
zero at the 5-10 per cent. level of probability, all other components 
of variation are significant. Furthermore, there is no evidence of 
overdominance in this non-interacting array, te D-H, is not 
significantly different from zero. This is in marked contrast to the 
analysis of the complete height data where there is significant inter- 
action and significant overdominance. The only other point of note 
is that all the components of variation are greater in 1953 than in 
1952 but only D and F, are significantly so. 
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The analysis of the non-interacting 3x3 diallel does not differ 
from the above in any important detail. Thus the dominance ratio 
is again not significantly different from complete dominance, i.e. 
H,/D = 0-90, and there is no significant residual interaction. 

Flowering time in 1952 shows no significant residual interaction 
(table 15) when tested against its own heterogeneity over duplicate 
blocks or against the pooled heterogeneities of all items. In 1953, 
on the other hand, we have highly significant residual interaction. 
This difference over seasons is confirmed by the joint analysis of the 
two seasons (table 12), where the heterogeneity of residual interaction 
is also significant. Thus these findings confirm those obtained from 
the scaling tests in both seasons. We have, therefore, four tests which 
agree in detecting no interaction in 1952, and three tests (that is all 
except the Wr/Vr regression), which have detected interaction in 
1953. Their non-detection by the Wr/Vr regression is, however, a 
valuable pointer to the type of interaction involved. There can be 
no doubt that it is of a different type from that found in height, which 
is comparable with the complementary gene interactions of classical 
genetics. It in no way appears to lead to spurious inflations or any 
other changes in the dominance ratio. Thus the unit slopes of the 
Wr/Vr regressions show that the estimated dominance ratio is constant 
over arrays, i.e. Wr—Vr = }(D—H,) for the F, generation and 
= 4(D—}H) for the F, and backcross generations, is constant, and 
hence it must be relatively uninfluenced by the non-allelic interactions 
whose incidence varies amongst the arrays (table 1c). The only 
comparable type of interaction in classical genetics which would fit 
these observations are duplicate genes. Why these interactions 
should be confined to the 1953 season, or alternatively, why they 
should remain undetected by all available methods in 1952 we cannot 
say as yet. But it is possible that the smaller total variation for this 
character in 1952, which is approximately }th of that in 1951 and 
1953 is at least partly responsible for this failure. On the other hand 
there may be a genuine absence of interaction in 1952, which must 
then be ascribed to genotype-environmental interactions which differ 
in the two seasons. 


6. CONCLUSIONS 


The present paper in the series on diallel crosses has two primary 
aims. Firstly, it aims to extend the theory and illustrate its application 
to F, and backcross generations derived from a set of diallel crosses ; 
and secondly, it aims to see how far conclusions and predictions from 
the analysis of parental and F, means are borne out by the fuller 
analyses made possible by the inclusion of these later generations. 

The conclusions of these investigations into the character height 
are easily summarised—in all essential details there is complete 
agreement between the assessment based on parents and F,s alone 
and that based on the analysis when extended to F,s and backcrosses 
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This close agreement is nowhere more striking than in relation to 
the genetical control of heterosis for this character. Every test 
designed to detect non-allelic interaction has given a positive result 
with the height data. Furthermore, the tests which allow one to 
trace the origin of the non-allelic interaction to particular crosses or 
arrays, ¢.g. the Wr/Vr regressions and the scaling tests, have in all 
cases picked out the same combinations of parental lines as those 
responsible for giving rise to the interactions. We may further note 
that the inclusion of F,s and backcrosses in these analyses has in no 
way altered our conclusions concerning the role of overdominance. 
In the joint analyses in both seasons we have found no evidence of 
overdominance after omitting from the analyses all crosses in which 
we have detected non-allelic interactions by means of the scaling 
tests and the F, Wr/Vr regressions. We can take it, therefore, that 
all the heterosis is the result of interactions of a complementary type 
between alleles at different loci. With such consistency over genera- 
tions and seasons the diallel analysis of parental and F, means provides 
an accurate and rapid method of assessing the potentialities of the 
various F, combinations as well as providing valuable information 
about the genetical control of heterosis. The ready availability of 
such information cannot but seriously affect future methods of utilising 
this heterosis for economic purposes. 

Turning to flowering time, we again find good agreement between 
the various methods of analysis, but for this character the agreement 
is confined to within seasons, there being significant differences between 
seasons for the genetical control of this character. Thus while all 
tests agree that non-allelic interaction is absent in 1952, the same 
tests consistently show its presence in 1953. The F, Wr/Vr regression 
alone of all the available methods of detecting non-allelic interaction 
has failed to detect non-allelic interaction for flowering time in 1953. 
Since this was the only test available in 1951 there is the possibility 
that non-allelic interaction played a significant role in the inheritance 
of flowering time in that season also. This failure of the F, Wr/Vr 
method has been traced to its inherent insensitivity to the type and 
magnitude of the non-allelic interaction found in 1953. Thus while 
this method is highly sensitive to the complementary type of inter- 
action, as found in our own height data, it is not so sensitive to duplicate 
gene interactions which we have now found in the flowering time 
data in 1953. For flowering time, therefore, the diallel analysis of 
parental and F, means provides an incomplete assessment of the 
potentialities of the various crosses. Furthermore, the high degree 
of instability of the genetical control of this character over seasons 
makes an assessment based on any one season of doubtful value, 
although, of course, the ability to detect these seasonal differences 
and further to partition the changes in genetical control amongst the 
components of variation is an important contribution of this method 
of analysis. 
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7. SUMMARY 


1. The theory of the diallel analysis of parental and F, means 
has now been extended to the F, and backcross generations derived 
from a diallel set of crosses. 

2. The joint analysis provides sufficient statistics to give least 
squares estimates of the components of variation, i.e. the additive, 
dominance and environmental effects, and their standard errors. 

3. A number of methods of detecting non-allelic interaction are 
given, including the regression of array covariance on array variance, 
which is applicable to any generation, the joint scaling tests and the 
homogeneity of the least squares estimates of the components of 
variation over statistics. 

4. If the within-family variances of the segregating generations 
are also included in the analyses we can detect linkage and estimate 
its effect. 

Illustrations of the analyses are drawn from our own data from 
Nicotiana rustica. This comprises an 8 x8 diallel, the parents and F,s 
of which have been grown in three consecutive seasons and the F,s 
and backcrosses in the last two. 

These analyses have shown that : 

(i) All tests agree that heterosis for the character height is the 
result of a complementary type of non-allelic interaction. 

(ii) The differences in magnitude between the components of 
variation for height in different seasons and different blocks within 
seasons is a genotype-environment interaction dependent on soil 
differences. 

(iii) The analysis of parents and F,s alone in any one season 
provides a satisfactory assessment of the genetical control of the 
character height that is completely borne out by subsequent genera- 
tions and seasons. 

(iv) There are significant differences in the genetical control of 
flowering time ‘n the two seasons 1952 and 1953, primarily relatable 
to differences in the weather. These differences involve not only 
variation in the magnitude of the components of variation but also 
differences involving the presence of duplicate gene interactions in 
1953 and their absence in 1952. 

(v) Linkage between alleles which exhibit reinforcing dominance 
relations is detectable in the flowering time data of 1952 and involves 
at least four factors. 
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In problems of quantitative inheritance it has been customary to 
make genetic analysis of populations using covariances between 
relatives, the method outlined in the classical paper by Fisher (1918). 
Associated with this approach is the notion of partitioning the total 
genotypic variance into additive and non-additive genetic components. 
The importance of this treatment of quantitative inheritance to the 
genetical theory of natural selection has been made clear by Fisher 
(1930), and the practical applications of these concepts to plant and 
animal improvement have been increasingly realised. 

In recent years techniques involving “ diallel crosses” have been 
used in problems which concern quantitative inheritance. The 
methods used and problems attacked by this technique have been 
diverse. Sprague and Tatum (1942), Henderson (1948, 1952), and 
Griffing (1953) have defined and applied the notions of general and 
specific combining ability to plant and animal experimental material 
using a variety of diallel crossing methods, but without an exact 
generalised genetic interpretation of the combining ability effects 
and variances. Hull (1946, 1952), Griffing (1950), Jinks (1954) and 
Hayman (1954a, 19546), again using diallel crossing systems, have 
given procedures for estimating other genetic parameters in terms of 
restricted gene models. 

The present paper is an attempt to give, in terms of population 
genetics, a generalised treatment of the major problems which occur 
with the use of diallel crosses, to integrate wherever possible the 
already existing techniques, and to show the relationship of the diallel 
crossing method to Fisher’s method of covariances between relatives 
as expressed in terms of additive and non-additive genetic variances. 

We shall use the term “ diallel crosses” to describe a procedure 
in which a set of p inbred lines are chosen and crosses among these 
lines are made. There are a maximum of f? possible crosses, which 
can be represented by a px matrix with elements x,; such that x;,; 
represents the i** inbred, x,;(i--j) represents the F, between the i'* 
and j'* inbreds, and x;; represents its reciprocal. Thus, the /? 
combinations can be divided conveniently into three groups: (1) the 


b(p 


p inbreds themselves ; (2) one group of nee F,’s ; (3) a group of 


P(p—1) reciprocal F,’s. 
2 
Experimental methods utilising diallel crosses may vary, depending 
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upon whether or not the inbreds and/or the reciprocal F,’s are 
included. With this basis for classification, there are four possible 
methods : (1) inbreds, one set of F,’s, and reciprocal F,’s are included 
(all p? combinations) ; (2) inbreds and one set of F,’s are included 


but reciprocal F,’s are not (Aes combinations) ; (3) one set of 
F,’s and reciprocals are included but not the inbreds (p(p—1) com- 


binations) ; (4) one set of F,’s, but neither inbreds nor reciprocal 
F,’s are included (Mea) combinations 


Theoretically the important difference among the methods is 
whether or not the inbreds are included. Since the term “ diallel ” 
has in the past been associated by other authors with the methods 
which include inbreds, we shall use the term “ modified diallel ”’ to 
designate methods (3) and (4), which do not include the parents. 
It is with these methods that we shall be most concerned. 

Once the experimental data has been obtained by use of one 
of the above methods, the problem is to estimate certain genetic 
parameters of the population from which the inbreds are derived. 
The problem may be considered in three major parts: (1) an 
examination of assumptions and conditions which are necessary if 
valid inductive inferences are to be made from the experimental 
material about the parent population ; (2) a specification of estimable 
genetic parameters and their genetic interpretation using a generalised 
gene model; (3) a description of the available methods by which 
experimental data are used to test hypotheses and estimate population 
parameters. 


1. ASSUMPTIONS AND CONDITIONS NECESSARY FOR 
VALID INDUCTIVE INFERENCES 

When a set of lines are used in a diallel crossing system for the 
purpose of estimating genetic parameters of the population from 
which the lines were derived, a chain of assumptions which relates 
the experimental material to the original population must be made. 
Starting with the experimental material, we must assume that the 
set of inbreds is a random sample from a population of inbred lines 
which in turn were derived from the parent population by means of 
an inbreeding system free from forces which change gene frequency. 

In addition to making these assumptions we must, before we can 
completely specify the assumptions and conditions necessary for valid 
inductive inferences, (1) determine the restrictions which must be 
placed on the parent population to ensure that it can be truly repre- 
sented by diallel crossing among its derived inbreds, and (2) determine 
which of the four possible diallel methods gives unbiased estimates 
of the population parameters. 

We shall investigate the first problem by considering the particular 
case in which the limiting frequencies of the homozygous genotypes 
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are derived by imposing a self-mating system on an arbitrary popula- 
tion. This population consists of genotypes involving two alleles at 
each of two loci which may or may not be linked, i.e. recombination 
value a<}. We shall start with arbitrary initial genotypic frequencies 
and obtain the limiting frequencies of homozygous genotypes using 
the methods described by Fisher (1949). The self-mating system is 
employed because it involves a minimum number of mating types 
and because the mating type frequencies are equivalent to the genotypic 
frequencies in any generation. We are generalising a solution given 
by Nelder (1952). 


TABLE 1 


Generation matrix with self-fertilisation relating the genotypic frequencies in one 
generation with those in the next generation * 
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Genotype /|Frequency| Yo So | to Up % | we | % | No | Ze 
A,B,/A,B, h Weeeee eae ce ee a 7 - 
A,B,/A,B, ry oO I o o oO fe) 3 } oe 
a |b 
A,B,/A,B, Sy o oO I o } o } oO » a adi 
ao | b 
A,B,/A,B, ty oO oO (9) I oO } o i" 7 rs 
A,B,/A,B, uy oO oO o Oo 4 o o te) a = 
A,B,/A,B, a, te) oO o o oO 4 to) o a < 
A,B,/A,B, W, oO o o oO oO oO } o | a 
A,B,/A,B, xy oO oO o o ts) (3) o } a « 
5? a* 
A,B,/A,B, Ii oO oO oO Oo c) oO oO oO = = 
a? 5? 
A,B,/A,B, 21 o o to) o te) o o o | siz 








* “4” is the recombination fraction such that a+6 = 1. 


The generation matrix (denoted as A) relating the frequencies 
in one generation to those in the next is given in table 1. Solution 
of the determinantal equation |A—AI | =o yields the latent roots 
and their corresponding principal components as given in table 2. 
Associated with each of the multiple roots A = 1, } are four pair wise 
orthogonal vectors which taken together form a 4x10 matrix of 
rank four. This matrix premultiplies its appropriate A matrix to 
give a 4X 10 zero matrix. 

Using these latent roots and principal components we obtain the 
limiting frequencies for the homozygous genotypes as n—>oo. These 
are given in table 3. 

If the initial genotypic frequencies in the parent population are 
those that would occur under equilibrium conditions with random 
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mating, and if the gene frequencies for A,, A, B, and B, are fy, po, 
9, and gq, respectively, then 
(freq. of A,A,B,B,) > pity (freq. of A,A,B,B,) — fogs, 


(freq. of A,A,B,B,) > 195, and (freq. of A,A,B,B,) —> foq3. 




















TABLE 2 
Latent roots and their associated principal components 
Principal components | 
Latent 
roots 
q r Ss t u v w x » zZ 
( A= I I I I I I I I I I 
ia J) B= I —! I —! 1} o| o|—1 fo) t) 
a | C= I 1 1 I Oo} 1 {—1 oO fe) 0 
D = | (1+2a) | (1+2a) | —(1+2a) | —(1+2a)| 0} 0} o| o | (b—a)|—(b—a) 
( E= oO oO oO f) 1 |-—1 Oo}; Oo o 0 
a J) F= oO oO oO oO oO oO 1 |—1! oO ty) 
_ | G= oO 0 fe) oO I I I I 2 2 
H = oO oO o oO —1I|-—!I I I oO 0 
ina I= 0 0 0 0 oo} o| o| o I I 
A= os jJ= 0 o 0 o o| o| o| o I -I 









































It is obvious that the genotypic frequencies occurring in the original 
population can be obtained by diallel crossing among the derived 
inbreds (if their limiting frequencies are taken into consideration). 


TABLE 3 


Limiting frequencies of homozygous genotypes 





Genotype Frequency 
B,/A,B - oe a: =. 
A,B,/A,B, fo tot uot dee + | |e | 2 
A,B,/A,B as otal [<. 
2B,/A,B, bts rete thet | Jot parm Zo 
A,B,/A,B .= =a Bee Ss 
BalaBs seo sot buethitet | | et [ta |e 
ABIAB, | Wo totHutint |] y4[ 1], 





2(1+2a) 


It is equally clear that with populations which cannot be con- 
sidered in equilibrium with random mating, diallel crossing among 
the inbreds will not generate the original population, and thus valid 


inferences are not possible using this technique. 
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With regard to the second problem, the author is indebted to 
Sir Ronald Fisher for pointing out that when random samples are 
drawn from the population of inbreds it is necessary to employ the 
“‘ modified diallel ” crossing systems to obtain unbiased estimates of 
the population parameters. We must defer to the appendix a demon- 
stration of this fact since we consider parameters not yet defined and 
estimating procedures not yet described. 


2. SPECIFICATION AND INTERPRETATION OF 
GENETIC PARAMETERS 


Since the diallel analysis is applicable only to random mating 
populations in equilibrium, the following discussion will be restricted 
to such populations. However, in other genetic respects we wish to 
consider a completely general situation. Thus, we want the genetic 
analysis to deal with an arbitrary number of alleles at each of an 
arbitrary number of loci. We want, in the final analysis, to interpret 
the genetic parameters in terms of a completely generalised gene 
model which will allow any magnitude of additive genetic, dominance, 
and epistatic effects. 

The genetic parameters of interest are the additive and non- 
additive genetic components of the parent population genotypic 
variance. We assume that the population phenotypic variance is 
partitionable into various genotypic and environmental components 
as follows : 


o*, = 0° +07, and since o?, = o?,-+o0?y,q, then o?p = 07, +07x, +07, 
where o*» = population phenotypic variance. 

o*,, = population genotypic variance. 

o*, = additive genetic variance. 


o*y, = non-additive genetic variance. 


o*x, may be partitioned into dominance and epistatic variance, as 
will be shown later. 

The concept of partitioning the total genotypic variance into 
various genetic compounds is due to Fisher (1918) and used by him 
in evaluating covariances between relatives. Fisher (1918) also 
originated the completely generalised genetic model which will be 
used by the present author in connection with these components. 
A model developed by Kempthorne (1955), which is adapted to 
random mating populations and which allows a complete orthogonal 
partitioning of the total epistatic variance, will also be used. 

The estimation of the additive and non-additive genetic com- 
ponents will be made from the experimental material in terms of 
general and specific combining ability variances. Therefore, it is 
necessary to start with random mating populations and (1) review 
the definitions of additive and non-additive genetic variances, (2) 
define the combining ability variances, and (3) show the relationship 








36 B. GRIFFING 


between the two sets of parameters in terms of a completely generalised 
gene model. 

In this section we shall use the summation convention in which 
repeated suffixes (in those terms which involve gene or genotypic 
frequencies) imply summation over these suffixes. For example, 


Pipidis = ZAp;,pjd;; and p,d;; = Xpjd;; etc. 
a t 


(a) Definition of additive and non-additive genetic effects 
and variances 


The definitions will be reviewed first for genotypes at one locus 
only, and then generalised to include genotypes at two or more loci. 
TABLE 4 


Genotypic values and frequencies for a random mating population 
of genotypes derived from alleles at one locus * 


Marginal 
A, Ag As ~ _ -- Aw frequencies 
a neue wane a : ; Se eee Se ee 
Ay | dy dy, dys _ aie = dim pi 
Ag | ds, dso dy5 ”- se = dom pa 
As ds, dsp ds, ad - = dsm bs 
A m din 1 diy 2 dy 3 a - bo dinm Pm 
Siemans os aoe ee : — eer ae | Sew seen 
Marginal 
frequency pr po Ps - - -- Pin I 








* Genotypic frequencies are the product of appropriate marginal frequencies ; i.e. 
frequency of d;; is pip;. 


A notation similar to that used by Kempthorne (1955) will be used. 

1. DEFINITIONS FOR GENOTYPES AT ONE Locus. Consider m alleles 
A,, Ag, ... A,, with frequencies f;, po, ... P, respectively. Genotypic 
values, measured as deviations from the weighted population mean, 
together with their associated frequencies are given in table 4. The 
genotypic value of A,A; is defined to be d;;, and such that p;p;d;; = 0, 
i,j == 1, @, ... Me 

We define the additive genetic effect of A; as a; = pj,d,; These 
effects are subject to the restriction p;a; = 0. We define the non- 
additive (dominance) effect of the A,;A; genotype as 8,; = d;;—a;— a, 
and these effects are subject to the following restrictions: 8;; = 3,, 
and /,5,; = 0 for all 7. The 6;; may be thought of as the interaction 
between the A; and A, alleles. 

In this way we derive the following model which represents the 
genotypic value in terms of additive and non-additive genetic effects 
d;; = a;+a;+8,;. The elements in this model have mean zero and 
are uncorrelated in the population. 
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The total genotypic variance is partitionable into variances due 

to additive and non-additive genetic components as follows : 

bifid? = 2p;0? +Pp8,7- 
These variances may be represented symbolically as o?, = o?, +07) 
where o*,, = p;p,d;? is the total genotypic variance, o?, = 2p,a,” is the 
additive genetic variance, and o7, = p,p,d;7 is the non-additive 
(dominance) genetic variance. 

We note that the a’s can be obtained by the least squares procedure. 
Assuming that the genotypic values can be expressed by the model 
d;; = a;+a;+8;; we choose a’s such that the sum of the weighted 
squared deviations from the linear model (a;+a,;) is a minimum. 
The a’s so defined, together with the associated restrictions, are 
identical to those defined above. 


TABLE 5 


Genotypic values and frequencies for a random mating population 
of genotypes derived from alleles at two loci * 

















A,B, A,B, A,B, - _ o- AwB, 

A,B, 11941 119449 11413 et ag ee im4in Pi 

A,B, 11491 11492 11493 ” — = im4an P92 

AmBy mi4n1 mi2na mi2n3 = Se _ amdan PmQn 
Marginal 
frequency pi P1492 Prds = = = PmQn I 

' 
° Np ol frequencies are the product of appropriate marginal frequencies ; i.e. 
frequency of jd, is pips9n9- 


2. DEFINITIONS FOR GENOTYPES AT TWO OR MORE LocI. We shall 
consider two genetic models under this heading: (a) a model in 
which loci effects are additive, and (6) a model in which the loci 
effects are not additive. 

(a) Additive loci effects. When the loci effects are additive the 
total genotypic variance is the sum of the separate genotypic variances 
for each locus. To illustrate, consider the two loci case in which 
there are alleles A,, A,, ... A,, with frequencies ~,, po, ... p,, at one 
locus, and alleles B,, B,, ... B,, with frequencies q,, go, ... g, at the 
other locus. 

Genotypic values, measured as deviations from the weighted 
population mean, together with their associated frequencies are given 
in table 5. The genotypic value of A;A;B,B, is defined to be ;;d,, 
and such that $,),9,9, ::¢,, = 0. In the following notation, suffixes 
placed before the constants pertain to the A locus and those after 
pertain to the B locus. 

c2 
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We define the effects at the A locus to be : 
1% = Pj i1dny = additive genetic effect of A,, 
9% = Pind iy = additive genetic effect of A,, 
9 = 9491 4441-14 — 5 = interaction (dominance) of A; and Aj. 


The restrictions on these constants are 
p; a = 0, ,8 = 58 and p; ;5 = 0. 


Similarly, we define the effects at the B locus to be : 
On = PiPids isn = additive genetic effect of B,, 
a, = PiPide ise = additive genetic effect of B,, 
82 = Pibs 944. —%,— 4, = interaction (dominance) of B, and Bi. 


The restrictions on these constants are 
95%, = 0, 8, = 8,, and g,8,, = 0. 


The model is thus derived to represent the genotypic value of 
A,A,B,B, in terms of additive and non-additive genetic effects and 
may be written 

iy, = 5a+ja+,8+a,+0,+8,,. 


After the total genotypic variance has been partitioned into 
additive and non-additive components for each locus the loci com- 
ponents may be summed, i.e. 


PP 19491 154? er = 2D; 27 +29,07 +P hs 38? +9.9:5,7, 
or symbolically o2, = o*,-+07y, where 
026 = £:P59491 134,,2 = total genotypic variance, 
o*, = 2p; ,a®+29,a,? = total additive genetic variance, and 
oy, = Op = Pip; 38? +4949, 5,,2 = total non-additive genetic variance. 


When loci effects are additive the extension to any number of 
loci each having an arbitrary number of alleles presents no difficulty. 

(b) Loci: effects not additive. Fisher (1918) introduced epistatic 
constants to account for the interaction of genotypes at one locus with 
those of a second locus. The epistatic constants and their frequencies 
for genotypes formed by alleles at each of two loci are given in table 6. 
The epistatic constants are defined as 


8€xr = 14441 —47— 52 — 8 —2,—2,—8,, 
and by rewriting this equation we arrive at the representation of the 
genotypic value of A;A,;B,B, by the following completely general gene 
model : 


9g, = a tya +48 ta,+a,+8,,+5€,2- 
The total genotypic variance is partitionable as follows : 
PiDiViI 1 ister? = Dj §P AQ pO AD Ds 138? +949 84? +PiP190491 i5€ RD" 
which may be written symbolically as 


_ a 
07g? = on? +ong? = on? +0p?+0/2. 
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oy,a” is broken down into variances due to two different types of 


interaction, namely : 
op” = p;p; ;8?+9,9,5,,2 = interaction between alleles at the same 
locus, and 
OP = PiPididi ii€xe = interaction between genotypes at different 
loci. 
Imposing the restrictions 


i9€kt = 42€ kts G4" KL = i9€tks PiPs i5€u1 =O ANA 94q, j5€4, = O 
on the 16 constants in table 6, results in four independent constants 


TABLE 6 


Epistatic constants and their frequencies for genotypes derived from 
two alleles at each of two loci * 

















Marginal 
B,B, B,B, BB, BB, frequency 
2 
A,A; 11€11 11€12 11€21 11€22 pi 
A,A, 12€11 12€12 12€21 12€22 Piha 
A2Ay 21€11 21€12 21€21 21€22 Pofr 
2 
AiA, 22°11 22€12 22€21 22€22 Pe*s 
Marginal 
frequency q\" 1192 97291 q2" I 








* The frequency of the constant je; is Pifj9e7 


which we arbitrarily choose to be ;;¢,,. The genotypic variance may 
then be rewritten in the following form given by Fisher (1918). 


Pibi9K 9 i4nr? = 2P; 547 +29,0,2 +piP; 8? +949 842 
I 


4P1P 29192 
+21) 290°(P1 11€22 +2 22€292)?+2h0°9192(91 22€11 +92 rates)*| 


where 


Z = (£191? 11€11 P1792” 11€22—f 2791" 22€11 +h 2742" 22€22)- 


|2t+-2p.%e00(01 11€11 +92 11€22)” +2P1h2919(P1 11€11 +P2 22€11)” 


(b) Definition of general and specific combining ability effects and variances 
and their relationship to additive and non-additive genetic effects and 
variances 

In an analysis of diallel crossing experiments it is convenient to 
redefine the genetic model in terms of general combining ability 

(g.c.a.) and specific combining ability (s.c.a.) effects and variances. 

To give a general definition of g.c.a. and s.c.a. effects we need 
only consider the representation of the zygotic array of a random 
mating population in equilibrium as the square of the gametic array. 
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This representation may take the form of a two-way multiplication 
table (for example table 5) in which the multiplying elements forming 
the row and column margins are the gametes together with their 
respective frequencies, and the products in the table represent the 
genotypes and their frequencies. We define, then, the g.c.a. effect 
associated with the i‘" gamete, say, as the weighted mean value of 
genotypes in the ,‘" row of the multiplication table. Thus, if the 
value of the genotype in the ,;‘" cell of the multiplication table is d;; 
(measured from the population mean) with frequency /;f;, the g.c.a. 
effect for the ,‘" gamete is g; = fjd,; The s.c.a. effect for the ,;"" 
genotype is 5;; = d;;—g;—g;. 

By use of these effects we can represent the ,;" genotypic value 
as a linear function of the average value associated with each of the 
gametes which combined to form the genotype plus the interaction 
effect between these gametes, i.e. d;; = g;-++g;+5;;. 

We have defined the combining ability effects in terms of gametes. 
When homozygous lines are used in diallel crossing techniques the 
combining ability effects are often defined in terms of the inbred 
lines. In this case the terms are equivalent since there is a one-to-one 
correspondence between inbreds and gametes. 

To show the exact relationships between combining ability and 
additive and non-additive genetic effects we shall first consider 
populations of genotypes involving only two loci and later generalise 
to populations of genotypes which involve any number of loci. 

1. POPULATIONS OF GENOTYPES INVOLVING Two LocI. We shall 
again consider the two situations in one of which the loci effects are 
additive and in the other not additive. 

(a) Additive loci effects. Let us consider a population of genotypes 
derivable from an arbitrary number of alleles at each of two loci as 
given in table 5. 

We define the g.c.a. effect associated with the A;B, gamete (or 
with the inbred A;A;B,B,) as ,2, = (1,9, :d,,- By substituting 


4, = a +ja+;8+a,+0,+8,, 


in the equation for ,g, and employing the restrictions associated with 
the elements in the model, we are able to determine the exact relation- 
ships between g.c.a. effects and additive genetic effects, which are as 
follows : 

iSk = iA TO, 


with restrictions = 0. 
421 = ja+a, | Pu vow 


The s.c.a. effects are defined as ,j5,, = ;j€,,—j;2,—42, and when 
appropriate substitutions are made we find that the s.c.a. effects 
are a function of the dominance effects, i.e. ;45,, = j8+8,,. The 
restrictions are jj5,) = jiS1xs Ska = siSkas ier = 5i8ix9 ANA Psd, is5k1 = O 
for each ik. 

The genotypic value of A,;A;B,B, in terms of g.c.a. and s.c.a. 














| 
| 











ST, ee 


DIALLEL CROSSES IN QUANTITATIVE INHERITANCE 4! 


effects is iter = 2p tii tisSke The total genotypic variance can be 
partitioned into g.c.a. and s.c.a. variances as follows : 


P1991 04a? = WGK Bie TOP isk 


Symbolically this can be represented as 0%, = 207... +07... and 
hence the following equivalence exists between combining ability 
variances on the one hand and additive and non-additive variances 
on the other : o?, = 207... and o?,) = o?...,, 

(b) Non-additive loci effects. Consider the more restricted population 
of genotypes derived from two alleles at each of two loci. The 
combining ability effects are defined as in (1a). However, the 
additional complication of dual epistacy necessitates the following 
genetic model: ,,d,, = ,a+,;a+,8+a,+a,+8,,+;;€,,, and, therefore, 
the g.c.a. effects now include a linear function of the epistatic effects 
as follows : 





I 
2, = ,a+a 
aes a . si ta” 


I 
&2 = 1a +a,— ——2Z 82 = 2a+ag+ 
, ? , 4P192 ™ ' , ree qa 
where Zz = (491? 11€11—P 19a 11€22—f0"1" 22€11 +h 2°92" 22€22)- 


Twice the g.c.a. variance contains not only all of the additive 
genetic variance but also a portion of the ieaing variance, 7.¢. 


181 = 19-+a,+ 





20 = 2p; ,a?+29,a 2° 
= Pe ah tegen oan” 
(z is a function of the linear xlinear interaction contrast in a 3? 
factorial experiment—see Cockerham (1954)). 
The s.c.a. variance (residual genotypic variance) includes all of 
the dominance and the remaining epistatic variance. 
Thus, when dual epistacy is taken into consideration, the partition- 
ing of the population genotypic variance is as follows : 


a, = o* 4 +o7) +o?, or a, ame e." +e." 


However, 2c, ..,2 = o,+4(linear x linear portion of o?,), and o,..,7 = 
? A 1 


g.c.a. 
o*,+(residual o,), where o?, is the total epistatic variance. 


2. POPULATIONS OF GENOTYPES INVOLVING MORE THAN TWO LOCI- 
To generalise the analysis to any number of alleles per locus and 
any number of loci, it is convenient to consider the gene model 
developed by Kempthorne (1955) which allows a complete orthogonal 
partitioning of the total epistatic variance. For definitions and an 
elaboration of the use of the model we refer the reader to Kempthorne 
(1955). We shall briefly consider the results when this model is 
applied to the concept of general and specific combining ability. 
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In populations of genotypes derived from an arbitrary number of 
alleles at each of two loci the following holds : 

(1) The genotypic value of A;A,B,B, as measured from the 
population mean is expressed as : 


yy = 54 +5a+58 +a, +0,+8,,+;(a0),+;(aa),+ (a0), 
+4( aa) » +; (08) gp +508) gy His (Ba) p+ is(Sa)  +55(88) x 


(2) The population genotypic variance is partitioned as : 
og = 04 +0"p +074, +0,p+0"pp» 


where o?, == additive genetic variance, o?, = dominance variance, 
o*,, = additive x additive epistatic variance, 7,» = additive x domin- 
ance epistatic variance, and o?pp) = dominance x dominance epistatic 
variance. 

(3) The g.c.a. effect associated with the A,B, gamete is 


iS, = 12 +a,+,(aa),, 


where ,a and a, are additive genetic effects, and ,(aa), is an additive x 
additive epistatic effect. 

(4) The g.c.a. barmsary IS Op ca? = Pj (a? +9,0,7 +P)9, {(aa),? which 
symbolically is o,....2 = $07, +}07,,. 

(5) The population genotypic variance can be partitioned in 
terms of the combining ability variances as follows : 





1." _ o*, +4o7 44 
Gai” Deed oD + dotn, +o*, D + opp 
Pq =a? +0? p+ O,+07,n +0" pp 


There is no difficulty in generalising to any number of loci each 
of which has an arbitrary number of alleles. The corresponding 
partitioning of the population genotypic variance can be written 
down as follows : 





a =o", +07 44 +4o7 naa +etc. 
— 2 2 : 2 2 

Csca. = a7) +4307, 4 +07 4p +O pp t+Fo7 aga toraapt+ Sxspp +o ppp + ete. 

0". = a, +07p+ Pra toanteppt+ Ora t+Osap tO app +o"ppp tete. 


We note that o,.,2 is equal to the covariance between parent and 
offspring in a random mating population at equilibrium. This 
illustrates the genetic connection between the diallel method and the 
method involving covariances between relatives. 


3. METHODS OF TESTING AND ESTIMATING GENETIC 
PARAMETERS FROM EXPERIMENTAL DATA 


We assume that the experimental material represents a random 
sample from some specified random mating population about which 
inferences are to be made. 

Analyses for the two different “ modified diallel” experimental 
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methods will be given. For convenience we shall designate the 
analysis as follows : 


Analysis (1) : For the experimental method in which ee") F,’s and 


their reciprocals are included. 
Analysis (2) : For the experimental method in which only one set of 
Men) F,’s are included. 

In both analyses we assume that a replicated experiment has been 
conducted which yields the necessary phenotypic mean values and 
a pooled error variance (with 2 degrees of freedom) associated with 
these mean values. We shall not consider the more complicated 
situation resulting from disproportionate sub-class numbers which 
has been excellently treated by Henderson (1948). 
re 


Analysis (|) (Involving cel F,’s and their reciprocals) 


We assume that the peat value of the ij‘* observation can 
be represented by the model : 
Kig = BE; +85 +55 +s +errOF, 
where 7, 7 = 1, 2, ... p, (i #7), and, 
» = population mean. 
g;(g;) = g.c.a. effect of the z**( j*") inbred line. 
5:3 = 8.C.a. effect associated with the cross of the i** and j" inbreds 
and such that s;; = 5;;. 
r;; = reciprocal effect associated with the crosses of the i** and 
j* inbreds and such that 7;; = —1;;. 
error = environmental error effect associated with the ij" observation. 
We assume that these error effects are not correlated with 
each other nor with the various genotypic effects and that 
they are normally distributed with mean zero and variance 
o,. o is a pooled error mean square obtained from the 
replicated experiment and has n degrees of freedom. 

The analysis of variance and expectations of mean squares are 
given in table 7. In determining the expectations of mean squares 
we assume that g,, 5;; and 7,; are random variables with zero means 
and variances o,’, o,? and o,?._ The construction of this analysis of 
variance is given by Yates (1947) and Kempthorne (1952). 

To test the hypothesis o,? = 0, 


use F = po waa [o(o—2) -— 2) n| degrees of freedom. 
e 
To test o,2 = 0, 
use F = vd with o(6—3) n| degrees of freedom. 
Me 2 


L 


To test o,2 = 0, 


use F = a with (o—1), op—3)| degrees of freedom. 
§ L 
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The estimation of the components of variance is accomplished by 
equating the observed to the expected mean squares and solving. 



































Thus, 
~>o Mr—Me ~., Ms—Me ra Mg—Ms 
oc? =. — ,o,° = ———- and o? = —2>_.. 
: ; 2(p—2) 
TABLE 7 
Analysis of variance for the experimental method in which F,’s 
and their reciprocals are included 
| 
Source DF Sum of squares * | brand ‘ae 
ae: es - ——|— ———— 
2+ 20,?+ 
R = as ee a | M Cc ° 
(g) p-i apoaye +X i) p(p—2) X?. § 2(p—2)a,? 
TL (%45 +%p4)*@— = (Xo +X.4)*) 
R(s) b(p—3) [; rh sti) 2(p—2) . ) | Ms o,?-+20,? 
2 SJ 2 
+————~ X 
| (p—1)(p—2) °° | 
R(r) P(p—1) 224(x;;—%;;)? | Mr o,"+20," 
2 ij 
icy 
Error from 
replicated n Me a,” 
experiment 
* where X.. = 22%;;. 
i: j 
TABLE 8 
Analysis of variance for the experimental method in which one set of F,’s 
only is included 
| . | Mean Expectation of 
Source | DI Sum of squares * | squares| mean squares 
R _ ee, hee ee M o,*+0,7+ 
(s) ii j=" b(p—2) xe | ‘ (p—2)o,* 
R(s) | P38) | aax,2-—— sex+ ——_*>_ x2] Ms | o,8+0,9 
ra eee (p= 1)(p—2) ** atte, 
i<j 
Error from 
replicated n | Me 0,” 
experiment | 




















* where X.. = 22x;;. 
ij 


i<j 


Analysis (2) (Involving one set of 2 F,’s only) 


We assume that the ij‘ observation may be represented by 


xij = M+8; +83 +5;5 + error, 








where i<j = 1, 2, ... p. The elements in the model may be defined 
in the same way as in analysis (1). 
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The analysis of variance which partitions the total sum of squares 
into orthogonal reductions is presented in table 8. This analysis is 
essentially the same as that first presented by Sprague and Tatum 
(1942) and later corrected by Federer (1951). The expectations of 
mean squares are found in the same way as in analysis (1). 

To test the hypothesis o,2 = 0, 


wie F = Ms with [Pe=a) ; n| degrees of freedom, 
Me 2 


and to test o,? = 0, 


use F = Be with |(o— 1), ee 3)) degrees of freedom. 
s 


The estimation of components of variance is accomplished as 
follows : 

‘ 4 Mg—Ms 

6,2 = Ms—Me, and 6? = —s P 


ps 





4. SUMMARY 


When a diallel crossing system is used in quantitative inheritance 
for the purpose of estimating genetic parameters of the population 
from which the inbreds were derived, the following assumptions and 
conditions must be met : (1) the parent population must be a random 
mating population in equilibrium ; (2) the experimental set of lines 
must be a random sample from a population of inbred lines which 
were derived from the parent population by the imposition of an 
inbreeding system free from forces which change gene frequencies ; 
(3) a “ moditied diallel ” crossing system must be used in which the 
lines themselves are not included in the experimental set. 

The additive and non-additive components of the parent genotypic 
variance are estimated by use of general and specific combining ability 
components of variance. The exact relationship existing between 
these two sets of parameters is given using a completely generalised 
gene model. 

When interpreted in terms of the classical method of covariances 
between relatives, the method of diallel crosses yields estimates 
equivalent to those obtained by covariance between parents and 
offspring. 

Analyses of variance are given which may be used for testing and 
estimating the population general and specific combining ability 
variances from experimental material. 


5. APPENDIX 


We wish to demonstrate by considering a particular case that use of the “‘ modified 
diallel ” crossing system gives rise to unbiased estimates of the population para- 
meters. We shall start with a random mating parent population in equilibrium 
consisting of genotypes involving two alleles at each of two loci, which may or 
may not be linked. 
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It has been shown that if an inbreeding system is imposed on the entire population 
(in the absence of forces which change gene frequency) the population of derived 
inbreds will have the following distribution : 











Genotype Designation Frequency 
A,A,B,B, Ly | Ai = tin 
A,A\B,B, Ly H 2 = Pde 
A2A,B,B, L; | fy = Pot 
A:A,B,B, Ly | 4 = Pods 








We assume: (1) the population of inbreds is indefinitely large and, therefore, 
the frequencies of the inbred lines are not changed by sampling without replace- 
ment ; (2) random samples of size are drawn from the inbred population, allowing 
repetition of lines ; (3) for each sample of lines a ‘‘ modified diallel ” crossing system 
is imposed, and for convenience the method which involves one set of crosses but 
not reciprocals will be used. 

Specifically it is necessary to show that (1) the frequency of genotypes generated 
by this experimental procedure, when averaged over all possible samples and when 
calculated after sample frequencies have been taken into consideration, will yield 
the genotypic frequencies of the original parent population, and (2) the expectations 
of G,? and G,* taken over all samples are the true population parameters o,* and 
o,*. (Symbolically this condition is indicated as E,(¢,*) = o,, and E,(@,*) = o,°). 


TABLE 9 


Frequencies of genotypes resulting from a modified diallel crossing system using the 
lines indicated in the column entitled “ Typical sample.” 
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The particular illustration with which we are concerned necessitates making 
a complete enumeration of samples of size four (this is the smallest possible number) 
and deriving the expected values of the statistics in question over all such samples. 

The sample frequency distribution may be obtained from the multinomial 
(ALi +feLetfebs3+f,L,).4 We shall denote a sample having L,, L;, L, and L; as 
Sixr It is convenient in tabulating frequencies to group the samples according 
to the partitions of the number four. Thus (4) corresponds to all samples S,,;; in 
which all four lines are L; (¢ = 1, 2, 3, or 4), (31) corresponds to all saniples § 
in which three lines are L; and one is L, (i, 7 = 1, 2, 3, or 4, but i #7), etc. 

Table 9 gives the frequencies necessary to show that the “ modified diallel ” 
crossing system does yield genotypic frequencies which, when averaged over all 
samples, are identical to the genotypic frequencies occurring in the parent popula- 
tion. To verify this, one must multiply the genotypic frequencies within each 


ttij 
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sample by the corresponding sample frequency and add over all samples. The 
frequency terms can be reduced from the fourth to the second degree if use is made 
of the identity (f,)* = 1. 

To demonstrate that the “‘ modified diallel’’ crossing system yields unbiased 
estimates of the population combining ability variances, it is necessary to perform 
an analysis of variance of the type given in table 8 on F,’s derived from each of 
the 35 different possible samples. For these analyses it is assumed that a genotypic 
value, as measured from the population mean, may be represented by the model 
dis = 8;+8;+5;3;, Where d,; is the genotypic value of the cross between the i** and 
j™ inbreds, g,(g;) is the g.c.a. effect of the i”(j*) inbred, and s,; is the s.c.a. effect 
of the ij cross. We then determine the composition of ¢,? and @,? for each sample 
in terms of squares and cross-products of the elements in the model, and finally 
show that when all samples are taken into consideration the following is obtained : 
E,(é,*) = 2fige? = o,*, and E,(¢,*) = Life® = a," 

To illustrate the procedure of evaluating ¢,? and ¢,* for a given sample, let us 
consider the class of samples designated S,,;; (i, j = 1, 2, 3, or 4, buti<j). We 
need perform only one generalised analysis of variance for this class of samples 
and then obtain particular solutions for each sample by assigning appropriate 
values for i and j. 

The mean values for F,’s derived from the sample of inbreds (S,;;;) by the 
imposition of the ‘‘ modified diallel ” crossing system may be recorded as follows : 


‘ 


. * fg 














i 
oT) Mad | a © ad 
t Mig XG xj. 
J * 43 xj 
a %j 
where : Xe = Lat&ot+San, 25, =torj 


Hee = 48s HIPs +See +25 yy 
Hye = 2B +485 — +25 qy+545 
and Xe. = 69, +68; +54 +4545 +545 


The problem is to determine the expressions for ¢,* and G,* in terms of the gg 
and s,, elements using the analysis of variance obtained from table 8 by putting 
el 

The mean square for R(g) is found to be 


Dy 
Mg =3 (12%. = =) = $[2(8s—2;) + (See—5yy) }*. 
The mean square for R(s) is 
x2 
Ms =} ( Eqs —4Ea% + +) = Hou—25y+5y)* 
a<b 3 


Using these mean squares we determine ¢,? and é,? as follows : 


&,* = 4(Mg—Ms) = 3[(¢s—83)* + (Sea—Ses) (Sts —544) + (8s —83) (Se — 53s) ] 
and &,? = Ms = $(54,—2543+54;)*. 


Thus, for each of the six possible samples which can be represented by S,;;; 
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TABLE 10 


Composition of c,* in terms of squares and cross-products of specific combining ability elements 
Sor each class of samples corresponding to the partitions of the number four 





Number 





Partition — . ct * Pn ll Square and cross-product terms 

‘ ua) f Siti ¢ a 7 ”" F None i 
“eo | ee | | ee a 
@) Siiss P 6 | fit ents Sig +543)? : ; 


(21) Sisie 12 I2f2fife | ¥(Sis—5is— Sie t5in)? 


(14) Sizer I 246: LS fe | 8U(Sis + Sin +5ir + Spe +551 +5n1)?—3 Se3(Siat 
SirtSju +551) —3 Sin(Sea+Sjx+5e1) — 
3 SéalSsa+Se1) —3 Spu(SpetSer) —3 SsaSea] 























TABLE 11 


Composition of G,* in terms of squares and cross-products of general and specific combining 
ability elements for each class of samples corresponding to the partitions of the number four 





Number 
«| Sample f Sample Ss d 
SS- tt 
Partition ty a a frequency quare and cross-product terms 
— | ——— Sn aUnnea: er" ee: 
(4) Siti 4 Si None a 


(31) Siti 12 4 fifi H(gs—gs) + (5-515) ]? 








(2?) Siiis 6 6f2f;* | $L(gs—Bs)? + (sis—Ses) (Sis—545) + (8s 283) 
(Sia 543). 
(21*) Siise 12 12fi7fiSe | Pola 8i(8i—8s—Su) +3(8s—8u)? +4 818+ 


2 S¢i(Sig+Siu—2 Spx) +54j(Sig—6 See +2 55x) 
+Sin (Sie +2 55x) +2[2 8: (Sis —S55x) 
—8i(S4e—2 Sig +2 Siu—Sjx) —Su(Ses +2 545 
—2 Sin—Sjx)]} 





oe Sisk 1 24 fifi | 1xt3(g2+g;?+en2+91°) —2(gigst+8isetgigr 
+858 +8581 +8n8 1) +545 (Set Ser +592 +551) 
+ Seu (Sea + 55+ See) +540 Sea(Spa-b Ser) 

+ Sy 552+ 5x1) +55 Set — 4 (Seger FT SeHS41 

+S 418sx) +2184 (Sis+ 54 +541 —-532—- 551 — Se) 
+85(S45—Ste—S4a+5je+551—Set) — 

— 8x (Sis —Sie+541—Se+551—Se) — 
8i(Sig+Sie—Sit+5jx—5Sj1—Ser) }} 
































DIALLEL CROSSES IN QUANTITATIVE INHERITANCE 49 


(4.¢. Syre2> Sigs» Siraa> Sees» Seoag, amd Sy544) we find that the values for ¢,? and 
cg,” are obtained by substituting the appropriate values for i and j into the generalised 
expressions given above. 

In a similar manner we can obtain generalised 6,2 and G,? expressions for each 
class of samples corresponding to the five partitions of the number four. These are 
given in tables 10 and 11. 

To show that &,? and ¢,* are unbiased estimates of the population parameters 
o,* and o,* (i.e. E.(é,?) = g,* and E,(@,*) = &,*) we may (1) expand tables 10 
and 11 to give a complete enumeration of all possible samples, (2) multiply the 
frequencies of square and cross-product terms within a given sample by the sample 
frequency, and (3) sum over all samples to obtain the weighted mean values for all 
square and cross-product terms. 

We are required to demonstrate that these expressions for E,(¢,*) and E,(é,) 
which contain both square and cross-product terms can be transformed into 
expressions containing squared elements only, and, in fact, to prove that 


E,(é,*) = 2fige and E,(,*) = ZF, si 
ij 


Let us consider E,(@,*) first. Since coefficients of terms involving g.c.a. effects 
only and coefficients of cross-product terms involving g.c.a. and s.c.a. effects are all 
zero, we need consider only square and cross-product terms involving s.c.a. effects. 
If to these terms we add the following five expressions each of which equals zero, 


2 f(Zfjsy)*, for i = 1, 2, 3, and 4) 
J L 
and (—) (2HA; 5u)* | 
J 
we find that E,(6,*) = 2/ifi 5°; = o,*. 


i, 
To determine E,(4,*) we must consider square and cross-product terms of 
g.c.a. effects only, of s.c.a. effects only, and cross-product terms involving both 
g.c.a. and s.c.a. effects. 
By use of the identity, Zaft "BX(ZS)* = 82 fSi8i8sn it is possible to show that 


the square and cross-product terms involving the g.c.a. effects only, when summed 
over all samples, yield simply Xf,g;? which is o,?. 

The over-all contribution to E,(@,?) of square and cross-product terms involving 
only s.c.a. effects can be shown to be zero, if to the coefficients we add the following 
expressions each of which is equal to zero : 


(ACh 54)" for i = 1, 2, 3, and 4 | 
| and hd 545)*. | 


Likewise, the over-all contribution to E,(j0*) of the cross-product terms 
involving both g.c.a. and s.c.a. effects can be shown to be zero, if to the existing 
coefficients we add the following expressions each of which is equal to zero : 


(- 2) Fie ZS su fori = 1, 2, 3, and 4) 
| a(S fad) Ohh Sii) | 


When all terms are taken into consideration we find that E,(¢,*) = fige? = 9,. 

This completes the demonstration that for inbred samples of size four, the 
“modified diallel *? crossing system yields unbiased estimates of the population 
parameters. 
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In most plant species in which reproduction is effected both by cross- 
and self-fertilisation, the relative contributions from these two kinds 
of matings depend to some extent on climatic and other environ- 
mental conditions and so vary from one generation to the next. 
In the present article, however, we shall confine our attention to the 
simplified but fundamental situation in which these contributions 
bear a constant ratio to one another throughout all generations. 

Consider a large population of plants reproducing in non-over- 
lapping generations in such a way that there is a constant probability s 
that any plant will be self-fertilised and a probability 1—s that it 
will cross with some plant chosen at random from the population. 
We shall suppose that all crosses are equally fertile and all genotypes 
equally viable. With this system of mating, gene frequency clearly 
remains constant. Suppose A and a are the genes present at some 
locus and f(A), p(a), [p(A)+p(a2) =1] their frequencies in the 
population. If »,(Aa) denotes the frequency of occurrence of the 
genotype Aa in the nth generation, we have 

busr(Aa) = 45p,(Aa) +2(1—s)p(A)p(2) 
and for the genotype AA, 
busr(AA) = sfp,(AA) +4p,(Aa)] +(1—s)p(A)*. 
Hence, if we put 
K, = pa(Aa)— #2—9 pcayp(a), 
then 
K, = (45)"K,. 

At equilibrium, K = 0, 


ie. py (Aa) = lim p,(Aa) = 489) p(ayp(a) 
and p, (AA) = lim p,(AA) = p(A)[1— 22—9 p(a)| 


s 
= p(A)*+ —* p(A)p(a). 


These equilibrium genotypic frequencies are those which would exist 
in a composite population comprising a random mating component 


2(1—s) 


supplying the proportion of the whole population and a 


5r 
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selfing component making up the remainder. Genotypic equilibrium is 
not attained in one generation as with random mating but is approached 
at a rate measured by —log,(4s). 

If B and 6 are two genes present at a second locus and » is the 
frequency of recombination between the two loci, the frequency of 
occurrence of a given chromosome or gene-combination such as AB 
in the genotypes of the (n-+-1)-th generation is 

buss(AB) = p,(AB) +4yLp,(Ab/aB) —p,(AB/ab)]. 
It can also be shown that 
bav1(Ab/aB) —Py.1(AB/ab) = }(s—2y)[p,(Ab/aB) —p,(AB/ab) ] 
+2(1—s)[p(A)p(B) —p,(AB)]. 


uy, = p(A)p(B) —p,(AB) 
Yn =P n(Ab/aB) —p,(AB/ad), 


()-2-() 


where Z stands for the matrix 


Hence, if we put 
and 


we have 





a me 
2(1—s) $(s—2y) 
The matrix Z has two positive and unequal latent roots A, ~ which 
satisfy the equation, 
—(2+s—2y)z+s(1—2y) = 0. 
The latent roots are therefore 


= }{2+s—279+-/[(2—s—2y)?+8y]}. 


Clearly, if neither 1—s nor_y is zero, A(> ) is less than unity and 


fim bP, (AB) = p(A)p(B) 
iim [p,(Ab/aB) —p,(AB/ab)] =o. 


A measure of the “a of approach to equilibrium is given by —log,A 
(see table 1). 
Corresponding with the two latent roots A and p, there exist two 
principal linear components of frequency, 
L, = 2(1—s)u,—(1—A)v, = A*.L, 
M,, = 2(1—s)u,—(I—p)2, = »".M, 
from which the values of u,, and v, can readily be determined for any 
positive integer n. 
When s = 1, 2.e. with complete selfing, the latent roots are 1 
and 4—y, the corresponding principal components of frequency being 
(1-+2y)u,—yv, = constant, 
and v, = (4—7)".v,. 
Hence, in this case, 


$.(AB) = p,.(AB/AB) = p,(AB) + —— aoe —fo(AB/ab)] 
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giving the limiting genotypic frequencies when complete homo- 
zygosity has been attained. In general, when s = 1, 


P.(AB) + p(A)p(B). 

With complete random mating, i.e. s = 0, the latent roots are 
1—y and o, the corresponding principal components of frequency 
being 

u,—dy ~~ (1—y)". [u,—3y . v,] 
u,—}0, =O (n>0). 
Hence 
b(A)p(B)—p, (AB) = 3[P,(Ab/aB) —p,(AB/ad)] 
= (1—y)""*[o(A) p(B) —p,(AB) —4p,(Ab/aB) +3p,(AB/ab.)] 


TABLE 1 


Rate of approach to equilibrium at two linked loci with recombination frequency y 
in a population with 100 s per cent. self-fertilisation and 100 (1—s) per cent. 
random mating 





Recombination frequency y 


















































\0°95 | ‘0044 0081 | ‘0114. | 0142 | -0166 | -0188 | -0207 | -0224 | 0239 | -0253 | 0266 


0°05 | O-10 | 0°15 | 0°20 | 0°25 | O-g0 | 0°35 | 0-40 | 0°45 | 0°50 | 0°55 | 0-60 

(0-00 | :0513 | "1054 | +1625 | +2231 | +2877 | +3567 | 4308 | +5108 | +5978 | -6932 | -7985 | -9163 
0°05 | *0499 | ‘1022 | +1572 | -2152 | +2764] +3412 | +4099 | +4831 | +5611 | -6444 | +7337 | ‘8291 
o-10 | 0484 -0989 | "1517 | +2069 | +2647 | -3253 | -3888 | -4553 | -5250 | +5978 | ‘6737 | *7524 
0°20 | 0452 | -0918 | +1398 | +1894 | +2403 | +2925 | +3459 | +4002 | +4553 | -5108 | +5664 | -6213 
0-30 =| -0416 | -0839 | +1269 | +1705 | +2144 | +2584 | -3023 | -3459 | -3888 | -4308 | +4715 | +5108 
{ 0°40 | 0376 | -0753 | «1129 | +1502 | +1870 | -22g1 | +2584 | -2925 | +3253 | +3567 | +3865 | -4147 
0°50 =| *0331 | 0657 | 0975 | +1285 | +1584 +1870 | +2144 | +2403 | -2647 | 2877 | «3092 | +3292 
0-60 | 0281 | 0551 | -0809 | +1054 | +1285 | +1502 | +1705 | +1894 | -2069 | -2231 | -2382 | -2523 
0°70 | :0224 | 0433 | 0628 | -0809 | 0975 | +1129 | +1270 | +1398 | +1517 | +1625 | +1725 | +1808 
080 | -0159 | -0303 | :0433 | 0551 | 0657 | 0753 | :0839 | -0g18 | -098q | -1054 | 1113 | +1165 
0°90 =| -0084 | ‘0159 | 0224 | 0281 | -0332 | -0377 | 0416 | -0453 | 0484 | 0512 | -0538 | -0537 
*0277 





With mixed selfing and random mating (s+1), the genotypic 
frequencies at equilibrium are given by 


p,,(AB/ab) = p,,(Ab/aB) = 
4(1—s)p(A)p(a)p(B)p(5)/[2—s +25 (1 —9)] 
P.(AB/Ab) = —p,,(AB/ab) +-4(1—s)p(A)p(a)p(B)p(6)/(2—s) 
p,,(AB/AB) = 
3p,,(AB/ab) +-p(A)p(B) —2(1—s)p(A)p(B)[ p(2) +0(4)]/(2—s). 


The equilibrium frequencies for the other genotypes are readily 

obtainable from these expressions on making appropriate gene 

substitutions. We see that when equilibrium is attained, linked factors 

are associated to an extent that depends not only on the proportion 

of selfing but also on the magnitude of the recombination fraction. 

In such a case, it might be possible to determine the recombination 
D2 
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fraction between two loci from an analysis of the population frequencies, 
. When » = 3, as for two loci on different chromosomes, 


_-:16s(1—s) 
po(AaBb)—p(Aa)p_(Bb) = =!) p(A)o(ap(B)000) 
p,,(AABb)—p,,(AA)p,(Bb) = Gay ayn MANOA (BOL) 

—» (AA)p (BB) = —4(!—5) _ 9(A)p(a)p(B)p(6). 
p(AABB)—p,(AA)p,.(BB) = — 4°" p(A)p(a)p(B)0(0) 


Two factors which segregate independently are thus not associated 
at random when equilibrium is reached. There is, in fact, a positive 
association between the homozygous (or heterozygous) states at loci 
on different chromosomes. This association may be represented as 


TABLE 2 


Association between the homozygous and heterozygous conditions at two independent 
loci in a population experiencing mixed selfing and random mating 








AA Aa aa 
BB I 2 1 
Bd —2 4 —2 
bb I —2 I 




















in table 2 which gives the excess of the equilibrium frequencies for 
genotypes at two independently segregating loci over the product of 
the equilibrium frequencies at the separate loci as a multiple of 


_45(1—5) _ ’ 
(4—s)(2—s)? P(A) p(a)p(B)p(). 


For given gene frequencies, this excess has its greatest value 
0: 1506 p(A)p(a)p(B)p(6) 
when s = 0°6946. (See fig. 1.) 

There are results analogous to these for randomly mating popula- 
tions of polysomic organisms where, for our present purposes, double 
reduction may be considered as equivalent to a small amount of 
selfing. It has been shown (Bennett, 1954) that at equilibrium in 
randomly mating populations of tetrasomic or hexasomic organisms 
the probability of occurrence of any chromosome or gene-combination, 
such as AB, equals the product of the appropriate gene frequencies, 
i.e. p,,(AB)=p(A) p(B). If there is no double reduction, linked loci 
will therefore be associated at random in the gametes or zygotes at 
equilibrium. However, when double reduction can occur, linked loci 
may not be associated at random in the gametes and zygotes at 
equilibrium though loci on different chromosomes certainly will be. 
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SUMMARY 


In a population that experiences mixed selfing and random mating, 
genes are associated at random on the chromosomes at equilibrium 
but there is a positive association between the genotypic state (homo- 
zygous or heterozygous) at different loci even when the two loci are 
on different chromosome, i.e. there is a tendency for heterozygosis to 
become concentrated in some members of the population. 


0:15 - 


0-10 


45 (I-s) 
(4 ~ s)(2-s)* 





0:05~ 
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Fic. 1.—Excess of the equilibrium frequency of double homozygotes for two independent 
loci over the product of the equilibrium frequencies of single homozygotes at these two 
individual loci expressed as a multiple of the gene frequency product p(A)p(a)p(B)p(d) 
and as a function of the amount of self-fertilisation. 
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|. INTRODUCTION 

Two closely-related species of Acacias, A. decurrens or green wattle 
and A. mollissima or black wattle, are grown in the Union of South 
Africa for tannin production. Green wattle, although generally 
considered hardier under adverse conditions, does not produce so 
high a quality of tannin. The small proportion of green wattle 
previously under cultivation has therefore now been largely replaced 
by black wattle. 

From the genetical point of view, however, green wattle has one 
marked advantage : the seed ripens some four months after flowering, 
as opposed to fourteen months in black wattle. Philp and Sherry 
(1946) therefore chose green wattle for their preliminary work on 
the degree of natural crossing in wattle and they were able to show 
that the species was largely cross-pollinated. Because of the close 
resemblance between the two species, they suggested that black wattle 
would also prove to be largely cross-pollinated. Selfed progeny from 
several black wattle trees have now been grown and scored for recessive 
seedling characters and it has been possible to calculate the approximate 
amount of natural crossing in many of these trees. 


2. METHODS 


Selfed and natural (open-pollinated) progenies were grown from 26 trees from 
various parts of the wattle-growing area. The method of bagging followed that of 
Philp and Sherry (loc. cit.), cellophane bags being used throughout. Bagging was 
confined to trees near the borders of plantations since the small crowns, inaccessibility, 
and the high proportion of bags damaged by wind made bagging of more central 
trees a hazardous and impracticable proposition. 

Seeds for sowing were made permeable by removing a small chip of the seed- 
coat from the distal end. This method of preparing the seed was found to give a 
higher and more even germination than the more usual treatment with hot water 
and this factor greatly facilitated the scoring of seedling characters. As far as 
possible, selfed and open-pollinated progeny were grown from seed collected in 
the same year, but in some instances additional natural seed collected in the following 
year was sown. 


3. FERTILITY AND GERMINATION 


The inflorescence is a panicle of some hundreds of globose head 
each consisting of 35-40 closely packed flowers. An average-size 
inflorescence consists of well over 20,000 flowers. Even with open 
pollination, only a small fraction of these flowers form pods, but it is 
obvious by inspection that selfed inflorescences produce, on the 
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average, fewer pods than open-pollinated inflorescences of the same 
size. 

Since counts of flowers and pods set are not feasible, the number 
of viable seeds per pod was used as a measure of fertility. Sections 
of young flower buds showed that the great majority of developing 
pods had two rows of seven ovules each. In a small percentage only 
thirteen ovules were counted but in no instance where this occurred 

TABLE 1 
Fertility and germination following self- and open-pollination 















































Self-pollinated | Open-pollinated 
] 
Tree number No. of Mean | Seed | Germin- N. of | Mean Seed | Germin- 
ods | Per | set | ation | ods | Per set ation 
pes pod | per cent.| per cent.| poes | pod | per cent.| per cent. 
ees at SEE ~ = | _ I- = 
Sweetwaters 1 203 | 3°03 21:7, | -7Q't | 40 5:30 | 382 | 89°4 
99 2] 258{ 1:08 77 | 637 fo) 412 | 295 | 62-7 
99 3| 17| 2°76 19°77 | 90°7 40} 4°70 | 336 | go-7 
‘ 4 10 | 1:00 71 | 85°7 40 | 4°80 | g47 | 853 
” 5 | 15 | 2-02 157 | 95" 40 | 560 | 403 | Qg6-0 
” 6) 42] 1°88 13°4 | 95°3 49 | 6-30 | 452 | 8971 
Yarrow 46 167 2'10 147 | 898 50 | 7:50 | ° 53:3 | 913 
‘a 48 37 | 1°80 125 | 64-7 50| 6:20 | 4ggo | 83-4 
; 52| 41| 260 | 185 | 43-7 | 50| 5:20 | 376 | 49:2 
99 53 9 | 280 19°8 | 84-0 50| 460 | 328 | 95:7 
Mowbray 66 50 | 1°72 m9 | «81-7 20; 815 | 582 88-0 
a 67 Not recorded | 47:1 Not recorded 540 
Bloemendal 6 52 | 1°84 13772 | «O25 50 | 6-08 43°4 78-0 
; 7 4| 5°75 | 42 | 956 50 | 8-06 | 57-6 93°2 
; 9 38 | 218 15°60 86-7 50 | 8or 581 93°5 
‘ 14 5 | 5°40 38-6 82:0 50 | 916 65°4 80-9 
9 15 | 260} 4°63 B3°1 93°5 50 | 7°30 §2'1 g2°6 
im 17 19 4°26 30°4 | 764 50 | 8-90 63°6 88-o 
5 18 II 3:00 214 | 95:0 31 7°84 54°6 96-6 
” 20 39 | 3°20 22" 97°8 50} 9:24 | 66:0 95°7 
9 25 53 | 3°83 27°4 g2°6 36 | 5°00 35°7 100°0 
99 27} 53} 3°57 25°5 87°4 50 | 8-42 60-1 93°8 
» 2 10 | 4°30 30°7 88-4 50 | 6°84 48-8 go'9 
9 34 | 133] 2°69 19°2 89-7 50 | 6-62 47°3 95°4 
| | | | | 
| 1526 | 2-72 19°5 | 853 | 1027 | 6°75 48:2 | 86-7 
| | | | | | 





could it be determined with certainty that an ovule had not been 
lost in sectioning. Counts were also made of seeds and ovules in 
mature pods from go trees. Except in the case of one tree, the highest 
number of seeds, or seeds and undeveloped ovules together, in some 
3500 pods examined was fourteen. The exception was tree No. 44, 
Bloemendal, in which out of 50 pods examined two with sixteen seeds 
each were found. 

It would seem, therefore, that the great majority of pods have a 
potential seed set of fourteen and that the occurrence of sixteen seeds in 
tree No. 44 was an unusual and individual characteristic of this tree. 
Consequently fourteen has been used in the calculation of fertility, 
which is based on the seed set as a percentage of this figure (table 1). 








Ss eC eS 


weV 





NATURAL CROSSING IN BLACK WATTLE 59 


The three groups of trees, i.e. Sweetwaters, Yarrow and Mowbray, 
and Bloemendal differed significantly in their behaviour as regards 
both fertility and germination (fig.). 


45 


25 


65 


45 


25 


PERCENTAGE SEED SET 


65 


45 


25 


= 


A 


LN 





B 








1 j 1 1 j 





50 


60 70 80 90 100 
PERCENTAGE GERMINATION 


Fic.—Graph showing the relationship between seed set and germination following open- 
pollination (open circles) and self-pollination (black circles) in individual trees of 
black wattle. 
A. Sweetwaters group. B. Yarrow and Mowbray group. C. Bloemendal group. 


In all groups, self-pollination resulted in a marked reduction in 
the mean number of seeds per pod. The possibility that this reduction 
was due, in part at least, either to an environmental effect induced by 
bagging or to inadequate pollination under the bag could not be 
entirely ruled out. 


Other observations, however, suggested that the 
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reduction in fertility was not due to these causes, but was a real effect 
of self-pollination. When cross-pollinations were made under bags, 
heavy yields of pods and seeds resulted, while rubbing the inflorescence 
lightly with flowers from the same tree did not increase the yield of 
seed. In all three groups, there were significant differences between 
individual trees, but there was no correlation in the percentage seed 
set from self- and open-pollination on the same tree. 

In the Sweetwaters trees there was no significant difference in 
germination between self- and open-pollinated seeds from the same 
tree, but differences between trees occurred. Trees from Bloemendal 
showed no difference in germinating capacity either between self- 
and open-pollinated seed from the same tree or between trees. In 
the Yarrow and Mowbray group, self-pollination brought about a 
significant (P = 0-01) decrease in percentage germination and trees 
also differed in over-all germinating ability. In all groups there was 
a positive correlation between the germination of selfed and natural 
seeds, indicating that germinating ability is a factor largely controlled 
by the parent tree. 

The reasons for the differences between the groups was not clear. 
It may possibly have been due to strain differences, but in’ view of the 
small number of generations of wattle so far grown in South Africa 
and the free interchange of wattle seed between areas, this does not 
seem to be likely. 


4. SEGREGATION AND NATURAL CROSSING 


Of the twenty-six trees from which selfed progenies were grown, 
twenty-one segregated for one or more recessive characters. Families of 
open-pollinated seedlings were grown also and in every instance when 
the selfed family segregated for a recessive, the corresponding natural 
family segregated for the same character, but in a very different ratio. 
The degree of natural crossing in individual trees was calculated 
from the divergence in open-pollinated families from the ratio in 
which a character was shown to be segregating in the selfed family. 
Segregation of recessive seedling characters in both selfed and natural 
progenies and the calculated degree of natural crossing are shown in 
table 2. 


5. DESCRIPTION OF RECESSIVE CHARACTERS 


The cotyledons of black wattle seedlings are normally dark green, 
while the first seedling leaf is a deep bluish-green when fully expanded. 
Variation occurs among the blue-greens, some individuals being 
darker than others, but recessives are always distinct. The first 
seedling leaf following the cotyledons is unipinnate. Stems and 
undersides of the leaves are flushed red or purple with anthocyanin. 


Albino. Albino seedlings had pale cotyledons and the first seedling leaf was 
white with a pink tinge. A few seedlings developed a trace of green but all died 
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within three or four weeks of germination. There was a slight shortage of recessives 
which was homogeneous over all bagged branches x*j5) het. = 1-103. 

Light-green. Light-green seedlings were paler than normal and lacked the 
typical “‘ blueness ”’. 

Yellow-green. Yellow-green seedlings were paler than light-green and the seedling 
leaf had a definite yellowish cast. Segregation was clear-cut and scoring within 


TABLE 2 


Segregation of recessive characters and percentage natural crossing 








wv 








l 
: . Minimum 
Segregation : Segregation 
Tree number Character paid in selfed Ratio d P in saa natural 
families naan families hestpagcn, 
per cent. 
Yarrow 46 | Blue-green, albino 206 5I gis "10—*05 484 | 18 86 
Yarrow 48 | Blue-green, light 23 9 Bil ‘go — ‘80 124 2 94 
green 
Yarrow 53 | Blue-green, yellow- 10 II git <-‘o1 275 II 85 
green 
Bloemendal 20 | Blue-green, yellow- 96 34 gil -go — ‘80 891 | 51 78 
green 
Bloemendal 17 | Blue-green, yellow- 62 6 I5i1 *70—*50 794 | 12 76 
green 
Bloemendal 7 | Blue-green, yellow 18 4 git *70—*50 1015 | 36 86 
Bloemendal 32 | Blue-green, yellow 80 25 3:1 *7O0— +50 413 | 32 71 
*1212 | 114 66 
Bloemendal 32 | Blue-green 62 18 a2% *70—*50 396 | 42 62 
variegated *1o1g | 193 36 
Yarrow 53 | Blue-green, blotchy 12 9 git *20—"10 273] 13 82 
Bloemendal 14 | Blue-green, blotchy 27 5 3:1 *50—*30 804 5 97 
Bloemendal 4g | Blue-green, blotchy | 47 14 oe -go— 80 go8 | 20 gI 
Ahrens 66 | Blue-green, blotchy 61 4 15:1 -go0— ‘80 468 5 83 
Bloemendal 15 | Normal, withered 231 57 3:1 "05 —°02 974} 29 88 
tips 
Bloemendal 27 | Noimal, shed 77 33 $¢2 *30— +20 96 I 96 
pinnules 
Sweetwaters 1 | Normal, pointed 128 50 $31 *50—*30 12g 4 88 
| pinnules 
Sweetwaters 2 | Normal, pointed 104 | 26+7?] 3:1 *70—*50 107 6 79 
pinnules 
Sweetwaters 6 | Normal, pointed 223 51 S21 "02 —-O1 129 9 74 
pinnules 
Yarrow 52 | Normal, pointed 31 9 3:1 "70 —*50 78 3 85 
pinnules 
Bloemendal 4g | Normal, bipinnate 58 6 I5i1 *5O— 30 898 5 gI 
Bloemendal 34 | Normal, bipinnate | 256 II 15:1 | :20—-10 1023 9 86 
Bloemendal 27 | Red stem, green | 56 17 gs1 ‘80—-+70 92 4 83 
stem 
































* Second collection of natural seed. 


each family easy, but comparison between families was difficult, partly because all 
the families segregating for yellow-green were not grown at the same time and 
partly because subtle differences, difficult to describe, occurred between families. 
Because of these differences it was thought to be unlikely that the same genes were 
concerned in all three families. Segregation in 15 : 1 as well as 3 : 1 ratios occurred. 

In the family from Yarrow 53 the ratio of blue-green to yellow-green was outside 
the generally accepted limits of probability for a 3 : 1 segregation, but gave a good 
fit with a 9: 7 expectation. The numbers in the family were small, however, and 
for the moment it is best regarded as a monofactorial difference. 

Yellow. The yellow recessives in the families from trees Bloemendal No. 7 and 
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No. 32 were very similar in appearance. These seedlings had yellowish cotyledons, 
but the seedling leaf, when it first unfolded, was scarcely distinguishable from normal. 
As the leaf grew older the colour changed rapidly to a clear yellow. Later leaves 
behaved in a similar manner and seedlings showing the character were much 
retarded in growth. 

Pale yellow. In the open-pollinated families from both tree No. 7 and tree 
No. 32 a type segregated which had not occurred in the selfed progeny. Three 
seedlings out of 1015 from tree No. 7 and four out of 1428 from tree No. 32 had 
normal green cotyledons, but a very pale yellow first seedling leaf. As the leaf 
matured it gradually turned green until at the end of a few days it was scarcely 
distinguishable from normal except for a somewhat slaty appearance. 

Although this character appeared in open-pollinated progeny only, some 
estimate of the ratio in which it was segregating may be made. In the selfed families 
from these two trees yellow behaved as a simple recessive. From the number of 
yellow recessives in the open-pollinated families it is possible to place the total 
number of selfed seedlings in each family as approximately 144 and 584 respectively. 
Assuming a 63:1 segregation of pale-yellow to other selfed seedlings the ratio is 
within the 0-05 per cent. limit of probability in both families and when the families 
are added together the ratio fits fairly well with expectation (P = 0:20—0:10). 
The non-appearance of pale-yellows in selfed families can be explained purely on 
the grounds of insufficient numbers. 

Variegated. This character was detected in the first seedling leaf, but showed 
up more clearly in the second to fourth leaves. The expression varied considerably 
and when patches of white were large, some deformation of the leaf occurred due, 
presumably, to unequal growth between green and white areas. 

Blotchy. Phenotypically, blotchy appeared to be the same character in all 
four families in which it segregated. Small brown spots appeared on the first 
seedling leaf, but it was not until four or five leaves had expanded that the character 
reached its full expression. As the leaves unfolded the ends of the pinne curled 
under, the leaf became paler, and the leaf surface developed numerous irregular 
brown blotches. Blotchy seedlings were markedly less vigorous than normal 
seedlings. The character segregated in both 3 : 1 and 15: 1 ratios. 

Withered tips. Tips of the pinnules, which were slightly more pointed than usual, 
died back to a depth of about one millimetre. Seedlings were stunted and weakly 
and many died in the nursery. Assuming the character to be monofactorial, there 
was a shortage of recessives which was homogeneous over all eleven bagged branches 
which made up the family (x%;19) het. 11-407). 

Shed pinnules. In this type the first seedling leaf was small and almost colourless 
with the pinnules curled and bunched together. Subsequently, the pinnules fell 
off without having expanded, leaving the bare midrib of the first leaf. A number 
of seedlings produced second and third leaves which were greatly reduced in size 
and of a pale pink colour. These leaves did not shed their pinnules. None of the 
seedlings survived more than a few weeks. 

Pointed pinnules. The pointed pinnule character did not become obvious until 
seedlings had formed some four or five leaves. Pinnules which were more widely 
spaced than usual tapered rapidly to a point in contrast to the normal rounded 
ends. In the family from Sweetwaters No. 6 there was a significant shortage of 
the recessive type which was consistent over all bagged branches (x%;33] het. 1-658). 

Bipinnate. The bipinnate character appeared as subdivision, to a greater or 
lesser degree, of the pinnules of the normally unipinnate first leaf. The expression 
of the character was very variable. At one extreme all the pinnules of the first 
leaf were sub-divided to give a completely bipinnate leaf, while at the other extreme 
only a few small subdivisions on a single pinnule occurred. Intermediates between 
these extremes were most frequent, one or two pinnules being fully sub-divided 
while the rest of the pinnules remained entice. Subsequently, it has been found 
that the factors responsible for bipinnate first seedling leaf also affected mature 
leaves. a few pinnules being fully or partially sub-divided. Mature leaves showing 
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this character are, therefore, partially tripinnate. At least two genes are responsible 
for this character. 


Green stem. Green stem recessives appeared to be completely free from anthy- 
cyanin on the stem in contrast to the normal reddish colouration. 


6. TWO-FACTOR RATIOS AND LINKAGE 


Families from four trees segregated for two recessive seedling 
characters in each (table 3). 

There was some indication of linkage between yellow-green and 
blotchy in the selfed progeny from Yarrow No. 53, but because of 























TABLE 3 
Two-factor segregation of seedling characters 
Blue-green | Blue-green Yellow-green Yellow-green 
Tree number normal blotchy normal blotchy P 
Yarrow 53 | 8 2 4 7 <o-o! 
Exp. 9:3:3:1 11°7 3°9 3°9 1'3 
Blue-green | Blue-green |Yellow-green 
normal variegated 
Bloemendal 32 61 18 25 ae 0°95—0-90 
Exp. 9:3:4 58°5 19°5 26°0 
Normal Normal Shed 
normal green stem pinnules 
Bloemendal 27 60 17 33 ade 0*50—0-30 
Exp. 9:3:4 61°2 20°4 27°2 
Normal Blotchy Normal Blotchy 
unipinnate | unipinnate bipinnate bipinnate 
Bloemendal 9 42 14 5 I 0*50—0-30 
Exp. 45:15:3:1 43°6 14°5 2°9 | I 
| 




















the low numbers in the family and the abnormal single factor segrega- 
tion of yellow-green, the confirmation or otherwise of linkage must 
await the collection of further data. 

In the other three trees, the two-factor ratios gave good agreement 
with expectation on independent assortment. Yellow was epistatic 
to variegated and shed pinnules to green-stem, resulting in 9:3: 4 
segregations. In tree No. g, blotchy segregated as 3 : 1 and bipinnate 
as 15:1 ratio giving two-factor expectation of 45:15:3:1, with 
which there was good agreement. 

In four families there was strong evidence of linkage between 
pointed pinnules and a second character which appeared at a later 
stage in development, namely dwarf habit. The parents of three 
of these families, i.e. Sweetwaters Nos. 1, 2, and 6, were open-pollinated 
seedlings from the same parent tree. It is not known whether these 
three are related to the fourth tree, Yarrow No. 52. Two-factor 
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segregations of pointed pinnules and dwarf habit are given in 
table 4. 

Single-factor ratios were upset by a higher mortality of pointed 
pinnule than normal seedlings ; analysis of the data showed that 
when homozygous, the pointed pinnule gene depresses viability by 
about 50 per cent., but there was no evidence of the dwarf gene 
depressing viability. There was no evidence of differential germination 
among the genetic classes. 

The four families were homogeneous for recombination value and 
so may be grouped together. Taking the four families, the recombina- 
tion value was 8-7-+-0-94 per cent. Since the recombination classes 
are poorly represented in some families, this value must be accepted 
with caution. It can, however, be said with some confidence that the 
genes pointed pinnules and dwarf habit are linked with rather less 
than 10 per cent. recombination. 

















TABLE 4 

Tree number Normal | Dwarf | Pointed Dwarf 

| ; pointed 
Sweetwaters 1 96 | 6 | 4 21 
» 2 78 | 5 5 10 
a. 6 171 9 I 19 
Yarrow 52 23 | 2 I 6 
Total . ; . | 368 | 22 II 56 











7. NATURAL CROSSING 


The percentage natural crossing which occurs with each of the 
trees was calculated on the assumption that all foreign pollen which 
effected fertilisation carried dominant genes in respect of the seedling 
characters for which the trees were heterozygous. The figure is, 
therefore, the minimum percentage natural crossing and takes no 
account of foreign pollen carrying the same recessive gene, a factor 
which may be of some importance in wattle. The calculated natural 
crossing is given in the last column of table 2. 

In most trees, when more than one recessive character segregated, 
the natural crossing calculated from the segregation of each character 
in the open-pollinated families showed fairly good agreement. From 
Bloemendal No. 32, however, a second collection of seed gave 
anomalous results for one character. This tree segregated for two 
recessives, namely yellow and variegated, both of which gave good 
3 : I ratios in the selfed family. 

From the first collection of open-pollinated seed, the numbers of 
normal to yellow and normal to variegated were 413 : 32 and 396 : 42 
respectively. As expected, the ratio of normal to recessives shows 
good agreement for both characters (x*,,; = 0-903). When the 
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proportion of normals to yellows in the first and second collections 
of seed is compared there is also good agreement, the figures being 
413 : 32 and 1212 : 114 respectively (x*,,; = 0°873). When, however, 
the number of yellows and variegated in the second collection are 
compared, there is wide divergence from the expected ratios, the 
figures being 1212 normal to 114 yellow and 1019 normal to 193 
variegated (x?;,; = 31°972). Similarly, there is a significant difference 
in the proportion of variegated to normal occurring in the first 
collection of seed and that of variegated to normal in the second 
collection. It is evident that variegated is segregating in different 
ratio from yellow in the second collection of seed and from both 
yellow and variegated in the first collection. 

The most likely explanation of this discrepancy is that a neigh- 
bouring tree to No. 32, also carrying the gene for variegation, had 
come into flower or the flowering time had coincided with No. 32 
in the second year, and that a number of the variegated recessives 
were, in fact, due to foreign pollen carrying the same recessive genes. 


8. DISCUSSION 

The estimated mean percentage natural crossing in black wattle 
is of very much the same order as that in green wattle, being 84 and 
85 per cent. respectively. Osborn (1931) gives counts of recessive 
seedling characters appearing in ten trees and, grouping the ten trees 
together, the total of normal to recessives was 23750 : 869. Assuming 
that the characters are all behaving as simple recessives, this gives 
a natural crossing of 86 per cent., again showing close agreement. 

It may be taken, therefore, that the mean minimum percentage 
crossing is in the neighbourhood of 85 per cent. In both species, 
there is wide variation between individual trees, the degree of cross- 
pollination ranging from 35 to 100 per cent. in green wattle (Philp 
and Sherry) and from 62 to 97 per cent. in black wattle. This range 
of variation may be connected with differences in self-fertility such as 
have been shown to occur, but on the somewhat scanty data available, 
there is no indication of any correlation between self-fertility and 
cross-pollination. 

Another possibility is that the apparent low degree of cross- 
pollination in some trees is due to neighbouring trees carrying the 
same recessive gene so that self- and cross-pollination are not completely 
distinguishable. Since wattle seed is heavy and is not dispersed over 
any great distance from the parent tree, the method of re-establishment 
by natural regeneration adopted in the majority of wattle plantations 
lends itself to the formation of groups of trees having recessive genes 
in common. Bateman (1947) has derived evidence from his studies 
of both insect- and wind-pollinated crops which suggests that if seed 
were not widely distributed from the parent tree, there would be a 
correlation between proximity and closeness of relationship. 

It is possible, therefore, that the assumption made when calculating 
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the degree of natural cross-pollination, i.e. that all foreign pollen 
effecting fertilisation carries a dominant in respect of the gene on 
which the calculations are based, may lead to considerable errors. 
A case in point is the anomaly encountered in tree Bloemendal 
No. 32. If the degree of natural cross-pollination in this tree were 
based on the segregation of the variegated recessive in the second 
collection of seed only, the calculated figure would be 36 per cent., 
whereas calculated from the segregation of the yellow recessive in 
both collections of seed or variegated in the first, the degree of natural 
crossing would be in the neighbourhood of 67 per cent. It is possible, 
therefore, that the percentage cross-pollination in wattle is considerably 
higher, in some instances at least, than that calculated from ratios of 
normal to recessives in natural seed. 

Variation. Where there is a high degree of cross-pollination, the 
general breeding behaviour of a population may be predicted with 
considerable accuracy. It is characteristic of such populations that 
they exhibit variability and that the intermediate phenotypes are 
favoured at the expense of the extremes (Darlington and Mather, 
1949). Because of the flow of variability maintained by outcrossing, 
extreme types are not lost to the population, but constantly reappear 
by segregation. Any alteration, whether natural or artificial, in the 
selective force, will favour more extreme types and cause the average 
to shift to meet the new conditions. Similarly, interference with the 
outbreeding mechanism will break down the relational balance 
established by natural selection and cause a new, and probably less 
efficient, balance to be set up. 

As is to be expected from the high percentage natural crossing, 
both green and black wattles conform closely to this pattern. All 
characters examined, including tannin content of the bark, showed 
wide variation, with the majority of individuals close to the mean. 
Some of this variation is undoubtedly due to environment, but from 
the results of progeny trials, it is concluded that much of the variation 
must be explained on a genetic basis. Families consisting of the 
open-pollinated progeny from individual trees differ significantly 
from each other in vigour and also show striking differences in 
morphological characters such as stem form, branching habit, foliage, 
etc. From an economic point of view, some of these families are 
better than average, others are poorer. Under natural conditions, 
the progeny of both extremes will be pulled back towards the average, 
but with artificial selection for the economically more advantageous 
characters, it should be possible to shift the average closer to the 
desired type. 

Again, forced inbreeding has the expected effect. Self-pollination 
in wattle uncovers numerous recessive genes, the great majority of 
which are deleterious to a greater or lesser extent when homozygous. 
At the same time, the relational balance built up by natural selection 
is broken down and a new and less efficient balance forced upon the 
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plant. The result of this new balance is the typical inbreeding 
depression in vigour shown by the progeny of both green and black 
wattles when the trees are self-pollinated (Anon., 1953). 

It is evident that the breeding behaviour of the black wattle 
closely resembles that of green wattle and that both must be regarded 
as typically outcrossing species. Philp and Sherry (1946) have 
discussed in some detail the effect of the pollination mechanism on 
the formulation of a wattle-breeding programme and it is unnecessary 
to reiterate their conclusions here. 


9. SUMMARY 


1. In twenty-six trees of black wattle studied, the percentage seed 
set on self-pollination was consistently lower than that from open- 
pollination. 

2. The progeny of twenty-one of the trees segregated for recessive 
characters in both selfed and open-pollinated families. From the 
divergence from expectation of the ratios of dominant to recessives 
in open-pollinated families, the approximate degree of natural crossing 
was calculated. 

3. Cross-pollination ranged from 62-97 per cent. with a mean of 
84 per cent. and in this respect black wattle was very similar to green 
wattle. 

4. Linkage between two characters, pointed pinnules and dwarf 
habit, was found in four families the recombination value being 
approximately 8-7 per cent. 
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Tue genetic history of any population may be elucidated either by 
historical studies using written, oral and archeological evidence or 
by analysing the present day population to discover, if possible, the 
components of its composition. The two approaches are comple- 
mentary and the primary aim of the present paper is to consider the 
frequencies of the ABO blood group genes in peoples in different 
parts of Ireland and to correlate these frequencies with what is known 
of the history of the population. As that which is new is the ABO 
blood group survey, based on data collected in the course of the 
work of the National Blood Transfusion Association, this will be 
described at length and some account will also be given of the distribu- 
tion of Rhesus-positive and Rhesus-negative people. 


1. AN OUTLINE OF THE HISTORY OF THE 
POPULATION OF IRELAND 


Ireland has no internal geographical] barriers to speak of. The 
mountains are low, the rivers shallow. In general the most productive 
agricultural land is in the east, the most barren in the west, but there 
are few parts which are not inhabitable on a basis of either a tillage, 
a grazing or a fishing economy. 

From archeological evidence it appears that Ireland was first 
inhabited in the immediate post-glacial period. Since these early 
mesolithic times the population has received successive additions by 
seaborne invasion. An inherent difficulty in interpreting archeo- 
logical evidence in terms of the racial composition of the population 
is the fact that this evidence records customs and ways of life and 
changes in these may be due to their natural evolution, perhaps 
stimulated by contact with peoples overseas, to large-scale invasions, 
or to invasions by a small number of people who, because of their 
supremacy, in, say, the art of war, come to impose their traditions 
on the original inhabitants. The evidence of skeletal remains, while 
theoretically of greater value, is so scanty that he would be rash 
indeed who attempted to deduce from it the genetic history of pre- 
historic man in Ireland. 

* Present address—School of Pathology, Trinity College, Dublin. 
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Turning to historic times it would appear that from a.D. 100 to 
A.D. 800 the country was ruled by a Gaelic aristocracy and it may be 
that this rule was over peoples of many distinct races. Prior to the 
Norse invasions of about 800 a.p. there was a unity of civilisation, 
law and language under a Gaelic High Kingship. Ireland was not 
conquered by the Romans. 

The Norse invaders established colonies and fortresses round the 
coast and these were especially strong at Dublin, Wexford, Waterford, 
Cork and Limerick. For two hundred years they tried to subdue the 
country but their attempts ended when Brian Boru routed their 
Dublin garrison in 1014. With the centre of their power destroyed 
the Norsemen’s coastal settlements became isolated communities 
enjoying citizen rights under the High Kingship. 

After a period of unsettled Gaelic rule, Ireland experienced a tide 
of Anglo-Norman invasion and settlement from 1166 to about 1300. 
This flow was mainly into south and east Louth, Meath, the eastern 
half of Westmeath, the eastern fringes of Offaly and Leix, Kildare, 
Dublin, Wicklow, Carlow, Wexford, Waterford, Kilkenny and Cork ; 
and to a lesser degree into Limerick, south Tipperary and north Kerry. 
(Fig. 1 shows the county boundaries.) While much of the west was 
conquered it was never a settlement area and was but sparsely 
garrisoned. In Ulster the area east of Lough Neagh was settled by 
1250 but the area to the west, together with Roscommon and Leitrim 
never experienced much settlement. Donegal was never conquered 
and its resistance was, in part, due to help received from the gallo- 
glasses, mercenaries of mixed Gaelic and Norse descent from west 
Scotland. Many of these mercenaries settled in north-west Ireland, 
and during the fourteenth and fifteenth centuries groups of them 
offered their services to Irish lords and chiefs in many parts of the 
country. Neither the Anglo-Norman nor the later settlements 
penetrated the more mountainous regions of Ireland, such as are 
found in parts of Wicklow, Antrim and south-west Dublin. 

In the fourteenth century there was an eastward movement of the 
Irish, at the expense of the Anglo-Normans, especially into Westmeath, 
Offaly, Leix, Limerick, north Kerry and the country regions of Cork. 
In Ulster the Anglo-Normans were pushed eastwards to the Ards. 

By 1500 the King’s writ ran only in Dublin, south Louth, Kildare 
and Meath, although some other districts were still effectively ruled 
by Anglo-Norman lords. ‘This further contraction of the “ pale” 
was as much due to Anglo-Normans being assimilated into the Irish 
community as to their replacement by the Irish. They were becoming 
Hibernis ipsis Hiberniores. 

In the sixteenth century, under Mary, Leix and Offaly were 
planted with English but the numbers of these settlers were so small 
that their influence on the genetic composition of the population was 
probably insignificant. Under Elizabeth there were unsuccessful 
attempts to establish an English plantation in Munster and such slight 
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English influences as there are to-day in the country regions of Cork 
probably date from this time. 

The seventeenth century witnessed important changes in the 
composition of the population. Under “ the articles of Plantation ”’ 
of James I in 1609 Ulster received numerous Scottish and English 
Protestants. While this plantation did not include the counties of 
Monaghan, Antrim and Down, the last two were soon assimilated by 
Scots who had settled in the Ards earlier in the century. These 

















Fic. 1.—County boundaries and percentage frequency of group B gene. In this and 
subsequent figures the areas X, Y, Z are those of Hart’s survey and show the gene 
frequencies calculated from his data. The present data for Northern Ireland are 
represented by a small circle. 
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settlements were so successful that by 1641 Protestants owned six- 
sevenths of the land of Ulster although they constituted numerically 
a smaller fraction of the total population. In the more southern 
parts of Ireland relatively unimportant settlements were established 
under James I in north Wexford, Longford, Leitrim, south Offaly 
and south Carlow and the adjoining part of Wicklow. 

The Cromwellian settlements, partly reversed under Charles II, 
transplanted many land owners (“Irish Papists ”»—which included 
Anglo-Normans and English) and their retainers from many counties 
to Clare and Connaught with a corresponding influx of English into 
both town and country. 

Apart from these additions the composition of the population in 
certain regions of Ireland may have been affected by racially selective 
losses due to emigration and other causes. Many thousands were 
lost by massacre, deportation and emigration during the seventeenth 
century. During the potato famine of 1846-48 the population fell 
from 8 millions to about 6 millions and these losses were mainly from 
the poorer country regions. 

Recent years have seen the growth of central government and 
industries in Dublin and the development of the county towns. This 
has resulted in migrations to these towns from the country. 


2. GENERAL CONSIDERATION OF BLOOD TRANSFUSION DATA 


The blood group data acquired by a transfusion service can be 
biassed in various ways. Appeals for blood may be made to particular 
population groups (say factory workers, or soldiers) ; there may be a 
conscious or unconscious selection of group O volunteers (“ Universal 
Donors”) ; blood donation clinics may be set up only in the areas 
of densest population or greatest propinquity to the headquarters of 
the service—hence rural dwellers may be under-represented in the 
data of the area sampled. Resignations from the donor panels may 
not be at random, or lists of donors and their groups when “ brought 
up-to-date ” may exclude some original volunteers on a blood group 
basis. Depending on the laboratory techniques used, there may be a 
tendency to group certain ABs as Bs, and there may be some bias 
towards a “ safety-first’ attitude in designating a doubtful blood to 
whatever category has least clinical danger. This is, however, more 
likely in Rh grouping where a doubtful Rh-negative is likely to become 
Rh-positive in the blood bank so that it may not be given to an 
Rh-negative recipient, while in ante-natal clinics the opposite tendency 
will be at work and a doubtful blood will be classed as Rh-negative 
because this is the safest practice from the mother’s point-of-view. 


3. METHODS USED IN COLLECTING THE DATA 


When this survey began the National Blood Transfusion Association already 
had details of the blood groups of 10,000 newly-volunteered Irish donors. However, 
it was not certain that a small percentage of these did not already know their blood 
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groups due to a wartime blood grouping campaign in certain towns, and because 
there might have been differential recruitment of group Os all this data was 
abandoned. The following scheme was then followed. 

At donor sessions all new donors (i.e. those who had never given blood before 
and were now giving a full donation for the first time) were asked whether they 
knew their blood group. If they did not the following information was recorded : 
1. Mr, Miss or Mrs. 2. Surname. 3. Christian name. 4. Present address. 5. 
Occupation. 6. Birthplace in Ireland. After their blood had been to the central 
laboratory of the Transfusion Service, there was added 7. ABO group and 8. Rh 
group (D group only). In the records department this information was typed on a 
card of one of four colours according to the ABO group. This “survey card” 
was different from the one used for the master card index in the central record of 
the donor panel, and so could be removed from the records office without causing 
confusion. The master card was also coloured according to the ABO group and 
as the survey card was typed by the same typist at the same time there was little 
chance of error. 

The survey cards were then sorted by county of birthplace. ‘‘ Birthplace ” was 
chosen in preference to “‘ present residence *’ for the following reasons : 1. It move: 
the geographical location of the blood groups backwards in time by twenty to 
forty years. This will tend to reduce the blurring effect of the considerable country- 
to-town population movements of the last thirty years. 2. It gives a country-wide 
cover in a short time because many Irish blood donors live or work in different 
cities, towns or counties from those in which they were born. As there are many 
counties in which the National Blood Transfusion Service does not hold donor 
clinics because of either distance, thinness of population, or the fact that local blood 
transfusion schemes are already in operation, a survey based on “‘ present residence ”’ 
would mean special expeditions to particular counties or else would involve reliance 
on second-hand data from sources where the laboratory methods or the bias in 
collecting could not be fully checked. 


4. LABORATORY METHODS 


High titre selected grouping sera, stored at —20°C., were used. Particular 
care was taken to make sure that the anti-A gave strong agglutination with A, and 
A,B bloods. All grouping was done in test tubes and usually in batches of 60 to 
120 bloods. Cells and serum were grouped ; a different technician reading the 
agglutinogen and the agglutinin content. Every batch of grouping was controlled 
with known cells and serum. Any blood showing a discrepancy at first reading 
between cells and serum was re-grouped. Any batch of groupings showing more 
than 5 discrepancies in 100 between agglutinogens and agglutinins at first reading 
was repeated in toto. (Such occasions were very rare, being usually due either to 
agglutinins acting as hemolysins in warm weather, or to attempting to read the 
results too quickly om days when the work was rushed.) 

Such findings as cold auto-agglutination, rouleaux formation, ABa and Aaf 
were fully investigated until the correct interpretation or subgroup was assigned 
(A,Ba,, A,a,, etc.). All Oa bloods were rechecked with pooled group O (a8) 
serum to eliminate the possibility of weak As being missed by the anti-A in the 
original test. Such bloods were very rare indeed. The recorded results were 
checked twice to avoid clerical errors. 

Rh grouping was by saline anti-D serum incubated for one hour at 37° C. and 
read by observing the pattern on the bottom of the tube with a hand lens. All 
bloods showing “‘ Rh-negative ” on this first reading were re-grouped using a different 
sample of saline anti-D. 


5. THE BLOOD GROUP DATA 
The ABO and Rh blood group data for men and women in each 
of the twenty-six counties of the Republic is given in the appendix. 
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As only a small number of these donors were born in Northern Ireland 
our data for these counties has been summed and entered under the 
heading ‘‘ Six Counties ”’. 


6. THE ANALYSIS OF THE ABO BLOOD GROUP DATA 


The reasonableness of the county samples was tested by the usual 
method of comparing the observed with the expected number of 

















TABLE 1 
| 
Gene frequencies 
Total 
County number a ia ——— alt; Anialedatcigaiedl 
of donors O A B 
| 
Carlow . ; : 688 | 21 22°25 0°05 | 0°714 | 0°209 | 0:077 
Cavan . = : 944 26 25°23 0-02 | 0:747 | 0°177 | 0:076 
Clare . : : 187 4 3°01 nee 0807 | 0°132 | 0-061 
Cork . : j 533 10 12°22 0-32 | 0°776 | 0°146 | 0-078 
Donegal : ; 142 I 5:01 2°39 | 0°726 | 0°177 | 0:097 
Dublin . : | 8,182 253 228-09 2°02 | 0729 | 0-201 | 0-069 
Galway . ; Avi 765 17 18°95 0-16 | 0:764 | 0°157 | 0:079 
Kerry . ‘ : 258 3 5°63 0-97 | 0'769 | 0 144 | 0-087 
Kildare . ; : 581 15 14°88 0:00 | 0°738 | 0-198 | 0-065 
Kilkenny ; ‘ 931 27 29°21 0-12 | 0°722 | 0-200 | 0:078 
Leitrim . ‘ ‘ 113 i 2°43 ees 0°768 | 0:163 | 0-069 
lex . : : 674 17 | 14°76 0:27 | 0°76Q | 0:164 | 0-067 
Limerick . ; 271 5 8°54 10g =| 0°729 | 0°188 | 0-083 
Longford , 2 236 | 8 5°95 0°55 =| 0°765 | 0-150 | 0:085 
Louth . ; : 1,650 31 39°58 142 | 0°750 | 0°186 | 0-064 
Mayo . : . 426 9 8-29 0:05 | 0°786 | 0148 | 0-066 
Meath . i ; 241 2 5°47 1°65 | 0°738 | 0-209 | 0:054 
Monaghan . . 160 3 5°50 0-86 | 0°743 | 0°155 | o-102 
Offaly . : . 156 5 5°41 0-02 | 0:706 | 0:212 | 0-081 
Roscommon . : 552 15 19°27 0-70 | 0:726 | 0:176 | 0-099 
Sligo. : . 157 6 3°77 a 0°750 | 0:184 | 0:066 
Tipperary. : 421 II 10°12 0:06 | 0°754 | 0°179 | 0°066 
Waterford. : 1,209 44 34°24 2:07. | 0-731 | 0°196 | 0°073 
Westmeath . ; 447 9 g'i2 0-00 | 0°763 | 0-171 | 0-066 
Wexford : , 438 II 16-99 1°49 | 0°689 | 0°226 | 0-085 
Wicklow ‘ . 1,037 27 30°95 0°37 | 0°717 | 0°213 | 0-070 
Six Counties . , 495 5 11°89 3°03 | 0°748 | 0°189 | 0-063 
Total data. : 21,894 592 587-48 0:03 | 0°738 | 0-190 | 0-072 





























ABs. For these calculations the numbers of men and women, Rh- 
positive and Rh-negative, were added together. The resulting values 
of x? and the frequencies of the O, A and B genes are presented in 
table 1. 

The good agreement between the observed and expected numbers 
of ABs suggests that the county data can be assumed to be good 
random samples. An interesting feature of this data is that, whereas 
all previously published data for Dublin showed a deficiency of AB 
donors (Dawson, 1952), the present Dublin data shows good agree- 
ment between the expected and observed numbers. This may perhaps 
be taken as an indication that there was a slight, but significant, bias 
in the previous data which has not appeared in the present survey 
due to the greater care that has been taken to ensure random sampling. 
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The variation in the distribution of group B was investigated by 
calculating the heterogeneity x? for the relative numbers of O+A 
and B+AB donors in the various counties. As its value is 34-98 for 
26 degrees of freedom (P =0-1) no significant heterogeneity is 
established. The county frequencies of gene B are plotted in fig. 1. 
There is no evidence of any area of Ireland having a significantly 
higher B gene frequency than any other area—the counties appear 
to be randomly distributed with respect to their B gene frequencies. 
(Gene frequencies calculated from Hart’s survey of blood groups in 
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Fic. 2.—Percentage frequency of group A gene. 














certain areas of Northern Ireland (Hart, 1944) are included in this 
map and in figs. 2 and 3. Hart’s data was sorted according to 
** present residence ” of donors.) 

The calculation of the heterogeneity y? for the relative numbers 
of A+-AB and O-+B donors will give an indication of the significance 
or otherwise of the variations from county to county of the frequency 
of group A. The x? has a value of 100-81 for 26 degrees of freedom 
which leaves no room to doubt the significance of the heterogeneity 
in the distribution of group A. As the frequency of group B is so 
small and shows no significant heterogeneity in its distribution, the 
variation in group A is to be interpreted as variation in the relative 
frequencies of genes A and O. These frequencies are plotted in 
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Fic. 3.—Percentage frequency of group O gene. 


























Fic. 4.—Percentage frequency of Rhesus-negative donors. 
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figs. 2 and 3. These maps clearly reveal an important difference 
between the eastern and western counties in that the counties in the 
west have higher O and lower A gene frequencies than those in the 
east. The frequencies are intermediate in the central counties. The 
general pattern is so clear that it would appear to be unnecessary to 
establish it by further statistical analysis. 


7. THE ANALYSIS OF THE RHESUS BLOOD GROUP DATA 


The analysis of the relative numbers of Rh-positive (D-positive) 
and Rh-negative (D-negative) donors (all ABO groups, men and 
women, summed) reveals significant heterogeneity between counties. 
The heterogeneity x? for 26 degrees of freedom is 44:01, P = 0-02. 
The frequency of Rh-negative donors in each county is shown in 
fig. 4, and it is apparent that almost all the counties with a high 
frequency form a block in the central region of the eastern half of 
the country. This block of counties is called Group 1 and includes 
Carlow, Cavan, Dublin, Kildare, Kilkenny, Louth, Meath, Monaghan, 
Westmeath and Wicklow. Group 2 includes all the other counties. 
This division separates those counties with a proportion of Rh-negative 
donors greater than 0-155 (Group 1) from those with a proportion 
less than 0-155 (Group 2) ; the only exceptions being Meath (0-149) 
in Group 1 and Sligo (0-166) and Kerry (0-182) in Group 2. These 
three exceptions fall into that third of the counties for which one 
has the fewest observations. 

When an area is subdivided by direct inspection of the blood 
group frequencies, precise statistical analysis of the heterogeneity of 
these frequencies between and within the subdivisions is logically 
hazardous. However, once heterogeneity is established in the data, 
it is justifiable to attempt to interpret this in terms of a spatial pattern 
of blood group frequencies. From the present data we can probably 
go no further than saying that it appears that in the central eastern 
counties of Ireland (Group 1) there is a significantly higher frequency 
of Rh-negative donors than in the other counties. 

The Rh-positive and Rh-negative frequencies in the total data of 
21,894 donors from all counties are 0-8387 and o-1613 respectively. 

The frequencies of Rh-positives and Rh-negatives in the 8,182 
donors born in Dublin are 0-8293 and 0-1707. These agree well 
with the previously published frequencies of 0-8376 and o- 1624 based 
on 4,058 women resident in Dublin (Stewart, 1947) x? = 1°35, 
P = 0°25. 


8. THE RELATIONSHIPS BETWEEN SEX, ABO AND 
RHESUS (Dd) BLOOD GROUPS 


The relationship between sex and Rh-positive and Rh-negative 
grouping was investigated by a 2X2 contingency table for the total 
data. The x? obtained is 1-02, 1 degree of freedom, P = 0-3. Thus 
there is no significant association between sex and Rh group. 
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The corresponding value of x? for the relation between the Rh 
and ABO groups is 5:10, 2 degrees of freedom, P = 0:08. In this 
analysis the numbers of donors of groups B and AB were added, thus 
giving a 3 <2 contingency table for the total data. The correlation 
between these blood groups is doubtfully significant and it is interesting 
to notice that what association there is can be wholly attributed to 
the data for Dublin, where group A contains a lower proportion of 
Rh-negative than do the other groups (table 2). The 3 x 2 contingency 
table for this Dublin data gives a x? of 12:24, 2 degrees of freedom, 
P = 0-004. 

The Dublin data is consistent with that for the rest of the country 
in showing good agreement between the expected and the observed 
numbers of AB donors and in the absence of any associations between 
the ABO blood groups and sex and between the Rh grouping and 


TABLE 2 
Proportions of O, A and B+ AB donors that are Rh-negative 





; | Total data 
Dublin Total data excluding Dublin 














° ae ‘ 0°1822 | 0:1648 0°1549 

A x ‘ : 0°1503 O°1541 0°1565 
| 

B+AB_. : 0'1762 | 0'1640 0*1569 





sex. It may be that this abnormal association of Rh-positive and 
group A is a rare chance error. Whether this is so, or not, will be 
decided by collecting more data for Dublin. 

The analysis of the proportions of O, A and B+-AB donors among 
men and women in the total data gives a x? of 3:03, 2 degrees of 
freedom, P = 0-2. Thus there appears to be no significant association 
between sex and the ABO groups. 

From the published data from many countries it appears that it 
is exceptional to find any associations between sex, Rh and ABO 
blood groups. The absence of significant correlations in the present 
survey may thus be taken to confirm, in a general way, the validity of 
the methods that were used for collecting the data. 


9. DISCUSSION 


For Europe the broad pattern of ABO blood group frequencies 
is now well known. The main western peoples have a high group A 
frequency while the frequency of the much rarer group B is highest 
among the Slavs and other eastern Europeans. 

The highest group O incidence is in the peripheral north-west— 
Iceland, Scotland, Ireland and parts of Wales. There are other 
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peoples living on the southern edge of the continent, in Sardinia, 
Crete and the western Caucasus, who, together with the Berbers of 
North Africa, have also high group O frequencies. 

The frequency of the Rh groups shows little variation as far as 
D grouping is concerned, but a notable exception to this is found in 
the Basque people who have a very high proportion of D negatives, 
along with much group O and very little group B. 

Within the British Isles group A is common in the south of 
England (the gene frequency is between 0-25 and 0-30) but the 
A:O ratio decreases as one travels north into Scotland. In the 
Highlands of Scotland, a region in which there are strong Norse and 
Celtic elements, the A gene frequency falls to below o-20. The 
ABO frequencies in Wales have been shown by Mourant and Watkin 
(1952) to contain remarkable local divergences. There is a high 
group O area (gene frequency of over 0-75) in North Wales, high 
group B (over 0-15) in a small area in Carmarthenshire, and a high 
group A area (over 0:30) in the south of the Principality. 

Against this general pattern of blood group frequencies in Europe 
we can consider the interpretation of the results of the present survey 
of Ireland. 

In Ulster the effect of the vast English and Scottish plantations 
under James I is clearly reflected in the relatively high A gene 
frequency. In a comparative analysis of the blood group frequencies 
of people with English, Scottish and Irish surnames in Northern 
Ireland (Hart, 1944) it was shown that while the frequencies among 
Scottish and Irish were similar those among English tended to match 
the frequencies found in England. Thus the relatively high A gene 
frequency in the north is probably to be attributed primarily to 
English settlers. It is interesting to note that the A gene frequency 
found in the present small sample from Northern Ireland is somewhat 
lower than the frequencies found by Hart. This is probably due to 
the fact that the present survey includes only those born in the Six 
Counties and now resident in the south—they would be expected to 
be drawn primarily from among those who had connections with 
southern Ireland, where lower A gene frequencies are found. 

In the three southern provinces of Ireland the most striking feature 
of the county pattern of the ABO blood groups is the highly significant 
difference in the O and A gene frequencies between the eastern and 
western counties. This difference was first shown by Hooper (1947). 
The O gene occurs with highest frequency (0-76 to 0-81) and the A 
gene with lowest (0-13 to 0-16) in the western counties of Mayo, 
Galway, Clare, Kerry and Cork. The eastern counties with O gene 
frequencies between 0-69 and 0-74 and A gene frequencies between 
0-19 and 0-23 approximate closest to the frequencies found in England. 
In a large sample drawn from many parts of the United Kingdom the 
O and A gene frequencies were 0-68 and 0-26 (Dobson and Ikin, 
1946). These regions rich in A genes correspond almost exactly with 
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the areas of most intensive Anglo-Norman settlement, later supple- 
mented by settlements under Cromwell and to a much lesser degree 
under James I. 

It is, however, interesting that the A gene frequency in the eastern 
counties is as high as it is. While admitting the unreliability of 
deducing gene frequencies in past populations from present-day gene 
frequencies, we may tentatively accept Dobson and Ikin’s frequency 
of the O gene as being roughly its frequency among those coming to 
Ireland from England in historical times. If the original population 
of the whole of Ireland had an O gene frequency similar to that now 
found in the west (0-77) we should have to conclude that the present 
population of the eastern counties (0-73) was a mixture of Anglo- 
Norman and later settlers with original inhabitants in about the ratio 
of 4 to 5. This is in reasonable accord with much historical evidence, 
and in particular with the fact that, following the Anglo-Norman 
incursions, large numbers of merchants and farmers from England 
settled in the eastern counties of Ireland. 

The influence of the Norse settlements appears to be masked by 
the east-west gradient with the possible exception of Linierick, where 
the relatively high A gene frequency may be due to a Norse component 
in the population. 

If this historical interpretation is valid we must conclude that the 
present pattern of A, B and O gene frequencies in Ireland is primarily 
the result of invasions, incursions and settlements that have occurred 
during the last 1200 years. 

The pattern of the A and O gene frequencies is also apparent in 
much anthropometric data for Ireland. For comparison some of 
this data is reproduced from maps in Coon’s The Races of Europe, as 
fig. 5. For all these measurements those of the east of Ireland are 
more close to those of England than are those of the west. It is also 
significant that those regions of Ireland where Irish is most widely 
spoken are situated in the zone of highest O gene frequency. 

The present data amply confirms the generally accepted picture 
of a slightly but significantly higher B gene frequency in the extreme 
west of Europe than in those areas immediately adjoining. The data 
of Dobson and Ikin (1946) for the United Kingdom showed a frequency 
of 0-060, that of Fraser Roberts (1953) for the north of England 
0-062 while the present data gives an average frequency of 0-072 for 
the whole of Ireland. 

The present slight variation in the distribution of the Rhesus (D) 
blood groups may reflect a prehistoric pattern. 

In a study comparing the distribution of the TH sound among 
the languages and dialects of the peoples of Europe with the distribution 
of the frequency of the O gene, Darlington (1947) stated that the 
Celtic languages contain the TH sound. While this may have been 
true in the distant past, modern Irish possesses the dental plosive 
and not the fricative. The only exceptional region is south-west 
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Ireland where the alveolar plosive is found and Darlington, noting 
the absence of the TH sound from this region, suggested that it and 
Portugal (a non-TH area) may represent relic regions of a Megalithic 
people. However, it now appears impossible to link these regions 
together as Portugal has one of the lowest O gene frequencies in 
western Europe (0-60) while the frequency in south-west Ireland is 
higher than in any other region of Europe (0-77). 

In general it may be said that the gradient of blood group gene 
frequencies found in Ireland and the absence of obvious pockets of 
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Fic. 5.—The distribution of physical characteristics. 
(Taken from maps in Coon’s The Races of Europe.) 


exceptional frequencies reflect the absence of geographical barriers 
to population movements within the country. In this the type of 
gene frequency pattern is similar to that of England and in marked 
contrast to that of Wales. 

The present data confirms the similarity in O gene frequencies 
between western Ireland and Iceland, where its frequency is 0-75. 
This suggests a closer genetic relationship between the populations of 
these two areas than between Iceland and modern Scandinavia, 
which, in spite of cultural affinities, have widely different O gene 
frequencies. 

F 
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10. SUMMARY 

1. Records of 21,894 blood donors in Ireland were sorted according 
to the counties in which the donors were born. The analysis of the 
ABO gene frequencies has shown that while the B frequency is higher 
than in England and varies little, there is significant reciprocal 
variation in the frequencies of the O and A genes. The average 
frequency of the O gene in those counties in the west and south-west 
is 0:77 while in the eastern counties the frequency is only slightly 
greater than 0-70. 

2. In the three southern provinces the regions of highest A gene 
frequency correspond to those regions which have received most 
settlers from England—Anglo-Norman and Cromwellian settlements 
being the most important. In Ulster the high A gene frequency is 
primarily attributed to the English who entered the area under the 
scheme of plantation of James I. The distribution of the O and A 
gene frequencies is thus interpreted in terms of additions to the 
population of Ireland during the last 1200 years. 

3. The patterns of various anthropological measurements are in 
general agreement with the patterns of O and A gene frequencies. 

4. Doubtfully significant heterogeneity between counties in the 
frequency of Rh-positive donors is recorded. There may be a block 
of counties in the central eastern part of the country which has a 
higher percentage of Rh-negatives than elsewhere. 
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APPENDIX 
TABLE 1 
Total data (Men) 
Oo A AB 
County Totals 
Rh+ | Rh— | Rh+ | Rh— | Rh+ | Rh— | Rh+ | Rh— 

Carlow 212 42 120 39 14 12 5 488 
Cavan 356 69 | 167 53 82 15 II 6 759 
Clare . 81 12 28 3 13 I 3 to) 141 
Cork . 218 38 81 14 52 3 9 I 416 
Donegal 56 7 26 10 17 2 I te) 119 
Dublin 2395 | 537 | 1568 | 274 | 473 | 100 | 151 | 29 | 5527 
Galway 2 44 | 143 23 59 8 12 2 587 
Kerry 9 20 39 II 19 3 2 I 194 
Kildare 218 55 | 135 22 47 4 10 2 493 
Kilkenny 305 69! 189 37 68 II 24 o 703 
Leitrim 4 19 I 9 2 6 I & 
Leix 269 42 128 15 50 10 6 3 523 
Limerick 106 17 53 II 22 3 4 ro) 216 
Longford 95 II 38 7 23 3 6 fe) 183 
Louth 605 | 129 | 325 56 96 23 18 3 1255 
Mayo 160 24 66 10 22 8 7 to) 297 
Meath go 17 64 9 17 2 I I 201 
Monaghan . 56 10 30 6 15 2 2 I 122 
Offaly | 84 8 42 6 16 2 3 2 133 
Roscommon | 171 30 100 17 50 9 6 3 386 
Sligo . 49 12 29 6 II 2 2 I 112 
Tipperary | 161 27 96 10 30 7 6 2 339 
Waterford | 427 77.| 261 32 gI 12 33 5 938 
Westmeath . | 164 35 93 12 23 10 7 ra) 344 
Wexford | 157 27 | 107 18 42 5 9 I 306 
Wicklow |} 304 55 | 218 50 69 6 9 4 715 
Six Counties | 213 go | 112 21 34 5 3 oO 418 

Totals . 7361 | 1448 | 4277 | 773 | 1494| 272 | 363 73 | 16,061 
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TABLE 2 
Total data (Women) 























O A B AB 
County aa Sees ——n| OCD 
Rh+ | Rh— | Rh+ | Rh— | Rh+ | Rh— | Rh+ | Rh— 
Carlow ; . 87 10 68 9 20 2 3 I 200 
Cavan ; . 76 26 49 10 10 5 8 I 185 
Clare . ‘ ‘ 26 2 | II I 3 2 I oO 46 
Cork . ‘ ; 53 13 31 i 12 I fc) te) 117 
Donegal . . 14 oO | 4 2 2 I o o 23 
Dublin ‘ . | 1152 | 253 749 | 136 | 230 62 67 6 2655 
Galway ‘ : go 18 30 7 27 3 3 te) 178 
Kerry ‘ ‘ 35 5 9 5 8 2 oO 0 64 
Kildare ; ; 35 8 30 5 5 2 2 I 88 
Kilkenny. ; 95 39 69 12 21 II o 3 228 
Leitrim : , 14 2 4 6 I o o 3) 27 
Leix . ‘ ‘ 75 II 32 13 II I 7 I 151 
Limerick . . 18 5 18 3 9 I I oO 55 
Longford . ; 23 8 13 I 6 oO 2 fe) 53 
Louth ‘ ; 175 23 114 26 38 9 8 2 395 
Mayo . ‘ 73 6 28 4 15 I 2 oO 129 
Meath ‘ ‘ 21 5 II 2 I o oO o 4o 
Monaghan . . 20 4 5 oO 7 2 oO fe) 38 
Offaly ” , 13 3 6 o I o o oO 23 
Roscommon . 83 9 37 5 23 3 4 2 166 
Sligo . . , 2 2 10 3 3 oO 3 oO 45 
Tipperary . ‘ 44 7 18 3 4 3 2 I 82 
Waterford . ‘ 110 27 84 13 24 7 5 I 271 
Westmeath . : 54 7 19 6 14 I I I 103 
Wexford. : 24 3 32 4 7 I I to) 72 
Wicklow. ‘ 140 36 89 9 30 4 9 5 322 
Six Counties ° 30 8 23 4 9 I I I 77 
Totals . . | 2604 | 518 | 1593 | 2906 | 541 | 125 | 130 26 5833 


















































GENOME AND POLYPLOIDY IN THE GENUS TRILLIUM. 
VI. HYBRIDISATION AND SPECIATION BY CHROMO- 
SOME DOUBLING IN NATURE * 
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Received 1.V.55 
1. INTRODUCTION 


THERE are two forms in Trillium Hagae Miyabe et Tatewaki, the 
sterile triploid and the chromosome-doubled fertile hexaploid. These 
two forms were concluded from chromosome evidence to be the 
hybrids between diploid Trillium kamtschaticum Pallas and tetraploid 
T. Tschonoskii Maximowicz (Haga, 1937, 1951@, 1952; Haga and 
Kurabayashi, 1953). The triploids were found at almost all the 
localities where both parental types are growing close together. The 
chromosome constitution was revealed to be that expected of the F, 
without exception. This is obviously due to the frequent occurrence 
of natural hybridisation and the nearly complete defectiveness of the 
triploids in seed setting. On the other hand, the chromosome-doubled 
forms are completely fertile in seed setting. These facts suggest that 
hybridisation and chromosome doubling are at present occurring in 
nature. In other words, it is likely that a new hexaploid species is arising 
following the processes by which widespread older polyploid species 
of the genus originated. 

The chromosomes of 7. kamtschaticum revealed profound structural 
diversification and heterozygosity in the patterns of the hetero- 
chromatic segments after cold-treatment (Haga and Kurabayashi, 
1954). This finding suggests that chromosome variation might be 
present in the chromosome-doubled forms at the time of their 
emergence. This was proved to be the case, as has already been 
briefly reported (Haga, 19515). This article attempts to relate these 
observations with the problem of speciation in the genus. 

Cytological observations were carried out using the preparations 
of cold-treated chromosomes in young ovular tissues, using the method 
of Darlington and La Cour, 1940, with some modification in the 
method of cooling (Haga and Kurabayashi, 1953). This method 
differentiates the eu- and hetero-chromatic segments which are 
constant in number, size and position for a given chromosome, as 
was observed by Darlington and La Cour (1938, 1940, 1941), La 
Cour (1951), Haga (1944), Haga and Kurabayashi (1953, 1954). 

Some minor controversies in respect of the constancy of the 
patterns are due in my opinion to the varying degrees of differentiation. 

* Contributions from the Department of Biology, Faculty of Science, Kyushu University 
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This may range from the non-differentiated to the fully-differentiated 
state even in cells of the same cold-treated tissue (Kurabayashi, 1948, 
1952), a fundamental fact which appears not to have been clear in 
some other studies (¢f. Wilson and Boothroyd, 1940, 1944). Observa- 
tions reported in this article were based on patterns observed in the 
final stage of differentiation. It should be noted that the hetero- 
chromatic segments are seen only after the above treatment. 

All the text-figures are the classified alignments of chromosomes 
in single cells at metaphase. 
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Fic. 1.—Range of variation in the patterns of the cold-treated chromosomes of T. 
kamtschaticum, heterochromatic segments being shown in dotted line or in slit-like gap, 
euchromatic segments in solid black, and centromere diamond-shaped. Top, 8 types 
of chromosome D. Bottom, 8 types of chromosome E. X1500. (For the other 
chromosomes, A to C, see Haga and Kurabayashi, 1954.) 





2. CHROMOSOME CONSTITUTION OF THE PARENTS 


(a) ZT. kamtschaticum. Five chromosomes, A to E, constituting 
a gametic set of this species, revealed many types in respect of the 
patterns of the heterochromatic segments (cf. fig. 1). The types were 
indicated by numbering as, for example, A1, A2 and so forth in 
chromosome A. ‘The homozygous presence, for example, of Az is 
shown as A-2-2 and a pair heterozygous for A1 and Ag as A-1-2. 
This method of description was adopted for all the five chromosomes, 
the numbering of the types being the same as that used in the previous 
paper (Haga and Kurabayashi, 1954). 

Plants from homogeneous populations show one and the same 
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constitution, that is all five pairs are homozygous (fig. 2), whereas 
those from highly heterogeneous populations show constitutions 
differing from plant to plant, the plants being heterozygous for one 
to all five pairs (fig. 3). 
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Fic. 2.—T. kamtschaticum (2x) from Kusakai. 
Homozygous plant, Ks-49-2. X 1500. 


(b) 7. Tschonoskii. As to the patterns of the heterochromatic 
segments of the tetraploid parent only a single plant from Siraoi was 
studied and it was found to be homozygous for two sets of chromosomes 
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Fic. 3.—T. kamtschaticum (2x) from Siraoi. Plant heterozygous for 
A, B and D, Sr-49-7. X 1500. 


i? 
es 


which were distinguishable by the patterns of differentiation (fig. 4). 
As manifested in the meiotic chromosome pairing in the triploid 
hybrid, one of these two sets, K,, is intimately homologous with set 
K, of T. kamtschaticum, though these two sets are not morphologically 
identical in untreated as well as in cold-treated chromosomes. 
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The other set, T, of the tetraploid parent has also relations to the 
sets K, and K,, though the relation is distinctly not intimate in 
comparison with that between K, and K, (Haga, 1937, 1951@). 
However, it is impossible to sort out individually the members com- 
posing K, and T from the combined set of 20 chromosomes. Arbitrarily 
five pairs having a large number of heterochromatic segments were 
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Fic. 4.—T. Tschonoskii (4x) from Siraoi, homozygous allotetraploid, K,K, and TT being 
sorted out above and below. Chromosomes A to E are aligned from left to right in 
each of the four sets. X 1500. 


sorted out and designated as K,, and the remaining five pairs as T. 
Thus assortment of chromosomes by shape at mitosis into K, and T 
is less clear than the distinction manifested by the meiotic pairing. 
The patterns of differentiation of K, and T comprised in 59 triploid 
hybrids and in 2 chromosome-doubled ones were the same as one 
another and as those of 7. Tschonoski (figs. 4, 5 and 7). Thus there 
is a high stability of the linear chromosome organisation of the 
tetraploid species. 





-\ 


Fic. 5.—T. Hagae (3x) from Kusakai. Plant H-Ks-50-2. Top, K,, 
middle K,, bottom T. X 1500. 
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Fic. 6.—K, from 3 triploid hybrids, K, and T which are the same as shown in fig. 5 are 
omitted. Top, plant H-Sr-48-1 from Siraoi. Middle, plant H-Ut-48-3 from Utoma. 
Bottom, plant H-Ty-48-1 from Toyoni. X 1500. 
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3. CHROMOSOME CONSTITUTION OF THE HYBRIDS 
As mentioned above, two sets of chromosomes K, and T in the 
hybrids revealed no variation. On the other hand, K, from T. 
kamischaticum showed variations in the patterns as was to be predicted 
TABLE 1 
Chromosome constitution of Ky, in 59 triploid hybrids 
in nine populations 


(Figures under the chromosome symbol indicate the number of 
the chromosome types) 
























































Locality A B Cc D E Freq. Total 
1. HOMOGENEOUS 
Kusakai_. : 13 12 6 3 7 6 6 
Hagino : - ig | 2 6 3 2 mu Ir 
Shadai ; : 13 7 6 9 2 7 i] 
2. HETEROGENEOUS | 
Siraoi ° , 13 6 6 3 2 ? 
13 6 6 7 2 I 
14 6 6 7 2 7 9 
Utoma 5 2 3 5 I I 
SI - | 3 5 3 7 
s. 4 2 | 3 5 4 I 
Pie | F 5 3 I 
7 & 1 9 5 4 I 
6 eo. 1 & 5 5 I 
6 | ey 36 5 3 I 
mo | -m | 5 7 3 I 
17 | 2 | 3 5 3 I 9 
Minami-Sat. , 10 | 2 3 5 I r 
10 | 2 3 5 2 I 
10 | 3 5 5 3 I 
16 | 2 3 5 3 Z 
6 | 38 5 5 4 I 5 
Toyoni 5 I 3 5 3 I 
5 | 2 3 5 8 I 
5 2 8 5 3 I 
| 8 3 5 3 r 
| 10 | r | 3 5 3 2 
14 | “2 5 8 I 
i SS 5 8 I 
21 | I 3 5 8 I 
23. | 2 3 5 I I 10 
3. SINGLE SPECIMENS | 
Samani_. : 2 | 2 6 6 3 I I 
Shoya : ; 8 | 3 | 3 5 3 I I 
| 














from the intraspecific chromosomal variations in 7. kamtschaticum 
(figs. 5 and 6). 

Chromosome constitutions of K, in 59 triploid hybrids are presented 
in table 1, the constitutions of K, and T being omitted. At Kusakai, 
all 50 plants of 7. kamtschaticum showed one and the same homozygous 
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Fic, 7.—T. Hagae (6x) homozygous allohexaploid from Sizunai. X 1500. 
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constitution ; A-13:13, B-12-12, C-6-6, D-3:3, and E-7-7 (fig. 2). 
The constitutions of K, of all six hybrids from this locality were 
revealed to be the same as the gametic set of 7. kamtschaticum at this 
locality ; A-13, B-12, C-6, D-3 and E-7 (fig. 5). This is a natural 
consequence of the complete homozygosity of the diploid parent 
population at this locality. On the other hand, at localities where 


TABLE 2 


Frequencies in T. kamtschaticum populations of the chromosome types 
comprised in the triploid hybrids 















































Population * Chromosome type and its actual frequency (_) 
Kusakai . - | Arg (100) | Br2 (100) C6 (100) Dg (100) E7 (100) 
—- cs See ts fe en ee 
Siraoi. - | Arg (44) B6 (60) C6 (100) D3 (27) E2 (100) 
Ar4 (56) | Bz (18) ae D7 (73) ve 
Nisikioka ** . | Ag8 (30) Bez (30) C6 (30) D7 (30) E3 (30) 
Samani . . | Ag (14) | Be (63) Ce (8) | Ds (26) Er (9) 
A5 (19) Bg (33) Cg (31) | D6 (39) E3 (67) 
A6 (7) | Brq (0) C5 (22) D7 (14) E4 (17) 
A8 (3) | C6 (16) Es (1) 
A16 (12) C7 (2) 
Ar7 (2) Cg (1) 
Toyoni ** ~ | Ags (11) | Br (23) | Cg (155) D5 (34) Er (12) 
Aro (4) | Be(100) | C8 (o) sxe E2 (0) 
Aig (1) Bg (59) | ee me E8 (5) 
A18 (8) | 
Aart (5) ‘ | Er (12) 
A323 (10) | E2 (0) 
vk aoa ks me ee = 
Minami-Sat ** | Aro (6) | Be (23) — C3 (33) | Ds (34) Eg (20) 
| Ar6 (4) —|_-Bg (13) | C5 (1) ee E4 (13) 








* Shadai is located close to Nisikioka, and Samani to Utoma. The number of plants 
observed is as follows: Kusakai = 50, Siraoi = 50, Nisikioka = 15, Samani = 50, 
Toyoni = 86, and Minami-Satunai = 36. 

The number of chromosomes observed for each of the chromosomes, A to E, is twice 
the number of plants observed. (ex Haga and Kurabayashi, 1954, and unpublished 
data (**) accumulated by the Trillium Research Group of the Hokkaido University). 


T. kamtschaticum populations are large, dense and highly heterozygous, 
the constitution of K, in the hybrids varies from plant to plant even 
within a population, that is, at Utoma, Minami-Satunai and Toyoni 
(fig. 6 and table 1). Thus it is evident that the hybrids with different 
K, constitutions were raised through matings of the 7. kamtschaticum 
gametes of various constitutions with those of T. Tschonoskii. Further, 
we can attribute the uniformity of the chromosome complements of 
the hybrids in certain populations to the complete or very high 
homozygosity of the 7. kamtschaticum population. Such cases are found 
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at Hagino and Shadai in addition to the case found at Kusakai 
(table 1). Generally speaking, K, comprised in the hybrids is found 
to consist of chromosomes which are frequent in the 7. kamtschaticum 
population at the locality where the hybrids were sampled (tables 1 
and 2). In other words, hybridisation takes place where both the 
parents are growing close together, and the hybrids grow up to 
maturity in the localities native to their parents (fig. 8). 





N 
HOKKAIDO 14 


100 KM 





LOCALITIES 





KUSAKAI 9 UTOMA 





| 
2 INO 10 SAMANI 
3 TANABU I) SHOYA 
4 HAGINO 12 TOYONI 
Ss AO 
HONSHU 6 simaon 
6 SHADAI 3 MINAMI- 
SATUNAI 
7 BIRATORI 
lo 8 SIZUNAI 14 AKKESI 











Fic. 8.—A map of the northern districts of Japan, showing localities where the hybrids 
T. Hagae (3x) were found growing, except at Ino (2) and Tanabu (3) where no hybrid 
was recorded. 

Chromosome constitutions of the doubled hybrids were analysed 

in two plants, one from Sizunai and the other from Utoma (fig. 8). 

These two plants were homozygous for all three sets K,, K, and T 

(fig. 7). However, these two doubled plants were different in respect 

to the constitution of K, as shown below. 

Sizunai A-5'°5 B-4-4 C-66 D-5:5 E-2:2 
Utoma A-5°5 B22 C66 D-5:5 E-3:3 

All the chromosome types were those found to be highly frequent 

in T. kamtschaticum populations at Sizunai and Samani, the latter 

locality being close to Utoma. At Utoma T. kamtschaticum and T. 

Tschonoskii grow in abundance but unfortunately these parental 
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plants at this locality have not been analysed cytologically as yet. 
The chromosome doubling manifested in the hexaploid may have 
occurred in one of the following three ways ; firstly, by the union 
of unreduced gametes produced in F;, sterile triploid hybrids, secondly, 
by the union of unreduced gametes produced separately in both 
parents, and thirdly, by chromosome-doubling in the F, fertilised 
egg cells. The first possibility seems nearly negligible as is evident 
in the meiotic chromosome behaviour in F, triploid hybrids (cf. Haga, 
1951; Haga and Takizawa, 1953). ‘The second alternative is not 
utterly impossible, since unreduced pollen grains are produced with a 
frequency as high as 17 per cent. in 7. kamtschaticum under extreme 
conditions (Haga and Kayano, in press). But it is not likely since 
no hybrid was of the composition K,K,K,T or K,K,K,TT, which, 
if the second process is really occurring, will be of more frequent 
occurrence than the hybrids K,K,K,K,TT. Thus chromosome- 
doubling in the fertilised F, egg cell seems the most probable. This 
process further explains why the two sets K,K, are homozygous in 
the doubled hybrids, despite the fact that they were sampled in the 
highly heterozygous 7. kamitschaticum populations. 


4. POLYPLOIDY 

Speciation by polyploidisation takes place by means of the following 
three steps in a serial succession. First, genetic differentiation of the 
genomes which are to be combined in a new polyploid species, second, 
hybridisation of parental plants having differentiated genomes, and, 
third, chromosome-doubling which makes the hybrids constant 
through generations. Three genomes, K,, K, and T, are those 
differentiated from each other, K, being intimately homologous with 
K,, and T homologous with both, though the latter relation is not 
intimate. The triploid hybrids, K,K,T, show meiotic pairing 5,,+5, 
and nearly complete sterility in seed setting, no F, or back-cross 
progeny being found in nature. As the chromosome-doubled hybrid, 
K,K,K,K,TT, forms 15,, by the differential affinity between the 
three genomes, the hybrids are perfectly fertile and capable of produc- 
ing offspring. 

Experimental crosses between the parental plants are easily 
attained irrespective of the direction of crossing, and as many good 
seeds are produced as are borne on the parent plants. Unfortunately, 
artificial germination of the seeds produced either by the hybrids or 
by the parental plants has not been successful. However, a method 
of germination by Barton (1944) was successful. 

In nature, F, hybrids are found at almost every locality where 
both parental species are found growing together (fig. 8). The 
frequency of these F, hybrids will depend on various ecological factors : 
on numerical proportions and densities of the two parental species, 
and time relation in flowering of the parental plants. For example, 
even within the species 7. kamtschaticum the average date of flowering 
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in the Trillium Reservatory Garden of the Hokkaido University, 
Sapporo, varied among the plants sampled from different localities. 
In 1953 it varied from May 4th to 5th in plants from Saruru (close 
to Shoya, cf. 11 in map) to May 13th to 14th in plants from Tomikawa 
close to Biratori, cf 7 in map) (¢f. Hiraizumi, 1954). Thus the 
frequency of hybridisation will vary at different localities and in 
different niches. 
TABLE 3 
Frequencies of parental plants and their hybrids in a summer-green forest 
at Siraoi (ex Haga 1942). 














Year 1937 1939 1940 1941 Total | Per cent. 
T. kamtschaticum ‘ 216 150 941 706 2016 565 
T. Hagae (3x) . . 34 18 74 43 169 4°7 
T. Tschonoskii . ‘ 224 185 494 481 1380 387 
Total . 3 474 349 1509 1233 3565 99°9 


























With the assistance of several students, collections of Trillium 
plants were carried out at Siraoi in successive years 1937-41 (except 
1938) with no preference for particular species or hybrids. The 
purpose of the collecting in these years was to stock the Trillium 
Reservatory Garden, therefore the records are not critical from an 
ecological point of view. However it will be helpful to see how 
frequently hybridisation is occurring at this locality. According to 
the records, the occurrence of F, hybrids is 4-7 per cent. of the total 
number of plants (table 3). 








TABLE 4 

Chromosome types of K, which have reached fixation 
Population Chromosome types Species n 
Kusakai_. Alg3 Bi2 C6 D3 E7 T. kamtschaticum | 50 
Alg Bi2 C6 D3 E7 T. Hagae (3x) 6 
Hagino , Alg3 Ba C6 D3 E2 T. Hagae (3x) II 
Siraoi ° ine ws C6 cas E2 T. kamtschaticum | 50 
TS B6 * C6 sus E2 T. Hagae (3x) 9 
Shadai .| Arg B7 C6 D7 E2 T. Hagae (3x) 7 
Nisikioka . A38 Ba C6 D7 E3 T. kamtschaticum | 15 
































* Cf. table 2, Siraoi. 


5. POPULATION STRUCTURE 
Chromosome pairs of 7. kamtschaticum are homozygous in small 
isolated populations, individual plants being composed of the same 
homozygous five pairs of chromosomes. This circumstance is mani- 
fested directly by the analysis of 7. kamtschaticum populations and 
indirectly by the surveys of the constitutions of K, in the hybrids. 
According to these lines of evidence four populations of 7. kamt- 


schaticum seemed to have reached fixation in all five chromosome pairs 
(table 4). 
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Here it is important to note that chromosome constitutions of 
homozygous or nearly homozygous 7. kamtschaticum populations are 
different from one another even though they are located close together, 
that is, populations at Shadai, Nisikioka, Siraoi and Hagino. Indeed, 
Shadai and Nisikioka are separated by only a few hundred metres ; 
Shadai is located about 4 kilometres east of Siraoi and Hagino is 
located about 8 kilometres west of the latter locality (fig. 8). Inter- 
vening areas between these localities are occupied by dense forests, 
by meadows, or by agricultural land where Trillium plants are now 
entirely absent or almost so. The fixation of a definite chromosome 
constitution at these localities is considered to have been brought 
about by the cooperation of genetic drift and severe selection in small 
populations (Haga and Kurabayashi, 1954). Then chromosome- 
doubled hybrids, raised in isolated 7. kamtschaticum populations such 
as are described above, would be different in their genetic constitutions 
from the very time of origin. 

Large 7. kamtschaticum populations, on the other hand, or those 
that have recently become small and discontinuous, are highly hetero- 
zygous for many chromosome types. In some extreme cases, every 
plant of a population differs from all the other members of the 
population in its chromosome constitution (Haga and Kurabayashi, 
1954). Itis not surprising therefore that the doubled hybrids appearing 
among such populations of 7. kamtschaticum vary from plant to plant. 
In other words, the genetic variations in the allopolyploid species 
originate from the genetic variations already present within, as well 
as between, the populations of the lower-ploid progenitors. Of 
course, the variation will proceed further and seems to accumulate 
from the time of origin of the polyploid species (cf. Haga, 1940 ; 
Stebbins, 1950). 


6. SUMMARY 


1. The genome of the diploid 7. kamitschaticum, K,, shows great 
variation. Many different types of chromosomes are distinguished 
in all five members, A to E, by the patterns of heterochromatic 
segments showing nucleic acid starvation. Populations of this species 
are either (i) completely homozygous for a certain definite type in 
all of the five chromosomes ; or (ii) heterozygous for many types in 
some or all of the five chromosomes. 

2. The two genomes of the tetraploid 7. Tschonoski, K, and T, 
were distinguishable in all five chromosome types but they revealed 
no heterozygosity or even variation between themselves which was 
studied in detail by the patterns of heterochromatic segments in the 
cold-treated chromosomes. 

g. Natural hybrids 7. Hagae (3x = 15, K,K,T and 6x = 30, 
K,K,K,K,TT) between T. kamtschaticum and T. Tschonoskii are 
found at nearly all the localities where the parental species are growing 
together. 
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4. Triploid hybrids from a homozygous T. kamtschaticum population 
possess genome K, identical with that of the homozygous diploid 
progenitor. Those from the heterozygous population show K, con- 
stitutions differing from plant to plant in accordance with the hetero- 
zygosity of the population of the diploid parent. 

5. Generally, the K, chromosome types of the triploid hybrids 
were the same ones found in the neighbouring 7. kamtschaticum 
populations. 

6. Both the chromosome-doubled hybrids sampled in two highly 
heterozygous 7. kamtschaticum populations proved to be homozygous 
for K,K,. But the K, constitutions of these two plants were different. 
Chromosome doubling has thus taken place after fertilisation in the 
F, and has been followed by inbreeding. In this way new constant 
hexaploid plants are arising at present in nature and the genetic 
variations may be involved in allopolyploid species from the time that 
they arise. 
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I. INTRODUCTION 


IN interchange heterozygotes the rings and chains of four chromosomes 
seen at first meiotic metaphase fall into two classes, arranged dis- 
junctionally and non-disjunctionally. Disjunctional arrangements 
appear as zig-zag configurations with successive centromeres directed 
to opposite poles thus, 


AB CD fig. 1). 
Naa 4 Nai (see fig. 1) 


This arrangement results in balanced viable gametes. Non-dis- 
junctional arrangements do not show this alternation of centromere 
orientation, and appear as open ring or chain configurations : ¢.g. 
AB—— --B BA AD 
or | 
AD DC BC. cD 


In plants these configurations result in the production of inviable 
gametes. 

Interchange heterozygotes in Pisum (Pellew and Sansome, 1931) 
and in maize (Burnham, 1934) give disjunctional and non-disjunctional 
arrangements with approximately equal frequencies and, as would 
then be expected, these plants show about 50 per cent. pollen and 
ovule sterility. In other species, however, plants heterozygous for 
an interchange show less than 50 per cent. sterility ; e.g. Campanula 
(Darlington, 1937) shows only 30 per cent. sterility, disjunction 
occurring in about 70 per cent. of cases. The higher disjunction 
of Campanula is accounted for by its larger chromosomes and by the 
high degree of terminalisation of its chiasmata (Darlington /.c.). 

Genotypic control of disjunction may be inferred from the differ- 
ences between species, but hitherto it has been analysed no further. 
The present account is concerned with differences in the disjunction 
properties of interchange configurations within a single species, 
differences which therefore afford an opportunity of further investiga- 
tion by breeding experiments. 











2. MATERIAL AND METHOD 


In 1951 an F, plant from a cross between two inbred lines of rye 
was found to be heterozygous for two independent interchanges. 
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These two interchange configurations were distinguishable at diplotene 
(fig. 1) as: 
Interchange A, involving the nucleolar organiser chromosome, and 
Interchange B, involving unequal lengths of chromosome, so that 
the ring or chain of four consisted of one short, one long, 
and two chromosomes of normal length. 
The subsequent history of the plants used is shown in fig. 2. In 





Fic. 1.—Meiosis in plants heterozygous for both of the interchanges. 


(a) Diplotene 

(6) MI showing two non-disjunctional rings of 4 

(c) MI showing non-disjunctional and disjunctional rings of 4 

(d) MI showing disjunctional ring and chain of 4. Two chromosomes are unpaired. 
X 1080. 


each generation seed was obtained by self pollination of plants hetero- 
zygous for both interchanges, each family being the progeny of a 
single plant. 

The F, plant was heterozygous for both interchanges. In 
subsequent generations, as a result of segregation, some plants were 
structural homozygotes, some were heterozygous for one of the inter- 
changes, and some were heterozygous for both. 

The interchange heterozygotes were classified into types A, B 
and AB at diplotene, and the number of disjunctional and non- 
disjunctional arrangements were then scored in 20 first metaphase 
cells of each of these plants. 














————————EEE 





Interchange A A and B B 
Class D N |No.of} DD | DN | NN |No.of| D N_ |No. of 
plants plants plants 
| Family a, a a, a, a a a, 
106 18 14 2 16 24 14 4 5 7 I 
107 102 31 8 30 23 9 4 63 27 5 
108 89 | 31 6 58 | 34 5 5 49 | 10 | 3 
Total 209 76 16 104 8: 28 13 117 44 | 9 
Proportion Po qa ve Pafp Peto aI» ee Py % | 
qaPr | | 
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3. RESULTS AND ANALYSES 
Fs 


The observed numbers of disjunctional and non-disjunctional 
arrangements in each F3 family are set out in table 1. 


F, ie er eit 
‘ ak, 
ae” Kae Gm” 




















F; 106 107 108 
tee —_—[—_—_ oe 
1 2 22 25 28 29 
Fy, 93 94 95 96 97 98 


Fic. 2.—The history of the material. P2 and P8 are inbred lines. Each number represents 
a family in F,, F,; and F, which in turn represents the progeny of a single plant 
self-pollinated in the previous generation. 


To test whether all of the F; families are behaving alike we can 
find estimates of the disjunction frequency, and compare these 
estimates with the observed values in a x? analysis. 


TABLE 1 


The observed numbers of disjunctional (D) and non-disjunctional (N) arrangements in the F 
families. (Each entry is the total from 20 cells scored in each plant) 
























































The most efficient method of finding the estimates is by maximum 
likelihood. If p, equals the disjunction frequency of interchange A, 
q,(= 1—p,) the frequency of non-disjunction, and g, and gq, the 
frequencies for interchange B, then the frequencies of the seven classes 
in table 1 are as shown in the last row of the table. 

G2 
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The log likelihood expression is 
L = a, log pa tag log 9, +43 log pap) +44 log(h.9+ Gabo) 
+s log 9295 +46 log p, +47 log 9, 
= (a, +4,)log p,+(42+45)log 9, +44 log(p.9,+9aho) 
+(a3+4¢)log p, + (a5+4,)log 9, 
Partial differentiation with respect to p, and p, gives the maximum 
likelihood equations 


dL _ +4,  a,+4s 4 4(9,—Po) 








win =o 
dp, Pa Va Palo + Vals 
dL a ag +g ™3 as -+ay + 44(7,—f,) =0 
dp, py Y Palo tIalo 


which reduce to 
2(4,;+a,+a3+4,+4;5)* p,° 
— (4, +@y +45 +44+45) (5a; +4, +6a3+34,+0,—a7)p,* 
+{(@, +43 +44) (40; +42 +543 +44 +44—47) + (a, +43) (42+45) 
+ (a@3-++4g) (@,+4@, +a3+45)}p, 
— (4,+4@3)(@,+24,+a4+4,.) = 0 

and 

bh = {(4, +4, +43 +44 +45)Pa— (a1 +43+44)}Pa 

2 (a, +a_+3+44+45)p,"— (3a; +42 +343 +204+05)p, +0, +45 





TABLE 2 


The maximum likelihood estimates of the disjunction and 
non-disjunction frequencies in the Fy 








Family Pa qa bo % 
106 0°5616 0°4384 0*4650 0°5350 
107 0°7461 0°2539 0°6744 0°3256 
108 0°7367 0°2633 08213 0°1787 
Estimate over 
all families 0°7136 02864 06944 0°3056 























The equation for p, can be solved by iteration, and that for p, by 
substitution of the value of p,. The estimates from the F, families, 
separately and pooled, are set out in table 2. 

The estimates of p, and p, derived from the F, totals are used to 
obtain expectations for the various classes in each of the three families. 
We can then calculate a x? testing the goodness of fit of the frequencies 
observed with the overall estimate of disjunction. This x? has 
ten degrees of freedom (see appendix) and will, of course, be influenced 
by any differences that may exist between the F; families, as well 
as by variation within them. 

Using the estimate obtained from each individual family a x? 
for two degrees of freedom can be found testing the goodness of fit 
within each of the three families. The total of these three x? is itself 








elf 





a 
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a x® for six degrees of freedom, and it will test whether departures 
from expectation within families are small enough to be attributed 
to sampling error. If significantly large it will reveal interaction 
between the interchanges within the family, i.e. whether one inter- 
change affects the disjunction of the other in plants heterozygous for 
both of them. 

The difference between the x?) and the x?, is a x? for four degrees 
of freedom which tests for heterogeneity of p values between the three 
families. 

The analysis of the F; is set out below. 














Item x? N ly 
Heterogeneity between families : 40°4118 4 -<0-001 
Interaction within families : ; 4°8332 6 0°50-0°70 
Total . * ‘ . 45°2450 10 

















It will be seen from this analysis that : 
1. There is no significant interaction effect within families. 
2. The families are highly heterogeneous: they show different 
frequencies of disjunction. 
The significant heterogeneity of F, families strongly suggests 
segregation at meiosis in the F, for genes controlling the disjunction 
of the interchanges. 


Fa 

As a further test of the hypothesis that the F, heterogeneity can be 
explained in terms of a genetic segregation six F, families were grown 
and analysed. These six families were grouped into three pairs 
according to their ancestry (see fig. 2). The observed numbers of 
disjunctional and non-disjunctional arrangements in the F, families 
are given in table 3. Table 4 shows the values of p, and , estimated 
from these observations. 

The analysis of the F, is made in the following way :— 


(1) The overall estimates are used in each of the six families to 
obtain a total x? for twenty degrees of freedom (see appendix). 
This x? tests overall goodness of fit. 

(2) Using the estimates from the paired family totals and calculat- 
ing expectations for each of the two families making up the 
pair, two x? for six degrees of freedom (families 93 and 94, 
and 95 and 96) and one x? for four degrees of freedom 
(families 97 and 98) can be obtained. These x? test for 
variation between the members of a pair of families and 
within the families. Together these three x? make up a 
x? for sixteen degrees of freedom. 
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(3) The estimates of p, and p, obtained from each individual 
family are used to obtain the interaction x? as in the analysis 
TABLE 3 


Numbers of disjunctional and non-disjunctional arrangements in the 
F, families (20 cells in each plant) 







































































| | | 
a A A and B B 
\. Class eee Px fer Oe = * 
Family \. p | N |No.of| DD| DN | NN |Noof| D | N. |Noof 
Wy @, | a, |plants| ay a% | Gs |plants| ag a, | plants 
ees: Sen: Wee cnet --— i NOW nad ce ial 
93 37 | 22 3 | 10 | 62 | 15 10 29 | 1 2 
o4 73 | 7). | oo) ey 3 24 | 16 2 
ee : ae - a ¥ | = — S| Eee el Seer | — 
ggandg4 | 110 | 69] 9 | 143 | 87 | 17 | #9 | 88] 97 4 
Pa pa re | | tit] 
95 95 | 19| 6 | o5 | 393/ 6 | 7 | 45] 8] 8 
96 61 | 17 4 | 36 | 16 3 3 41 18 3 
SeTRGL PN INOS IK lee EBS ee Go, 
95 and 96 156 | 36 | 10 | 131 49 9 10 86 | 26 6 
97 46 | 14 3 28 7 2 2 0 
98 65 | 14 4 67 | 28 2 5 o 
’ Senet TOON: DRL pees ae aici 
97 and 98 III | 28 | 7 | 95 | 35 4 7 0 
| 
seis | _ sina Ipeaeihchesuleia | eclaastageliadaes aieadipaiaeaialal —_ 
| | 
Total . - | 377 | 133 26 369 | 171 | 30 go 139 53 10 
' 
TABLE 4 


The maximum likelihood values of pg and p, estimated 
JSrom each row of table 3 








Family Pa Po 

93 0:6986 0°7716 

94 0°6478 0*7140 

95 0°8303 0*8401 

96 0*8009 0°7325 

97 0+7863 08845 

98 08210 0°8505 

93 and 94 0'6746 0*7604 
95 and 96 0*8236 0°7947 
97 and 98 08019 0°8737 
Over all families 0+7654 0°7852 

















of the F;. Four of these x? items have two degrees of 
freedom (families 93, 94, 95 and 96) and two of them have 
one degree of freedom (families 97 and 98). The total of 
these items is therefore a x? for ten degrees of freedom 
testing for interaction between the interchanges in double 
heterozygotes. 
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(4) The difference between (1) and (2) above is a x? for four 
degrees of freedom testing for heterogeneity of p values 
between the three family pairs (analogous to heterogeneity 
between the three F; families). 

(5) The difference between (2) and (3) above is a x? for six 
degrees of freedom testing for heterogeneity between 
members of the same family pair. (A test of genetic hetero- 
geneity within the F, families, z.e. of segregation at meiosis 
in the F, plants.) 


The results of this analysis are set out below. 








Item x? N P 
1 Grand total . : ‘ : 61-0249 20 
2 Family pairs total , A 26°1213 16 te 
1-2 Heterogeneity between family 34°9030 4 <0-001 
pairs 
3 Interaction within families ‘ 18-6678 10 | 0°05-0°02 
2-3 Heterogeneity within family pairs 7°4535 6 0°3-0°2 














It can be seen that :— 


(1) There is heterogeneity between the three pairs of families. 
This heterogeneity confirms that observed at F3. 

(2) There is no sign of heterogeneity between the families within 
pairs. 

(3) There is some doubt as to the presence of interaction within 
families. When the individual family x? are examined 
(see table 5) it can be seen that families 93 and 94 make 
up nearly the whole of the interaction item. 











TABLE 5 
The interaction x* of the F, families 

Family x? N P 
93 4°8735 2 0*1-0°05 
94 7°1378 | 2 005-002 
95 1-9865 | 2 ‘xa 
96 1°7616 2 
97 2°7136 t 
98 0:2148 I 




















There is no evidence of interaction within families 95, 96, 97 and 
98 but we cannot neglect the possibility of interaction in families 93 
and 94. Further evidence must be sought in future years. 

That the difference between F, families is a result of segregation 
of genes controlling disjunction is confirmed by the variation between 
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family pairs in the F,. Fig. 3 shows that the disjunction values of 
parents (F;) and offspring (F,) are closely correlated. It is clear, 
therefore, that the frequency of chromosome disjunction as opposed 
to non-disjunction is subject to genotypic control. 
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DISJUNCTION FREQUENCY : F3 FAMILIES 


Fic. 3.—Graph showing correlation of disjunction frequencies between the three F, 
families and their F, offspring. Data taken from tables 2 and 4, the figures used in 
the F, being the estimates from the pooled pairs of families. 


4. DISCUSSION 


The above account shows that variation in disjunction frequency 
of the same chromosome interchanges in rye is due in part to variation 
in the genotype. Heritable variation in disjunction within a species 
may well also account for the variation in pollen and ovule sterility 
reported by Burnham (1934) in maize “cultures” heterozygous for 
interchange T8-9. 

There is evidence suggesting that disjunction in the permanent 
interchange hybrids of QOenothera is genotypically controlled. The 
species have been shown to have characteristic disjunction frequencies. 
Cleland (1926), for example, reports that just over 20 per cent. of 
metaphases in O. muricata show deviation from complete disjunction, 
while in O. pratinicola (Kulkarni, 1929) only 3 to 7 per cent. of meta- 
phases show such deviations. 

Heritable variation in chromosome behaviour is subject to the 
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action of natural selection in just the same way as the more con- 
ventional morphological characters (see Rees, 1955). In interchange 
heterozygotes natural selection would clearly favour those genotypes 
which lead to a high proportion of disjunctional arrangements and 
to a relatively high fertility. Such selection may well have played 
an important part in populations of two kinds, firstly, where inter- 
changes float in the population, some individuals being homozygotes 
and some heterozygotes (e.g. Campanula, Darlington and Gairdner, 
1937), and secondly, where the interchanges are fixed in the hetero- 
zygous condition as in Oenothera and Rhoeo. In both these types of 
population the heterotic advantages of structural heterozygosity are 
counterbalanced by a loss of fertility. Selection may well have 
minimised this infertility by adjustment of the degree of disjunction at 
meiosis. 


5. SUMMARY 


1. The frequency of genetical disjunction of interchange con- 
figurations in rye is shown to be genotypically controlled. 

2. There is some evidence that in plants of some families the 
disjunction properties of two interchange configurations are not 
independent. 


3. Since the disjunctional properties of interchange heterozygotes 
are heritable, selection would be effective in adjusting the degree of 
disjunction to achieve greater fertility. Such selection may well 
have played an important part in the establishment of interchange 
heterozygotes in populations of species in Oenothera and Campanula. 


APPENDIX 


The data from each family in F, and F, is a combination of three groups. Two 
of the groups stem from the plants heterozygous for the A and the B interchanges 
respectively, and each of these provides two classes of observation, viz. disjunction 
and non-disjunction of the interchange in question. Thus each provides one degree 
of freedom for testing agreement with any expected frequency of disjunction. The 
third group, consisting of the double heterozygous individuals, provides three 
classes of observation, viz. disjunctional for both, for only one, and for neither of 
the interchanges. This third group gives us two degrees of freedom. 

There are thus four degrees of freedom appropriate to each row of tables 1 and 3. 
By using the values of /, and /, estimated from any one row and finding expected 
values for that row, we get a x* for two degrees of freedom, the other two having 
been used in estimating p, and p,. Similarly, by using estimates derived from the 
F, totals and finding expectations for the three F, families we get a x? for 3x4 
degrees of freedom, less the two degrees of freedom used in estimating the parameters, 
i.e. a x* for ten degrees of freedom. 

Since there are no observations on plants heterozygous for interchange B in 


families 97 and 98 of the F,, the number of degrees of freedom appropriate to these 
two families is (2+-1)—2 = 1. 
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|. INTRODUCTION 


Vittax and Mota (1953) treated a Triticum-Secale F, hybrid with 
colchicine solution and found in the progeny plants with 42 chromo- 
somes instead of 56 as expected. These plants resembled the wheat 
parent and the “ knob ” appendages characteristic of rye chromosomes 
at pachytene were missing. ‘They therefore believed that the rye 
chromosomes had been lost at the same time that the wheat chromo- 
somes had been doubled through the action of colchicine. 

During the summer of 1953 an F, plant with the somatic chromo- 
some number of 2n = 21, derived from a cross between Avena barbata 
(2n = 28) and Avena strigosa ssp. hirtula (2n = 14), was injected with 
O-I per cent. aqueous colchicine solution at a stage when meiosis 
was imminent in the pollen mother cells. The injection was made 
into the elongating flowering culms. The aim was to repeat the 
process of differential genome elimination which had been postulated 
by Villax and Mota (loc. cit.). 

A study of meivsis in the Avena hybrid was undertaken on pollen 
mother cells stained in aceto-carmine after fixation and storage in 
Carnoy (Darlington and La Cour, 1947) with added ferric chloride. 


2. CYTOLOGICAL OBSERVATIONS 


A study of meiosis revealed that the effect reported for the Triticum- 
Secale hybrid had not been reproduced, and that the great majority 
of pollen mother cells showed a chromosome complement and meiotic 
behaviour normal for this hybrid. 

Another abnormality was found, however, in the occurrence of a 
small proportion of binucleate pollen mother cells. In order to 
ascertain if this aberration was due to the colchicine treatment, an 
untreated hybrid of the same parentage was studied. This plant 
showed the same abnormality, and the data presented in this paper 
have been derived from it. 


(i) Cells 
Of 27 anthers examined, 23 were found to contain binucleate 
pollen mother cells, ranging in number from 1 to g per anther, and 
with a mean frequency of 1-2 per cent. of the total pollen mother cells 
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examined (table 1). In the preparations these cells were frequently 
found to be associated in a group, indicating an origin in the same 
TABLE 1 


Frequencies of binucleate P.M.Cs. in anthers which contain them compared with 
Srequencies of uninucleate P.M.Cs. at various stages of meiosis 





















































Stage of division Total cells 
Slide Cell Early ego « 
a. type prophase Diakinesis} MI Al Int. | MII | AII | TII UN BN 
18 UN* 224 ae ee <2 Pe ne and re 224 ses 
BN (C)tf 4 ~ és we ine ve oes er 4 
19 UN 247 das ane re ae ane nea dea 247 ose 
BN (C) 8 ies aE mc A es 8 gs ee ee 8 
21 UN 135 nes ous =a pe or a dus 135 ove | 
BN (C) 3 ins ‘iy A ee 1 aa. ee ee 3 | 
24 UN er 257 150 acs iia ads Res Raa 407 per 
BN (C) ove 3 2 oes soe | eee eee Hee 5 | 
17 UN are 282 55 1 22 es ite ae | 37° “as 
BN (C) 9 ies Fe a a aes “ii | 9 
I UN es Per 448 2 ans cs ce wl 45° wel 
BN (C) 2 ore 3 ade ise Ae rie Fe 5 | 
| | 
6 | UN BY okay = 514 4 13 eas er | ae 
BNC) | sk 255 I des es ne aes és3 | 1 
22 UN ~ a 470 tis ee a ves ce 1 ee aid 
BN(C) | 3 ve I Sc a | I 
| 
26 | UN | | 459 I | 46c 
BN (C) 7 as 7 
25 UN | | 243 8 8 | | 259 od 
BN (C)_ | | 2 I 1 | 4 
| | | | 
27 UN ae | xe 196 23 Ct sn as ve 225 aed 
BN (C) xa i Be Sasa eas | ae rag ete 3 
| | 
23 UN se a eG 198 16 ‘yn ree ah vo) id 
BN (C) a ee 2 is; Sealy” tka = ae 2 
| | H 
14 | UN or bare 39 4 297 | 55 re acc he nA 
| BN (C) an, ao 3 cD aot Ge ata ae 4 
16 | UN | rae, bes = ce 182 | 12 3 | 372 569 a 
| BN (C) ae | ibe I I a eres a ade 4 
| | | | | | 076 | 60 
* UN = Uninucleate. + BN (C) = Central nucleus of binucleate cell. 
Note.—All peripheral nuclei in all slides are in early prophase. 
part of the anther loculus. The binucleate cells were frequently 
somewhat retarded in their stage of division in relation to the normal 
cells of the same anther, and their lagging was more marked as the 
normal cells approached completion of meiosis (table 1). Onset of 
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prophase appeared to be synchronous in the two types of cell. No 
cells with more than one nucleus were seen at the second division. 


(ii) Nuclei 
In the binucleate cells the nuclei were disposed, one in a central 
or sub-central position, the other in a sub-central or peripheral 
position (plate, figs. 1-4; and fig. 1). 
Both nuclei appear to commence prophase together (plate, fig. 1) 
but a lack of synchronisation soon becomes apparent. As the central 
nucleus proceeds through its stages, the lagging of the peripheral 





Fic. 1.—Binucleate pollen mother cell. Note lack of synchronisation in the division of 
the two nuclei, and the general condition of their chromosomes. Central nucleus at 
metaphase 13;; 6;, 6,; peripheral nucleus existing as long unspiralised prophase 


threads. x 800. 
nucleus becomes very marked (plate, figs. 2-4 and table 1) and it 
does not progress beyond an early prophase condition (plate, fig. 4). 

The central nucleus was seen in favourable preparations of meta- 
phase to be triploid (2n = 21) (fig. 1), indicating that the peripheral 
nucleus was not merely a part removed from the normal complement. 
In the majority of cells, however, chromosomes of the central nucleus 
were clumped together (plate, fig. 2) and frequently had a “ sticky ” 
appearance. There was no evidence of failure of separation at 
anaphase and apart from these two disturbances to metaphase, the 
central nuclei behaved and completed their first division in a normal 
manner. 

The peripheral nuclei, on the other hand, showed a highly 
pathological appearance. Plate figs. 1-4, show some of them in 
association with central nuclei at various stages of the first meiotic 








112 J. W. H. HOLDEN AND M. MOTA 


division. At early prophase (plate, fig. 1) the peripheral possesses 
no recognisable nucleolus, while that in the central nucleus is well 
developed. They are also characterised by reduced staining capacity, 
poor spiralisation and fragmentation of their chromosomes. The 
degree of expression of these latter abnormalities varies from cell to 
cell, ranging from long, extremely thin, poorly stained threads (plate, 
fig. 2) to nuclei where the chromosomes showed contraction, 
improved staining and extreme fragmentation. In occasional nuclei 
these contracted fragments were seen to be paired, and while this 
may be attributed to the zygotene pairing of homologous chromosome 
fragments, it is more probably the result of reproduction. Occasional 
small ring formations are interpreted as the result of reunion of 
fragments. The determination of the number of chromosomes in 
these nuclei was impossible. It is clear that the fragments will 
disappear into the cytoplasm of the pollen mother cell or its daughters. 
One cell was seen where the central daughter nuclei had reached 
first telophase and where the peripheral fragments were. being divided 
into two groups by the new cell wall. 


3. DISCUSSION 


There are a number of reports in the literature of the occurrence 
of binucleate pollen mother cells. These reports are classified in 
table 2 according to the phenotypic similarity of the two nuclei and 
to their genotypic relationship so far as this can be determined. 

In only two cases is there any parallel to the marked phenotypic 
differences found between the nuclei in our material. 

Rees (1955) found binucleate pollen mother cells in inbred rye. 
There, the extra nuclei vary in their chromosome number from one 
acentric fragment to a complete diploid complement. He attributes 
their presence to two different causes, both operating at pre-meiotic 
mitosis, (i) the micronuclei to spindle errors and (ii) the complete 
nuclei to failure of cell wall formation. 

Kihara and Lilienfeld (1934) report binucleate pollen mother 
cells in a J. aegilopoides - Ae. squarrosa hybrid. They attribute the 
binucleate condition to the migration of metaphase chromosomes or 
of whole or parts of nuclei when in meiotic prophase. They find a 
variation in the behaviour of the “ migrant” nuclei. Those existing 
as long poorly staining threads, or bundles of threads, eventually 
degenerate, and are resorbed into the cytoplasm. These nuclei they 
suppose have migrated while in the prophase condition. In others, 
the “migrant”? nucleus exists as normally spiralised metaphase 
chromosomes which are thought to have migrated in that condition, 
or as entire prophase nuclei which have come to lie in contact with 
the legitimate nucleus, and have subsequently developed in a normal 
manner. No indication of chromosomal or nuclear migration has 
been seen in our material, and another explanation of the origin of 
the binucleate condition will be advanced later. However, this 
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suggestion of a relationship between nuclear behaviour and position 
in the cell is similar to that found by us. 

A notable feature of the peripheral nuclei is their lack of a nucleolus. 
The point arises whether they are complete nuclei, or micronuclei 
deficient in nucleolar organisers. In the latter event their lack of 
nucleoli and subsequent degeneration could be attributed to their 
genetic deficiency. This explanation of their differentiation must be 


TABLE 2 


Classification of reported cases of binucleate pollen mother cells according to 
genetic similarity and phenotypic behaviour of the nuclei 








GENOTYPE 
Presumed identical | Known to be different 
a —_— — - - - _— ' 
l 
Lactuca. Gates and Rees, 1921. Triticum-Aegilops. Kihara and 
<ea. Randolph and McClintock, 1926. | _ Lilienfeld, 1934. 
Raphanus-Brassica. Karpechenko, 1927. | Pyrus. Crane and Thomas, 1939. 
Iris. Inaryama, 1929 
Triticum-Aegilops. Kagawa, 1929. 
Xea. Beadle, 1930. | 
| Prunus. Darlington, 1930. 
B | ca | Crepis. Hollingshead, 1930. 
> & | Triticum-Aegilops. Kihara and Lilienfeld, | 
1 1934. 
a Lolium. Jenkin and Thomas, 1939. (1, 2, | 
* 
- . 
Secale. Miintzing and Prakken, 1941. 
Triticum-Secale. Nakajima, 1954. (1) 
Allium. Darlington and Haque, 1955. 
Secale. Rees, 1955. (2) 
as 
BE | Triticum-Aegilops. Kihara and Lilienfeld, | Secale. Rees, 1955. 
#5 | | 1934. | 
S25 | Avena. H. and M. 
a? 

















* 1. Imperfect synchronisation. 
2. Reduced spiralisation of one chromosome complement. 
3. Some asynapsis in one nucleus. 


rejected, however, since all the peripheral nuclei lack a nucleolus 
and therefore it would entail the operation at pre-meiotic mitosis of 
a mechanism for the production of regular genetic differences between 
sister nuclei. 

Since it has been possible in a few cells to determine that the 
central nucleus consists of the full triploid complement, we make the 
assumption that the peripheral nucleus is present as a result of failure 
of cell wall formation at pre-meiotic mitosis. 

When two nuclei are present in a common cytoplasm it is possible 
that competition would occur between them during prophase for 
the materials necessary for chromosome synthesis. If the proposed 
mode of origin of the binucleate cells is correct, competitive advantage 
(in a genetic sense) cannot be postulated as the prime cause of the 
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difference in behaviour of the two nuclei, for they are genetically 
identical. The prime cause must therefore lie outside the nucleus 
and in the cytoplasm. 

This view is supported by the second notable feature of these cells, 
that is, the consistent relationship between position and behaviour 
of the nuclei. As already described, the degenerate nuclei are farther 
from the cell centre than the normal nuclei. In fact they are usually 
peripheral. 

Koller (1947) referring to the synchronisation of meiosis or mitosis 
in binucleate or multi-nucleate cells makes the point that the 
simultaneous division of several independent nuclei can only be 
understood by assuming a distribution throughout the cell of the 
substance, or its precursors, necessary for chromosome synthesis. He 
also suggests that the demands of the individual nuclei guarantee a 
stable equilibrium in the even distribution of the product of the 
“* centres of synthesis”. While this may well be true, it cannot apply 
in all cases. 

Cytoplasmic differentiation has been shown by La Cour (1949) 
to be closely associated with nuclear differentiation in form and 
function in pollen grains of Scilla sibirica and Luzula purpurea. ‘The 
behaviour of the nuclei in the Scilla tetrad, with the generative nucle1 
arranged at the apices of the inner walls, and the vegetative nuclei 
in a peripheral position, bears a close resemblance to the orientation 
of the nuclei in our material (see fig. 2). The Luzula pollen mother 
cell (La Cour, oc. cit.) illustrates an intermediate developmental stage 
between the binucleate Avena pollen mother cell and the differentiated 
tetrad as in Scilla, and shows with its differentiated nuclei but undivided 
cytoplasm an even more remarkable similarity to the phenomenon 
described here. In Luzula the four deeply staining generative nuclei 
lie in a central position, with the vegetative nuclei towards the 
periphery. Darlington and Mather (1949) infer from the Scilla tetrad 
that the concentration of “ generative-nucleus-forming substances ” 
is central in the pollen mother cell. The differences described here 
between the two nuclei in the Avena pollen mother cell support this 
view, since the function of the generative nucleus on the one hand, 
and the central pollen mother cell nucleus on the other, is to divide 
again, and, whatever factors are necessary for the division of the 
former, are likely to be necessary for the division of the latter also. 

It is suggested, therefore, that a cytoplasmic gradient is responsible 
for the difference between the two nuclei of the pollen mother cell, 
just as it is responsible for the differentiation of the pollen grain three 
cell-generations later. While the gradient in the two types of cell 
may be in respect of different factors, both operating to produce a 
similar effect on nuclear behaviour, it is also possible that the factors 
may be identical. Mather (1948) has shown how the nuclear genotype 
may alter the nuclear phenotype through the medium of the cytoplasm 
and has pointed out the lag in time that frequently occurs between 
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the presence of a nucleus in the cytoplasm and its effect on itself or a 
subsequent nucleus in a later cell generation. In the case of the 
binucleate cells described, we may have an example of an exceptionally 
early effect of the normal pollen grain differentiation process, brought 
on by the abnormal binucleate condition. 

Fig. 2 (d) shows the spatial relationships of the differentiated areas 
of cytoplasm in the tetrad of Scilla (after La Cour, 1949) and their 
probable relationship to those postulated for the pollen mother cell, 2 (a). 














Fic. 2.—Relationship between cytoplasmic differentiation in pollen grain (d) (adapted 
—_ La Cour, 1949) and that postulated for tetrad (c), dyad (b) and pollen mother 
cell (a). 


Degenerate nuclei : broken line. 

Generative nuclei : black. 

Concentration of “* generative-nucleus-forming substance ” : stippled. 

Arrows in (¢) indicate direction of gradient and direction of migration of tetrad nuclei. 


According to Sax and Husted (1936) the primary microspore 
nucleus at the time of its division lies near the inner wall, while they 
and Geitler (1935) have shown that in a large majority of species the 
generative nucleus also lies in this position. The fact that the primary 
nucleus migrates to this position and organises its chromosomes and 
spindle eccentrically within the cell is strong indication that the 
differentiation gradient retains the same orientation as that in the 
pollen mother cell, despite two intervening nuclear and cell divisions. 

There remains the apparent inconsistency of the dyad nuclei, 
which divide in a position remote from the focus of the gradient, 
where they might be expected to suffer the same degeneration as 
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the vegetative nucleus of the pollen grain and the peripheral nucleus 
of the pollen mother cell. However, two factors are worth considera- 
tion in this connection, namely (1) the haploid nucleus of the dyad 
consists of chromosomes which have already divided and which 
bring into the nucleus their own charge of chromatin material, and 
(ii) the interphase stage of this nucleus is usually very short and its 
need for physiologic activity, apart from the reproduction of the 
centromeres, would therefore seem to be nil. Thus as a result of its 
special internal self sufficiency it is independent of the cytoplasmic 
gradient. 

Since we have no reason to suppose that the peripheral nucleus 
lacks a nucleolar organiser, the absence of a nucleolus from this nucleus 
is revealing, for the nucleolar material is believed to be intimately 
associated with chromosome synthesis. It is possible that the forces 
operating to maintain the gradient with its focus at the cell centre 
have also been responsible for the transference of nucleolar material 
from the peripheral nucleus at the preceding telophase. — 

The differentiation of the nuclei in the binucleate pollen mother 
cells would therefore seem to be a premature function of the 
phenomenon which is later to differentiate the nuclei in the pollen 
grain. 


4. SUMMARY 


1. A small proportion of binucleate pollen mother cells was found 
in the F, hybrid Avena barbata x A. strigosa ssp. hirtula, 2n = 21. 

2. A marked difference in the behaviour of the two nuclei at 
first meiotic division was found to be associated with their positions 
in the cell. 

3. The central nucleus was essentially normal in behaviour ; it 
was triploid. 

4. The peripheral nucleus contained no nucleolus. Its chromo- 
somes showed reduced staining capacity, poor spiralisation, and 
fragmentation. It degenerated and disappeared. 

5. The binucleate condition is attributed to failure of cell wall 
formation at pre-meiotic mitosis and the two nuclei are therefore 
considered to be identical sisters. 

6. Cytoplasmic gradients are supposed to determine differentiation 
of the genetically identical nuclei in the pollen grain. It is proposed 
that the same cause is responsible for the differentiation of the nuclei 
in these binucleate cells, for there is a marked similarity in their 
behaviour to that of the nuclei in the pollen grain. 
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Plate 


Illustrating differentiation of nuclei in binucleate pollen mother cells of the hybrid 


A, barbata X A, strigosa ssp. hirtula F,. Aceto-carmine. 


Fic, 


Fic. 


1.—Early prophase. Nucleolus present in central nucleus and lacking in peripheral. 
x 1100, 

2.—Central nucleus at diakinesis, configurations clumped. Peripheral at early prophase 
and showing fragmentation. 1300. 


Fic. 3.—Central nucleus at anaphase-telophase. Univalents dividing on plate. Extreme 
fragmentation of peripheral. x 1510. 

Fic. 4.—Central nucleus at telophase. Peripheral still in an early prophase condition and 
showing extreme fragmentation. x 880. 
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THE GENUS NICOTIANA. Origins, Relationship and Evolution of its Species in the 
Light of their Distribution, Morphology and Cytogenetics. By T. H. Goodspeed. 
Published by the Chronica Botanica Co., Waltham, Mass. (Wm. Dawson and Sons, 
Ltd., London, W.C. 2) as Volume 16 of Chronica Botanica. Buckram, pp. xxii+536 
118 plates and illustrations. $12.50. 

This work contains a detailed survey of the geographical distribution, 
morphology and cytology (including chromosome morphology) of the 
species of the genus WVicotiana. The evolutionary relationships of the various 
species are discussed in the light of these facts and in the light of cytological 
studies of interspecific hybrids. The work brings within a single volume an 
enormous mass of information, enabling the many ramifications of the 
genus to be seen as a whole. 

Under the joint authorship of T. H. Goodspeed, H. M. Wheeler and 
P. C. Hutchison, a complete taxonomic account is given of all the 60 species 
recognised, with whole-page line drawings to illustrate their morphology. 
Although the book will at once become the standard work for the taxonomist, 
it will at the same time be of value to the cytologist and all interested in 
evolution, particularly in the mechanisms of species formation. In this 
respect, two sections of the genus show aneuploid series in their chromosome 
numbers, one at the diploid and one at the tetraploid level. All the other 
sections show great uniformity of chromosome number, and in them new 
species seem to have arisen, some by allopolyploidy and some by gradual 
evolutionary divergence without change of chromosome number. The 
parentage of several of the amphidiploid species, deduced from morpho- 
logical and geographical evidence and from the evidence of chromosome 
numbers and behaviour, is confirmed by details of chromosome morphology. 
Thus the amphidiploid Nicotiana tabacum has a haploid chromosome 
complement made up of 10 with median, 5 with sub-median and g with 
sub-terminal centromeres. The corresponding figures for its diploid 
progenitors, NV. otophora and WN. sylvestris are 7, 0, 5 and 3, 5, 4 respectively. 
An interesting feature of WV. otophora is the marked difference in size between 
the chromosomes with median centromeres and those with sub-terminal 
centromeres, compared with some related forms in which chromosome 
size is more uniform. The diagram (fig. 57, p. 310) which sums up the 
postulated evolutionary relationships of the various species is misleading 
with regard to the parentage of V. tabacum. Many taxonomists will welcome 
the broad concept of a species implied by the statement (p. 145) that 
“these tetraploid races” are “ insufficiently distinct from NV. suaveolens 
morphologically to deserve specific rank ”’. 

Criticism of such an outstanding publication seems out-of-place, but 
it is disappointing to find no reference to the classic work of East and his 
collaborators in determining the genetic basis of incompatibility of pollen 
and style in Nicotiana. Indeed there is hardly any discussion of incom- 
patibility in Nicotiana, apart from a general statement (p. 317) that some 
species are self-sterile and some self-fertile. This omission seems remarkable, 
in view of the undoubted evolutionary significance of self-sterility, and the 
fact that the genetic basis of the most common type of incompatibility 
was first discovered in this genus. Nevertheless the book is a monumental 
work. H. L. K. WuirEeHouse. 
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GENETICS AND METABOLISM. By Robert P. Wagner and H. K. Mitchell. John Wiley, 

New York, and Chapman and Hall, London. 1955. Pp. vii+-444. 60s. 

Many research workers will find this book indispensable. It is however 
in some respects uneven. There might be more attention to non-American 
work. There might also be more care in stating what has been discovered 
rather than in referring the reader (in the manner so common among our 
periodical abstractors and summarisers) to a bibliographic reference. And 
there might be a little more consideration of plants, especially in regard 
to development and to the cell. Such amendments would make the book 
appeal to a wider circle of students. C. D. Dar.incTon. 


THE DISTRIBUTION AND ABUNDANCE OF ANIMALS. By H. G. Andrewartha and 

L. C. Birch. University of Chicago Press. Pp. xv + 782. £5, 12s. 6d. 

This is the best book on the subject that has so far been produced, 
and an outstanding contribution to ecology. It is up-to-date, original and 
extremely well documented: the bibliography, of over 1100 references, 
is a compilation of great value. It is particularly unfortunate, therefore, 
that so important a work should have been published in America and, 
therefore, at a price which virtually places it beyond the reach of private 
purchasers in other countries. It will be widely used for many years to 
come, but most biologists will be driven to consult it in libraries. 

There is really no justification for this excessive price. The book is 
certainly very long, I think unnecessarily so, and it contains numerous 
line-drawings, graphs and tables, but only six photographs. These occupy 
less than half a page each, and they all show views in the interior of 
Australia. The lack of balance which this suggests is fortunately not 
apparent in the text. Both the authors work in Australia but, while good 
use is made of the wonderful opportunities for ecological study which 
that region provides, their world-wide choice of examples is commendable. 
It is, indeed, not clear why these six photographs have been chosen since, 
though instructive in themselves, hundreds of other environments are 
mentioned which merit illustration as much or more. It would have been 
better to omit these pictures if they are to any degree responsible for the 
excessive cost of the work. 

The chief purpose of this book is to explain, as far as possible, the 
distribution and numbers of animals in nature. With that end in view, 
it analyses the environment, and the reaction of animals to it, in great 
detail. In the course of that task, it disposes of several widely held concepts 
which have long confused ecological discussions. For instance, it challenges 
the importance of “‘ competition” in regulating numbers, along the lines 
suggested by Nicholson and many others. Andrewartha and Birch rightly 
point out that animals in nature are remarkably rare relative to their 
requirements, and that it is not usual for them to use up their food, possible 
places for shelter or other requirements and, therefore, that competition 
for them does not normally control population-size (p. 23). They further 
stress, and rightly so it seems to me, that the widely accepted sub-division 
of the environment into “‘ density dependent ” and “ density independent ” 
components is misleading, because the latter type does not exist: there 
appears to be no evidence that the density of the population is independent 
of any aspect of the enviroment (p. 17). 

Too much importance seems to be attached to the chance loss of genes 
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from small populations, a concept elaborated in particular by Sewall 
Wright. It is, indeed, pointed out that to be effective the numbers must 
be very small, so much so that the survival of the group, if isolated, is 
threatened, but that limitation is not observed in some of the instances 
quoted here. For example, it is stressed that the smallness of local popula- 
tions “ provides a possibility of diversity through both selection and the 
chance loss of genes ” (p. 716). The moth Rhyacia alpicola, cited by Huxley 
from this point of view, is still suggested as a possible example of these two 
effects ; though its local populations, certainly comprising many thousands, 
are vastly too big for such random loss of genes to contribute to the observed 
differences between them. 

It is unfortunate that no strict definition of polymorphism is given 
and that the term is used so widely that the special features of that condition, 
which is of great evolutionary importance when strictly applied, are lost. 
Indeed the statement (p. 671) that “ all Mendelian populations are more or 
less polymorphic ” appears to be meaningless. 

The subject of fluctuations in numbers is well handled, so is their effect 
of increasing the rate of evolution. It is not made clear, however, that 
organisms are capable of generating, autonomously, as it were, their own 
cycle of numerical fluctuations. 

The very large number of quotations and statements from other authors 
seem to be rendered with commendable accuracy. There are a few failures 
in this respect, but they are not to be wondered at in a work of such 
magnitude. Thus on p. 715 the two populations of the moth Thera 
juniperata in Britain, isolated by absence of the food plant, are mentioned 
and it is said that “‘ each population is now recognised as a distinct species 
(Huxley, 1942, p. 185) ”. This is not so, nor does Huxley say it : his words 
are “ the British forms of this species are restricted to two areas, one in the 
north, the other in the south, and as a result subspeciation has occurred ”’. 

The extensive mass of information in this book never gets out of control, 
and the vast array of examples seems always quoted to good purpose. The 
short section on the Genetic Aspects of Ecology is the least satisfying, and 
it provides but an incomplete survey of that aspect of the subject. But 
it is not as a genetic work that this book is to be judged, useful as it will 
prove to geneticists. It has succeeded in the task of analysing the environ- 
mental factors which govern the distribution and numbers of animals 
more thoroughly than has any other work, and it is, therefore, an important 
contribution to biology. E. B. Forp. 
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ABSTRACTS of Papers read at the HUNDRED AND NINETEENTH 
MEETING of the Society, held on 18th and 19th NOVEMBER 1955, at 
the GALTON LABORATORY, University College, Gower Street, 
London 


EVIDENCE FOR AN EXCHANGE HYPOTHESIS FOR ‘‘ CHROMATID ”’ 
ABERRATIONS INDUCED BY X-RAYS 
S. H. REVELL 
Chester Beatty Research Institute, Fulham Road, London 
A recent paper (S. H. Revell, Proceedings of the Radiobiology Symposium, Liege, 1954, 
Butterworth) suggested a new interpretation for ‘‘ chromatid ” changes in Vicia 
faba root tip cells: that all aberrations of this category are actually the result of 
chromatid exchanges, structurally similar at metaphase to meiotic exchanges. 
Qualitative evidence only was then given for the hypothesis. 

In the present paper quantitative data from X-ray experiments with Vicia 
are offered in support of the same hypothesis. It is shown that there is no early 
excess of true chromatid breaks, such as would be difficult to accommodate within 
the new hypothesis; and that certain types of aberration are observed in the 
proportions which the hypothesis predicts. 

Some implications of this interpretation are discussed. 


DIFFERENTIATION OF CELLS WITHIN CLONES OF SALMONELLA : 
SPONTANEOUS PHENOTYPIC VARIATION IN RESPECT OF 
MOTILITY 

C. QUADLING and B. A. D. STOCKER 
The Lister Institute, London 

Morphologically, bacteria of the genus Salmonella consist of short rods which 
may be motile by means of flagella (wild type), or non-motile, lacking flagella. 
Mutation occurs from the motile to the non-motile condition (and vice versa). 
In some strains only a characteristic proportion of the individuals are motile. It 
has now been shown that motile cells occur at very low frequencies in some strains 
formerly considered non-motile (circa 10-5 in one case). These rare motile cells 
(isolated by micro-manipulation) are not mutant; their progeny after several 
generations are all or nearly all non-motile, and indistinguishable from the parent 
strain. 

The distribution of such motile cells amongst replicate cultures (studied by a 
pour-plate technique in semi-solid agar) suggests that small groups of motile cells 
arise by division of single differentiated cells. The incidence of “ differentiation 
events ” causing the appearance of such small groups appears to be random amongst 
the population. It seems that each cell of a strain behaving in this manner has a 
genetically determined low probability of undergoing such an event. The event 
probably consists of a transient ability to synthesise new flagella. The “‘ penetrance ” 
of a character (motility) can thus be studied amongst clonal cell populations. 


A METHOD FOR THE DETECTION OF AUTOSOMAL LETHALS 
INDUCED BY HIGH FREQUENCY RADIATION IN MICE 
J. B. S. HALDANE 
Dept. of Biometry, University College, London 
A mouse from an irradiated multiple recessive stock is crossed with a multiple 
dominant, and lethals are detected by shortages of recessives in F,. An F, of n 
mice is searched for a lethal if it contains less than r recessives at any locus. r is 
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chosen so that the probability of suspecting a lethal which is not present is less 
than one in 300. The length of chromosomal map swept per locus is about 
set®) cM. With n = 30, r = 1, this is about 16 cM, but each such length is 
only searched in } of the F,. Thus with 8 markers one could detect about 2 per cent. 
of autosomal lethals by scoring 30 mice. This method should detect lethals with 
about the same frequency, for a given dose of radiation, as Russell detected visibles 
by irradiating dominants at 7 loci, and crossing to recessives. 


RESOLVING POWER OF GENETIC ANALYSIS IN TERMS OF DNA 


G. PONTECORVO and J. A. ROPER 
Dept. of Genetics, University of Glasgow 


The resolving power of recombination in Aspergillus nidulans (total map 500 units) 
is now down to 10~* and in Drosophila (total map 280 units) to 10-5. The DNA 
per haploid nucleus is 0-044 picograms in Aspergillus and 0-085 in Drosophila, 
corresponding to 4107 and 7:7X107 nucleotide pairs. The minimum recom- 
bination fraction measured in Aspergillus—1 x 10~* between two alleles of the pro 1 
region (Forbes)—corresponds to 8 nucleotide pairs or 27A, and in Drosophila— 
8x 10-® between two alleles of the w region (Mackendrick)—corresponds to 
240 nucleotide pairs. In Aspergillus the chromonema—if it is DNA—could well 
be made up of only one Watson-Crick duplex fibre and cannot be made up of 
more than about 8. The ratio of the smallest recombination fraction within one 
region in Aspergillus (ad 8) to the total recombination for that region is 1 : 150. 
In the w region of Drosophila it is 1 : 62. This supports the idea that the mutational 
sites separable by crossing over within each region are in the order of hundreds. 
(By “‘ region ” we mean the series of mutational sites at which mutant alleles behave 
as allelic to one another.) These results agree with those of Benzer with phage T4. 


SIGNIFICANCE TESTS FOR HERITABILITIES AND ANALOGOUS 
PARAMETERS ESTIMATED FROM VARIANCE COMPONENTS 
B. WOOLF 
Institute of Animal Genetics, Edinburgh 

In a standard breeding experiment, s sires are each mated to d different dams, 
each of whom has n offspring on whom some trait is measured. Assuming variability 
among the offspring is made up of three additive components, S due to sires, D 
from dams and W among full sibs, the analysis of variance is 
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Derived estimates of heritability are 4S/(S+D+W) from sires, 4D/(S+D-+W) 
from dams and 2(S+D)/(S+D+W) from both parents. There is no accepted 
test of significance for such parameters. A simple computational routine is given 
for the large-sample formula of Osborne and Paterson (1952), and an easily 
evaluated transformation is described which stabilises the variance and should be 
more nearly normal in distribution. The tests can be used when d and n are not 
constant. The methods are quite general, being applicable to analogous parameters 
estimated from three or more mean squares in a “ nested classification.” 
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HOMOSTYLY IN THE PRIMROSE—A PROGRESS REPORT 


JACK L. CROSBY 
Botany Dept., Durham Colleges in the University of Durham 


An account will be given of some of the results which have so far been obtained 
from the author’s research over the past ten years on homostyle populations in 
Somerset, Buckinghamshire, and County Durham. These will include the results 
of studies on incompatibility reaction, on inbreeding depression, and on the genetical 
structure of natural populations, together with detailed analyses of fertilisation in 
natural populations and of spatial variation in homostyle frequency. 


THE EFFECTS OF SELECTION AND LINKAGE ON 
INBREEDING PROGRESS 


E. C. R. REEVE 
Institute of Animal Genetics, Edinburgh 


Selection against homozygosis at a given locus, due to lower survival of the 
homozygotes, will reduce the rate of progress towards fixation during inbreeding 
if there is selection within lines only, or may prevent fixation altogether if selection 
between lines also occurs. A simple general solution will be given for the effects 
of selfing when the homozygous phases of the two alleles at a locus are under different 
selection pressures and the intensities of selection within and between lines are also 
different. This makes possible a comparison of the relative importance of selection 
within and selection between lines in holding up inbreeding progress. 

Selection against homozygosis at one locus will slow down inbreeding progress 
at other loci which are linked to it but not directly affected by selection. A general 
formula for the magnitude of this linkage effect with inbreeding by selfing will be 
given, and some particular cases will be discussed. 

When inbreeding is by sib-mating, this problem presents much greater mathe- 
matical difficulties. Some preliminary results will be briefly discussed. 


GENETIC CORRELATION BETWEEN BODY SIZE AND 
EGG PRODUCTION IN DROSOPHILA MELANOGASTER 


FORBES W. ROBERTSON 
Institute of Animal Genetics, Edinburgh 


Progeny tests and selection experiments in wild stocks of Drosophila melanogaster 
imply a lack of genetic correlation between body size and one of the principal 
components of fitness, namely rate of egg production. However, there is also 
evidence of a phenotypic correlation between the two characters among individuals 
of the population, reared under favourable conditions. This cannot be attributed 
merely to environmental variation since there is no correlation among genetically 
identical flies from a cross between inbred lines, reared under similar conditions. 
Also the phenotypic correlation due to segregation can be greatly increased by 
rearing larve under standard sub-optimal conditions. There is again no correlation 
in the genetically identical flies. Estimates of environmental variance and the 
additive genetic variance leave a considerable fraction of the phenotypic variation 
of flies from the wild stock unaccounted for. It appears that such non-additive 
effects are responsible for the observed conrelation between size and egg production 
and that they are particularly sensitive to nutritional variation for their joint 
expression. The data are considered in relation to the results of selection and 
methods of detecting genetic variation in fitness. 
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THE ABO BLOOD GROUPS AND SECRETOR STATUS IN 
DUODENAL ULCER 


C. A. CLARKE, R. B. McCONNELL 
31 Rodney Street, Liverpool 
and 


P. M. SHEPPARD 
Dept. of Zoology, University Museum, Oxford 


Recent work has suggested an association between blood group O and duodenal 
ulcer in several areas of Europe. The finding has been obtained by comparing 
the blood groups of ulcer patients with a control series of unaffected people living 
in the same area. 

Such controls can be unsatisfactory in that a population of mixed origin may 
contain elements with a high frequency both of group O and duodenal ulcers without 
the two being causally connected. 

Data are presented using family studies which are not subject to this criticism. 
Up to the present (176 sibships) there is no significant increase in the frequency of 
O in the propositi compared with their unaffected sibs. 

The secretor status has also been investigated in these sibships. Here there is 
a significant increase in non-secretors among the duodenal ulcer patients compared 
with controls. The implications of this finding are discussed. ‘ 


FATHERLESS FAMILIES OF LEBISTES RETICULATUS 


H. SPURWAY ; 
Dept. of Biometry, University College, London 


Sixteen related Lebistes reticulatus between the ages of 184 and 769 days of age 
have produced fatherless litters totalling more than 270 babies and embryos. With 
two exceptions (1 male and 1 hermaphrodite) the individuals (about 100) surviving 
to be sexed have been females, and 3 : 1 segregation ratios of three domestic mutants 
have been observed among them. Fatherless females can themselves produce 
fatherless litters. With two exceptions the sired offspring both from fatherless 
females, and from the mothers of fatherless females, occur with a normal sex ratio. 
Females which are the offspring of crosses between two stocks produce larger 
fatherless litters than do females from the original stock, but these large litters 
contain many miscarriages. The one fatherless male was sterile. The herma- 
phrodite produced g fatherless litters, several of which were entirely still-born. 
‘** She ” performed male-courtship behaviour to females but these did not produce 
offspring. After death, well developed testis-tissue was found in “ her.’’ All other 
15 mothers were morphologically perfect females and no testicular tissue or sperm 
were found in the two examined histologically. 


PRIMARY SEX RATIOS OF THE SEA-WEED FLY 
COELOPA FRIGIDA (F.) AND THE HOUSE MOUSE 


U. PHILIP 
Dept. of Zoology 
and 


D. A. EVAN 
Dept. of Statistics, King’s College, Newcastle-on-Tyne 


Zygotic ratios, unmodified by mortality subsequent to fertilisation, are rarely 
preserved in crosses. It is therefore of interest to analyse these ratios statistically 
when they are encountered. The sex ratio of the sea-weed fly Coelopa frigida was 
found, in an extensive experiment, to show hardly any variability between the 
families with low mortality. The chance of getting a more homogeneous sample 
is less than 5 x 10-, clearly not in agreement with the binomial distribution. An 
analysis of the data by McDowell and Lord on the house mouse also shows an 
unexpected homogeneity though not as strikingly as in Coelopa. 
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EFFECT OF GAMMA IRRADIATION ON THE AFFINITIES 
OF LOLIUM PERENNE AND FESTUCA PRATENSIS 


J. D. H. REUSCH 
Dept. of Agricultural Botany, University College of Wales, Aberystwyth 


Attempts to increase the crossability of Lolium perenne and Festuca pratensis by 
irradiating the gametes with gamma rays from a Co® source gave some unexpected 
results, which, however, have been consistent for two consecutive years. 

Using Lolium perenne as female and crossing with pollen from irradiated Festuca 
pratensis, the percentages of seedset and germination were determined in relation 
to irradiation dose. 

The results showed that up to about 500 r there was a decrease in seedset and 
germination. This was followed by a sharp rise to about 1500 r with a gradual 
decrease towards 3000 r. 

The possible significance of these results will be discussed. 


EXPERIMENTAL EVIDENCE CONCERNING ‘‘ AFFINITY ”’ 


M. E. WALLACE 
Dept. of Genetics, Cambridge 


In 1953 Dr Donald Michie and the writer presented the concept of, and observa- 
tions supporting, a new phenomenon involving non-random assortment of 
independent markers, found in the House Mouse; this concept was termed 
“* affinity’ (Nature, 171, 26). An association of a chromosome III marker with 
three chromosome V markers, reported by the writer at that time, has been thoroughly 
investigated in order to test whether it can be explained on the hypothesis of 
“ affinity.” The results are significant and are given with a tentative map of 
linkage group V the position of the point in V responsible for the association. The 
possibility that this is the centromere is briefly discussed. 





NOTICE OF CONGRESS 
FIRST INTERNATIONAL CONGRESS OF HUMAN GENETICS 


Tuis congress will be held in Copenhagen, rst to 6th August 1956. 

Provisional programme and information are sent on request. 
Address: The Secretariat of the Congress, University Institute for 
Human Genetics, 14, Tagensvej, Copenhagen, N., Denmark. 





PROFESSOR @JVIND WINGE 


Professor Winge’s 7oth birthday is due on May 19th, 1956. At the same time 
he completes his long service to the Carlsberg Laboratory. Geneticists all over the 
world will wish to offer their tribute of congratulations on this occasion to one 
whom we have admired for the immense versatility and fruitfulness of his successive 
researches. To no other living man do we owe the elucidation of so many unfore- 
seen novelties, from his early work on polyploidy, sex chromosomes and sex-linkage, 
to his later pregnant researches in the Saccharomycetes. Each topic to which he has 
given his mind has turned out to be of importance for the development of our science. 

&. D. D., BR. A. F. 











